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 I2 R1 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,539,774 G→A A578T (GCC→ACC)  feoB → ferrous iron transport protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,539,7740GA99.6% 841.7 / NA 234A578T (GCC→ACC) feoBferrous iron transport protein B
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (107/125);  total (107/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863
                                                                                          |                                                                                         
 cgtcggtcAGCCGGGTGATCACCCCGGTCTGCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCa                                                                                        <  2:2748105/93‑1 (MQ=255)
  gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAt                                                                                       <  2:5564336/93‑1 (MQ=255)
   tccgtcAGCCGGGTGACCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  1:1657142/93‑1 (MQ=255)
   cccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  2:206909/92‑1 (MQ=255)
       tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGa                                                                                  >  2:168614/1‑93 (MQ=255)
       tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGa                                                                                  >  1:1585414/1‑93 (MQ=255)
         cgCCGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                <  1:2906189/92‑1 (MQ=255)
         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                <  1:3365347/93‑1 (MQ=255)
            cGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             <  1:3923827/93‑1 (MQ=255)
            cGGGGGATCACCCCGGTCTTCAAGACAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             >  1:1430992/1‑93 (MQ=255)
             gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaag                                                                            >  2:6126964/1‑93 (MQ=255)
              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           >  1:3204496/1‑93 (MQ=255)
              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           <  1:3825738/93‑1 (MQ=255)
               gTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTg                                                                          >  2:5004962/1‑93 (MQ=255)
                tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAgggg                                                                         >  2:5846871/1‑93 (MQ=255)
                tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTgg                                                                         >  1:3676322/1‑93 (MQ=255)
                tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTgg                                                                         <  1:4877771/93‑1 (MQ=255)
                 gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGt                                                                        >  1:1971746/1‑93 (MQ=255)
                 gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTAACATGGCGAAAGAAGTGGt                                                                        >  2:1252779/1‑93 (MQ=255)
                  aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTa                                                                       >  2:3004865/1‑93 (MQ=255)
                    cACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGt                                                                     <  1:3871146/93‑1 (MQ=255)
                     aCCCCGGTCTTCAAGCCAATTGGCGTGCATGATGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    <  1:5685038/93‑1 (MQ=255)
                      ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTgg                                                                   >  2:950790/1‑93 (MQ=255)
                        ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGt                                                                 <  1:6133846/93‑1 (MQ=255)
                        ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGt                                                                 >  1:3717854/1‑93 (MQ=255)
                         cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTa                                                                >  1:4623975/1‑93 (MQ=255)
                         cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTa                                                                <  1:4834689/93‑1 (MQ=255)
                           gTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACg                                                              >  1:1125894/1‑93 (MQ=255)
                            tCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGc                                                             >  2:5333960/1‑93 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                                                           >  2:5595481/1‑93 (MQ=255)
                                cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCaa                                                         <  1:1342999/93‑1 (MQ=255)
                                 aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                        <  1:1048761/93‑1 (MQ=255)
                                 aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                        >  2:5115455/1‑93 (MQ=255)
                                  aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       <  2:5923730/93‑1 (MQ=255)
                                  aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       <  2:4839013/93‑1 (MQ=255)
                                  aGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       <  2:623982/93‑1 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      <  2:2295833/93‑1 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      >  1:976627/1‑93 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      >  1:3455170/1‑93 (MQ=255)
                                    ccAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc                                                     >  1:4329453/1‑93 (MQ=255)
                                    ccAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc                                                     <  1:2425474/93‑1 (MQ=255)
                                     cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    >  1:3226513/1‑93 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   >  1:2326805/1‑93 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   <  2:4308881/93‑1 (MQ=255)
                                       aTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  >  1:4703488/1‑93 (MQ=255)
                                          ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac                                               <  2:5026821/93‑1 (MQ=255)
                                              tGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGCAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg                                           >  2:5291993/1‑93 (MQ=255)
                                               gCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc                                          <  2:2464878/93‑1 (MQ=255)
                                                cATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCa                                         >  1:188214/1‑93 (MQ=255)
                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        >  1:78089/1‑93 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    >  1:2211504/1‑93 (MQ=255)
                                                      gATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                                   <  2:4155016/93‑1 (MQ=255)
                                                       aTAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAata                                  <  1:3723839/93‑1 (MQ=255)
                                                        tAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAatat                                 >  2:91062/1‑93 (MQ=255)
                                                         aaCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtt                                <  2:1971746/93‑1 (MQ=255)
                                                         aaCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtt                                <  2:1318235/93‑1 (MQ=255)
                                                          aCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGAGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTc                               >  1:2523319/1‑93 (MQ=255)
                                                           cTGGCAGGCAACGGTTGTCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              <  2:1542839/93‑1 (MQ=255)
                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              >  2:4234993/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  2:2501289/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  1:5311507/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            <  2:5432945/93‑1 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            <  2:1127400/93‑1 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            <  2:2025784/93‑1 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  2:444616/1‑93 (MQ=255)
                                                              gcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGa                           <  1:5838505/93‑1 (MQ=255)
                                                                aggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACg                         >  2:6053920/1‑93 (MQ=255)
                                                                aggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACg                         >  1:5971195/1‑93 (MQ=255)
                                                                 ggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGa                        >  2:1134337/1‑93 (MQ=255)
                                                                  gcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGaa                       >  2:4864567/1‑93 (MQ=255)
                                                                  gcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGaa                       <  2:4011142/93‑1 (MQ=255)
                                                                   caACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAag                      >  2:3721456/1‑93 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     >  1:1805393/1‑93 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     >  1:5361401/1‑93 (MQ=255)
                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    >  1:4998630/1‑93 (MQ=255)
                                                                      cGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt                   <  2:4685644/93‑1 (MQ=255)
                                                                       ggTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGtt                  <  2:1467709/93‑1 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 >  1:1177072/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 >  1:2835031/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 <  1:2605629/93‑1 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  2:1378701/93‑1 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  2:32800/93‑1 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  2:2342941/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  1:5945753/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  2:495917/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  1:2805601/93‑1 (MQ=255)
                                                                              cTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCg           <  1:3522088/93‑1 (MQ=255)
                                                                              cTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCg           <  2:5133257/93‑1 (MQ=255)
                                                                               tGTTTACAGTTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          <  1:1748429/93‑1 (MQ=255)
                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          >  2:3304569/1‑93 (MQ=255)
                                                                                 tttACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa        >  1:860925/1‑93 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACCCCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  1:2797055/93‑1 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  2:5452693/93‑1 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  1:3044036/93‑1 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  2:381342/93‑1 (MQ=255)
                                                                                     cAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAt    <  1:428428/93‑1 (MQ=255)
                                                                                      aGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAtt   >  1:1516316/1‑93 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  >  1:1268023/1‑93 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  <  1:1922581/93‑1 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  <  1:6262966/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ 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
CTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTAAC  >  NC_000913/3539681‑3539866
                                                                                             |                                                                                            
CTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTA                                                                                              >  SRR4302533.3286029/1‑94 (MQ=60)
    CGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAT                                                                                          >  SRR4302533.1795638/1‑94 (MQ=60)
    CGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAT                                                                                          <  SRR4302533.5985597/94‑1 (MQ=60)
    CGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAT                                                                                          <  SRR4302533.4455800/94‑1 (MQ=60)
     GTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATG                                                                                         <  SRR4302533.5660630/94‑1 (MQ=60)
     GTCCGCCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATG                                                                                         <  SRR4302533.3818555/94‑1 (MQ=60)
      TCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGG                                                                                        <  SRR4302533.6052715/94‑1 (MQ=60)
      TCCGTCAGCCGGGTGACCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGG                                                                                        <  SRR4302533.1667169/94‑1 (MQ=60)
       CCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATCACTGGCAGGCAACGCTTGGTCTGTTTACAGGTACCATGGC                                                                                       >  SRR4302533.1983703/1‑94 (MQ=60)
        CGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCG                                                                                      <  SRR4302533.6332091/94‑1 (MQ=60)
        CGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCG                                                                                      >  SRR4302533.3525223/1‑94 (MQ=60)
         GTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGA                                                                                     >  SRR4302533.1594965/1‑94 (MQ=60)
         GTCAGCCGGGAGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGA                                                                                     >  SRR4302533.1846632/1‑94 (MQ=60)
           CAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAA                                                                                   <  SRR4302533.3418023/94‑1 (MQ=60)
           CAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAA                                                                                   <  SRR4302533.789329/94‑1 (MQ=60)
            CGCCGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAG                                                                                  <  SRR4302533.2924032/94‑1 (MQ=60)
            AGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAG                                                                                  <  SRR4302533.3386135/94‑1 (MQ=60)
              CCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAA                                                                                <  SRR4302533.2172280/94‑1 (MQ=60)
              CCGGGGGATCACCCCGGTCTTCAAGACAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAA                                                                                >  SRR4302533.1439641/1‑94 (MQ=60)
               cggatcatcccccctGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAG                                                                               <  SRR4302533.82678/79‑1 (MQ=60)
               CGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAG                                                                               <  SRR4302533.3947784/94‑1 (MQ=60)
                GGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGT                                                                              >  SRR4302533.3224222/1‑94 (MQ=60)
                GGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGT                                                                              >  SRR4302533.3878510/1‑94 (MQ=60)
                 GGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTG                                                                             <  SRR4302533.3849127/94‑1 (MQ=60)
                  GTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGG                                                                            >  SRR4302533.2689103/1‑94 (MQ=60)
                  GTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGG                                                                            >  SRR4302533.3698917/1‑94 (MQ=60)
                   TGATCCCCCCTGTCTGCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGTCGAAAGAAGTGGT                                                                           <  SRR4302533.3230125/94‑1 (MQ=60)
                   TGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGT                                                                           >  SRR4302533.1983813/1‑94 (MQ=60)
                   TGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGT                                                                           <  SRR4302533.4907322/94‑1 (MQ=60)
                   TGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGT                                                                           >  SRR4302533.2338902/1‑94 (MQ=60)
                    GATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTA                                                                          <  SRR4302533.1639762/94‑1 (MQ=60)
                     ATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAG                                                                         >  SRR4302533.3297695/1‑94 (MQ=60)
                     ATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAG                                                                         >  SRR4302533.4852780/1‑94 (MQ=60)
                     ATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAG                                                                         <  SRR4302533.5168396/94‑1 (MQ=60)
                       CACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTG                                                                       <  SRR4302533.5351442/94‑1 (MQ=60)
                       CACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTG                                                                       >  SRR4302533.641998/1‑94 (MQ=60)
                       CACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTG                                                                       <  SRR4302533.3894774/94‑1 (MQ=60)
                        CCCCCGGTCTGCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGG                                                                      <  SRR4302533.4935660/94‑1 (MQ=60)
                        ACCCCGGTCTTCAAGCCAATTGGCGTGCATGATGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGG                                                                      <  SRR4302533.5719376/94‑1 (MQ=60)
                        ACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGG                                                                      <  SRR4302533.1196946/94‑1 (MQ=60)
                         CCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGG                                                                     >  SRR4302533.1588567/1‑94 (MQ=60)
                          CCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGT                                                                    <  SRR4302533.4016666/94‑1 (MQ=60)
                          CCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGT                                                                    >  SRR4302533.3740673/1‑94 (MQ=60)
                          CCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGACAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGT                                                                    >  SRR4302533.3397849/1‑94 (MQ=60)
                           CCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTA                                                                   <  SRR4302533.6170802/94‑1 (MQ=60)
                           CCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTA                                                                   >  SRR4302533.4652003/1‑94 (MQ=60)
                            CGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTAC                                                                  <  SRR4302533.4669535/94‑1 (MQ=60)
                            CGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTAC                                                                  <  SRR4302533.4863992/94‑1 (MQ=60)
                             GGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACG                                                                 >  SRR4302533.1132580/1‑94 (MQ=60)
                              GTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGC                                                                >  SRR4302533.1164456/1‑94 (MQ=60)
                               TCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCT                                                               >  SRR4302533.4319184/1‑94 (MQ=60)
                               TCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCT                                                               >  SRR4302533.2903777/1‑94 (MQ=60)
                                TTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTC                                                              <  SRR4302533.4866528/94‑1 (MQ=60)
                                 TTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCA                                                             >  SRR4302533.3744170/1‑94 (MQ=60)
                                   CAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAAC                                                           <  SRR4302533.1351083/94‑1 (MQ=60)
                                    AAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACA                                                          <  SRR4302533.1054927/94‑1 (MQ=60)
                                    AAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACA                                                          >  SRR4302533.2429595/1‑94 (MQ=60)
                                     AGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAC                                                         <  SRR4302533.3673407/94‑1 (MQ=60)
                                     AGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAC                                                         >  SRR4302533.982356/1‑94 (MQ=60)
                                     AGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAC                                                         >  SRR4302533.3476496/1‑94 (MQ=60)
                                      GCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACC                                                        >  SRR4302533.4355537/1‑94 (MQ=60)
                                       CCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCC                                                       >  SRR4302533.3246406/1‑94 (MQ=60)
                                       CCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCC                                                       <  SRR4302533.2440447/94‑1 (MQ=60)
                                        CAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCT                                                      >  SRR4302533.2341210/1‑94 (MQ=60)
                                         AATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTC                                                     >  SRR4302533.4731999/1‑94 (MQ=60)
                                          ATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCT                                                    <  SRR4302533.1512347/94‑1 (MQ=60)
                                            TGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTAC                                                  <  SRR4302533.5758816/94‑1 (MQ=60)
                                             GGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACA                                                 <  SRR4302533.1317372/94‑1 (MQ=60)
                                             GGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACA                                                 <  SRR4302533.3953262/94‑1 (MQ=60)
                                                GTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCG                                              <  SRR4302533.1390918/94‑1 (MQ=60)
                                                GTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCG                                              >  SRR4302533.5814476/1‑94 (MQ=60)
                                                  GCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCA                                            >  SRR4302533.189384/1‑94 (MQ=60)
                                                   CATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG                                           <  SRR4302533.1061583/94‑1 (MQ=60)
                                                   CATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG                                           >  SRR4302533.78591/1‑94 (MQ=60)
                                                       AAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAA                                       >  SRR4302533.2225141/1‑94 (MQ=60)
                                                       AAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAA                                       >  SRR4302533.1650828/1‑94 (MQ=60)
                                                        AGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAAT                                      >  SRR4302533.3907270/1‑94 (MQ=60)
                                                        AGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAAT                                      <  SRR4302533.2859664/94‑1 (MQ=60)
                                                         GATAACTGGTAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATA                                     >  SRR4302533.489553/1‑94 (MQ=60)
                                                         GATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATA                                     >  SRR4302533.2808672/1‑94 (MQ=60)
                                                          ATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAT                                    <  SRR4302533.3746686/94‑1 (MQ=60)
                                                           TAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATT                                   <  SRR4302533.918678/94‑1 (MQ=60)
                                                            AACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGAGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTC                                  >  SRR4302533.2538864/1‑94 (MQ=60)
                                                             ACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCA                                 <  SRR4302533.1819625/94‑1 (MQ=60)
                                                               TGGCAGGCAACGGTTGGCCTGTTTACAGGTCCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACTCCCTCTACACCGCAGAAAATATTCAGG                               <  SRR4302533.1241235/94‑1 (MQ=60)
                                                               TGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGG                               <  SRR4302533.3722991/94‑1 (MQ=60)
                                                               TGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGG                               <  SRR4302533.2568042/94‑1 (MQ=60)
                                                               TGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGG                               >  SRR4302533.5343338/1‑94 (MQ=60)
                                                                GGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGA                              >  SRR4302533.284747/1‑94 (MQ=60)
                                                                GGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGA                              <  SRR4302533.1924596/94‑1 (MQ=60)
                                                                 GCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGAC                             <  SRR4302533.5873805/94‑1 (MQ=60)
                                                                  CAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACG                            >  SRR4302533.6007341/1‑94 (MQ=60)
                                                                  CAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACG                            <  SRR4302533.3615722/94‑1 (MQ=60)
                                                                      CAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGGCGAAGA                        <  SRR4302533.1676618/94‑1 (MQ=60)
                                                                      CAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGA                        >  SRR4302533.5393601/1‑94 (MQ=60)
                                                                      CAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGA                        >  SRR4302533.1816382/1‑94 (MQ=60)
                                                                       AACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAG                       >  SRR4302533.5028934/1‑94 (MQ=60)
                                                                       AACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAG                       <  SRR4302533.1466633/94‑1 (MQ=60)
                                                                        ACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGT                      <  SRR4302533.2816724/94‑1 (MQ=60)
                                                                          GGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTC                    >  SRR4302533.2852479/1‑94 (MQ=60)
                                                                          GGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTC                    <  SRR4302533.3116335/94‑1 (MQ=60)
                                                                          GGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTC                    >  SRR4302533.1184136/1‑94 (MQ=60)
                                                                          GGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGGGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTC                    >  SRR4302533.3027546/1‑94 (MQ=60)
                                                                           GTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCA                   <  SRR4302533.2621634/94‑1 (MQ=60)
                                                                           GTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCA                   <  SRR4302533.2131689/94‑1 (MQ=60)
                                                                             TGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACCCCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAAT                 <  SRR4302533.3085654/94‑1 (MQ=60)
                                                                             TGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAAT                 <  SRR4302533.608047/94‑1 (MQ=60)
                                                                             TGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAAT                 <  SRR4302533.2458529/94‑1 (MQ=60)
                                                                             TGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAAT                 <  SRR4302533.5955179/94‑1 (MQ=60)
                                                                             TGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAAT                 >  SRR4302533.3786073/1‑94 (MQ=60)
                                                                               GCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCC               <  SRR4302533.2822869/94‑1 (MQ=60)
                                                                               GCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCC               <  SRR4302533.5981710/94‑1 (MQ=60)
                                                                                CCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCG              >  SRR4302533.1000474/1‑94 (MQ=60)
                                                                                 CTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGG             <  SRR4302533.3543798/94‑1 (MQ=60)
                                                                                  TGTTTACAGTTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGC            <  SRR4302533.1759074/94‑1 (MQ=60)
                                                                                   GTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCA           >  SRR4302533.865969/1‑94 (MQ=60)
                                                                                   GTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCA           >  SRR4302533.4625772/1‑94 (MQ=60)
                                                                                      TACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAA        >  SRR4302533.2424409/1‑94 (MQ=60)
                                                                                       ACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACCCCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAT       <  SRR4302533.2814258/94‑1 (MQ=60)
                                                                                       ACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAT       <  SRR4302533.3062670/94‑1 (MQ=60)
                                                                                        CAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATT      <  SRR4302533.430983/94‑1 (MQ=60)
                                                                                        CAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATT      <  SRR4302533.5380069/94‑1 (MQ=60)
                                                                                        CAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATT      >  SRR4302533.1525446/1‑94 (MQ=60)
                                                                                        CAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATT      >  SRR4302533.6199387/1‑94 (MQ=60)
                                                                                         AGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT     >  SRR4302533.1275614/1‑94 (MQ=60)
                                                                                          GGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTA    >  SRR4302533.4743266/1‑94 (MQ=60)
                                                                                          GGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTA    <  SRR4302533.6300402/94‑1 (MQ=60)
                                                                                          GGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTA    <  SRR4302533.1934315/94‑1 (MQ=60)
                                                                                            TACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTAAC  >  SRR4302533.278569/1‑94 (MQ=60)
                                                                                            TACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTAAC  >  SRR4302533.3671382/1‑94 (MQ=60)
                                                                                            TACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTAAC  <  SRR4302533.2413379/94‑1 (MQ=60)
                                                                                            TACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTAAC  <  SRR4302533.2310870/94‑1 (MQ=60)
                                                                                            TACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTAAC  >  SRR4302533.972549/1‑94 (MQ=60)
                                                                                             |                                                                                            
CTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTTAAC  >  NC_000913/3539681‑3539866

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