Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,837,903 A→G K138R (AAA→AGA)  b1757 → putative thiosulfate sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,837,9030AG98.5% 439.0 / ‑2.9 138K138R (AAA→AGA) b1757putative thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (70/65);  total (72/66)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 147 total aligned reads displayed.

GACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837815‑1837992
                                                                                        |                                                                                         
gACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTg                                                                                       <  2:2210022/93‑1 (MQ=255)
  cGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTgac                                                                                     >  1:219306/1‑91 (MQ=255)
  cGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGcc                                                                                     <  1:759706/93‑1 (MQ=255)
   gCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCg                                                                                    >  2:4113614/1‑93 (MQ=255)
     cGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCa                                                                                  >  2:1217985/1‑93 (MQ=255)
     cGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCa                                                                                  <  2:1580423/93‑1 (MQ=255)
      gTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAt                                                                                 >  1:1124746/1‑93 (MQ=255)
       tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                >  2:3114000/1‑93 (MQ=255)
       tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                >  1:4737242/1‑93 (MQ=255)
       tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                <  2:4552723/93‑1 (MQ=255)
        cAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAttt                                                                               <  2:4232139/93‑1 (MQ=255)
         aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              >  2:1607918/1‑93 (MQ=255)
         aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              <  1:4673037/93‑1 (MQ=255)
         aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              >  1:371265/1‑93 (MQ=255)
          aaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTg                                                                             >  2:3389535/1‑93 (MQ=255)
          aaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTg                                                                             >  1:3865767/1‑93 (MQ=255)
          aaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCATCTGCAAAGACTGCCGCATTTTg                                                                             >  2:3389536/1‑93 (MQ=255)
           aaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGa                                                                            <  1:3611523/93‑1 (MQ=255)
              gcgACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagca                                                                         >  2:1973455/1‑93 (MQ=255)
                gACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagcagc                                                                       <  2:1152873/93‑1 (MQ=255)
                 aCTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCt                                                                      <  2:684320/93‑1 (MQ=255)
                  cTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTg                                                                     <  2:3700737/93‑1 (MQ=255)
                     cAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGtt                                                                  <  1:2774074/93‑1 (MQ=255)
                     cAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGtt                                                                  <  2:4566975/93‑1 (MQ=255)
                         aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAATACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  2:3470132/93‑1 (MQ=255)
                         aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  2:3329470/93‑1 (MQ=255)
                         aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  2:3329469/93‑1 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAGCCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             <  2:1496330/93‑1 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             >  2:4475384/1‑93 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             >  2:3444779/1‑93 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             <  2:1294059/93‑1 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             >  2:2464688/1‑93 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGCTTATcc                                                             >  2:387428/1‑93 (MQ=255)
                           gCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATccc                                                            >  2:691366/1‑92 (MQ=255)
                            cAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTCATCAGCTGGTTTATCCGc                                                           >  1:860704/1‑93 (MQ=255)
                                tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       >  1:3574473/1‑93 (MQ=255)
                                tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       <  2:3097490/93‑1 (MQ=255)
                                  tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTTCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  1:1880411/93‑1 (MQ=255)
                                  tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  1:2101995/93‑1 (MQ=255)
                                  tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  2:4273880/93‑1 (MQ=255)
                                   aaCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCt                                                    <  1:3343852/93‑1 (MQ=255)
                                    aCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGCCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTg                                                   >  1:55395/1‑93 (MQ=255)
                                     cGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGc                                                  <  1:3195575/93‑1 (MQ=255)
                                      gCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCa                                                 >  1:670776/1‑93 (MQ=255)
                                       cATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCAc                                                <  2:4669001/93‑1 (MQ=255)
                                        atatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACg                                               >  1:1765839/1‑93 (MQ=255)
                                         tatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGa                                              >  1:2001381/1‑93 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCATCTGGTTTATCCGCAATGGCTGCACGAcc                                            <  1:1280176/93‑1 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  1:2487699/1‑93 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  2:2562956/1‑93 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATCGCTGCACGAcc                                            <  2:871068/93‑1 (MQ=255)
                                            ctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCt                                           <  2:2891034/93‑1 (MQ=255)
                                            ctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCt                                           <  2:1899882/93‑1 (MQ=255)
                                             tcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTg                                          >  2:1194795/1‑93 (MQ=255)
                                             tcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTg                                          <  1:3076384/93‑1 (MQ=255)
                                              ccATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGc                                         <  1:2911887/93‑1 (MQ=255)
                                              ccATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGc                                         >  1:450789/1‑93 (MQ=255)
                                                aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:399185/1‑93 (MQ=255)
                                                aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:2855586/1‑93 (MQ=255)
                                                aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:4518907/1‑93 (MQ=255)
                                                   cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    >  2:3942064/1‑93 (MQ=255)
                                                   cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    >  1:2525975/1‑93 (MQ=255)
                                                    tGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAg                                   >  1:4169523/1‑93 (MQ=255)
                                                    tGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAg                                   >  1:1770557/1‑93 (MQ=255)
                                                     gAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAgg                                  >  1:2485609/1‑93 (MQ=255)
                                                     gAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAgg                                  >  2:644122/1‑93 (MQ=255)
                                                      aGTGACGCGCTAAGCGACCCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGt                                 <  2:3857845/93‑1 (MQ=255)
                                                       gTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTa                                <  1:655253/93‑1 (MQ=255)
                                                         gACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTaaa                              >  2:343349/1‑93 (MQ=255)
                                                           cgcgCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAaga                            <  2:860938/93‑1 (MQ=255)
                                                              gcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGAAACAAGGGAAAGAGGt                         >  2:271380/1‑93 (MQ=255)
                                                               cTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGtt                        <  1:2469464/93‑1 (MQ=255)
                                                                tAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       >  2:1080312/1‑93 (MQ=255)
                                                                tAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       >  2:2075205/1‑93 (MQ=255)
                                                                tAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       >  2:4103979/1‑93 (MQ=255)
                                                                 aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      <  1:1213865/93‑1 (MQ=255)
                                                                 aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      <  2:1042288/93‑1 (MQ=255)
                                                                  aGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACg                     <  1:3135208/93‑1 (MQ=255)
                                                                    cGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGc                   <  1:1188207/93‑1 (MQ=255)
                                                                      aaCCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGa                 >  1:2344513/1‑93 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  2:2961034/93‑1 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:2219428/93‑1 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:870963/93‑1 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             >  1:2314153/1‑93 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:3238789/93‑1 (MQ=255)
                                                                             ccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          <  1:3162196/93‑1 (MQ=255)
                                                                             ccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          <  1:4181688/93‑1 (MQ=255)
                                                                             ccGTCGGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          <  2:3508977/93‑1 (MQ=255)
                                                                              cGTCTGCACAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGcc         <  1:3519086/93‑1 (MQ=255)
                                                                              cGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGGTACGGCGAAAACTGcc         >  2:1478219/1‑93 (MQ=255)
                                                                               gTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCg        >  1:339225/1‑93 (MQ=255)
                                                                                tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       >  1:3315674/1‑93 (MQ=255)
                                                                                tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       >  2:44448/1‑93 (MQ=255)
                                                                                tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       <  2:2085097/93‑1 (MQ=255)
                                                                                 cTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGt      <  1:1506961/93‑1 (MQ=255)
                                                                                  tGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTg     <  2:1497450/93‑1 (MQ=255)
                                                                                    cAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGAc   <  1:1288195/93‑1 (MQ=255)
                                                                                     aaagaCTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACt  <  1:3028321/93‑1 (MQ=255)
                                                                                     aaagaCTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACt  >  2:227592/1‑93 (MQ=255)
                                                                                     aaagaCCGCCGCATGTGGAGCAGCGGGGTTATCCGCAAAGGCTGCACGACCCGCAACAAGGGAAAGAGGTTACGGCGAAACCGGCCGGCGAcg  >  2:3616432/1‑92 (MQ=255)
                                                                                        |                                                                                         
GACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837815‑1837992

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACTGGA  >  NC_000913/1837810‑1837995
                                                                                             |                                                                                            
ATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACATTCCCGTCTGCAAAG                                                                                              >  SRR4302544.1715692/1‑94 (MQ=60)
 TGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGA                                                                                             >  SRR4302544.1531866/1‑94 (MQ=60)
  GTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAC                                                                                            <  SRR4302544.2642286/94‑1 (MQ=60)
  GTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAC                                                                                            <  SRR4302544.1686861/94‑1 (MQ=60)
  GGCGACGTCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAC                                                                                            <  SRR4302544.1686860/94‑1 (MQ=60)
      ACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGAC                                                                                        >  SRR4302544.221571/1‑94 (MQ=60)
       CGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCG                                                                                       <  SRR4302544.767116/94‑1 (MQ=60)
          CGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAT                                                                                    >  SRR4302544.1135950/1‑94 (MQ=60)
          CGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAT                                                                                    >  SRR4302544.4157727/1‑94 (MQ=60)
           GTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATT                                                                                   >  SRR4302544.4783940/1‑94 (MQ=60)
             CAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTT                                                                                 >  SRR4302544.375057/1‑94 (MQ=60)
              AAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTG                                                                                >  SRR4302544.3137182/1‑94 (MQ=60)
              AAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTG                                                                                <  SRR4302544.4719224/94‑1 (MQ=60)
              AAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTG                                                                                >  SRR4302544.3904152/1‑94 (MQ=60)
                AACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAG                                                                              <  SRR4302544.3646852/94‑1 (MQ=60)
                   GCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAG                                                                           >  SRR4302544.2275255/1‑94 (MQ=60)
                          CAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTT                                                                    <  SRR4302544.2800755/94‑1 (MQ=60)
                            AAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGATCAGCTGGTTTAT                                                                  >  SRR4302544.3209187/1‑94 (MQ=60)
                              AAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCC                                                                >  SRR4302544.1264017/1‑94 (MQ=60)
                                GCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGA                                                              >  SRR4302544.869219/1‑94 (MQ=60)
                                GCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTCATCAGCTGGTTTATCCGC                                                              >  SRR4302544.869220/1‑94 (MQ=60)
                                  AGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAA                                                            >  SRR4302544.178615/1‑94 (MQ=60)
                                    GTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATG                                                          <  SRR4302544.4594728/94‑1 (MQ=60)
                                    GTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATG                                                          >  SRR4302544.3608752/1‑94 (MQ=60)
                                       TAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTTCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCT                                                       <  SRR4302544.1898592/94‑1 (MQ=60)
                                       TAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCT                                                       <  SRR4302544.2122223/94‑1 (MQ=60)
                                        AACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTG                                                      <  SRR4302544.3375740/94‑1 (MQ=60)
                                        AACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGCCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTG                                                      >  SRR4302544.55952/1‑94 (MQ=60)
                                          CGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCA                                                    >  SRR4302544.677405/1‑94 (MQ=60)
                                          CGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCA                                                    <  SRR4302544.3226028/94‑1 (MQ=60)
                                           GCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCAC                                                   >  SRR4302544.195122/1‑94 (MQ=60)
                                            CATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACG                                                  >  SRR4302544.1783119/1‑94 (MQ=60)
                                             ATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGA                                                 >  SRR4302544.989677/1‑94 (MQ=60)
                                             ATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGA                                                 >  SRR4302544.2020815/1‑94 (MQ=60)
                                               ATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACC                                               >  SRR4302544.2511527/1‑94 (MQ=60)
                                                TCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCATCTGGTTTATCCGCAATGGCTGCACGACCT                                              <  SRR4302544.1292789/94‑1 (MQ=60)
                                                 CTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTG                                             >  SRR4302544.303812/1‑94 (MQ=60)
                                                  TCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGC                                            >  SRR4302544.455261/1‑94 (MQ=60)
                                                  TCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGC                                            <  SRR4302544.3105699/94‑1 (MQ=60)
                                                  TCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGC                                            <  SRR4302544.3415828/94‑1 (MQ=60)
                                                   CCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCA                                           <  SRR4302544.2939898/94‑1 (MQ=60)
                                                    CATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAA                                          >  SRR4302544.403233/1‑94 (MQ=60)
                                                    CATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAA                                          >  SRR4302544.2883165/1‑94 (MQ=60)
                                                    CATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAA                                          >  SRR4302544.4564159/1‑94 (MQ=60)
                                                    CATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTCAGCAGCTGGTTTATCCGCAATCGCTGCACaccctgaca                                          >  SRR4302544.2883163/1‑85 (MQ=60)
                                                       CCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAA                                       >  SRR4302544.2550250/1‑94 (MQ=60)
                                                        CTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAG                                      >  SRR4302544.1377049/1‑94 (MQ=60)
                                                        CTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAG                                      >  SRR4302544.1787885/1‑94 (MQ=60)
                                                        CTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAG                                      >  SRR4302544.4210828/1‑94 (MQ=60)
                                                         TGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGG                                     >  SRR4302544.2509412/1‑94 (MQ=60)
                                                          GAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGT                                    <  SRR4302544.2185791/94‑1 (MQ=60)
                                                            GTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAA                                  <  SRR4302544.661738/94‑1 (MQ=60)
                                                                   GCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTT                           <  SRR4302544.674894/94‑1 (MQ=60)
                                                                    CTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTA                          <  SRR4302544.2493061/94‑1 (MQ=60)
                                                                    CTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTA                          <  SRR4302544.2493060/94‑1 (MQ=60)
                                                                      AAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACG                        <  SRR4302544.1225882/94‑1 (MQ=60)
                                                                       AGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGG                       <  SRR4302544.3164986/94‑1 (MQ=60)
                                                                       AGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGG                       <  SRR4302544.2522632/94‑1 (MQ=60)
                                                                         CGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCG                     <  SRR4302544.1200008/94‑1 (MQ=60)
                                                                          GAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGA                    >  SRR4302544.2366723/1‑94 (MQ=60)
                                                                             CCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGAAAAGAGGTTACGGCGAAAC                 >  SRR4302544.3832458/1‑94 (MQ=60)
                                                                              CTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACC                >  SRR4302544.2336052/1‑94 (MQ=60)
                                                                               TTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCT               <  SRR4302544.2240306/94‑1 (MQ=60)
                                                                               TTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCT               <  SRR4302544.879596/94‑1 (MQ=60)
                                                                               TTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCT               <  SRR4302544.3269787/94‑1 (MQ=60)
                                                                               TTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCT               <  SRR4302544.556430/94‑1 (MQ=60)
                                                                                  CCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCC            <  SRR4302544.4223099/94‑1 (MQ=60)
                                                                                  CCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCC            <  SRR4302544.3192269/94‑1 (MQ=60)
                                                                                   cgacgGCAAAGACTGCGGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCG           <  SRR4302544.4112252/89‑1 (MQ=60)
                                                                                   CGTCTGCACAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCG           <  SRR4302544.3552693/94‑1 (MQ=60)
                                                                                   CGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCG           <  SRR4302544.3552694/94‑1 (MQ=60)
                                                                                   CGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCG           >  SRR4302544.342740/1‑94 (MQ=60)
                                                                                    GTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGG          >  SRR4302544.3347250/1‑94 (MQ=60)
                                                                                      CTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTG        <  SRR4302544.1521900/94‑1 (MQ=60)
                                                                                      CTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTG        >  SRR4302544.3279814/1‑94 (MQ=60)
                                                                                         CAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT     >  SRR4302544.425252/1‑94 (MQ=60)
                                                                                         CAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT     <  SRR4302544.1300880/94‑1 (MQ=60)
                                                                                          AAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACTG    >  SRR4302544.1943581/1‑94 (MQ=60)
                                                                                          AAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACTG    <  SRR4302544.3057336/94‑1 (MQ=60)
                                                                                            AGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACTGGA  <  SRR4302544.703924/94‑1 (MQ=60)
                                                                                             |                                                                                            
ATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACTGGA  >  NC_000913/1837810‑1837995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑