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 F159 I0 R1 879 305.9 9588012 96.6% 9262019 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 581,203 A→T 100% intergenic (‑114/‑310) ompT ← / → pauD DLP12 prophage; outer membrane protease VII (outer membrane protein 3b)/tRNA‑OTHER

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273581,2030AT95.8% 497.0 / 9.3 141intergenic (‑114/‑310)ompT/pauDDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)/tRNA‑OTHER
Reads supporting (aligned to +/- strand):  ref base A (1/5);  new base T (72/63);  total (73/68)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 141 total aligned reads displayed.

GTCAGGACTATTCCCAGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCAT  >  CP009273/581058‑581345
                                                                                                                                                 |                                                                                                                                              
gTCAGGACTATTCCCAGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTcac                                                                                                                                            >  2:314775/1‑150 (MQ=255)
gTCAGGACTATTCCCAGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCAACTGTTTATATTTTGTTTAGTATTGAATGATTTTcac                                                                                                                                            >  2:314774/1‑150 (MQ=255)
 tCAGGACTATTCCCAGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTcaca                                                                                                                                           >  2:2176585/1‑150 (MQ=255)
             ccAGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTg                                                                                                                               <  2:2260106/150‑1 (MQ=255)
               aGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGtt                                                                                                                             <  1:3871685/150‑1 (MQ=255)
                  aGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGttttt                                                                                                                          <  1:4251187/150‑1 (MQ=255)
                  aGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGttttt                                                                                                                          <  1:4251209/150‑1 (MQ=255)
                  aGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGttttt                                                                                                                          >  2:3215092/1‑150 (MQ=255)
                  aGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGttttt                                                                                                                          >  2:3213270/1‑150 (MQ=255)
                   gTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAAGTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTg                                                                                                                         <  2:2113544/150‑1 (MQ=255)
                     ttCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGt                                                                                                                       >  1:1267555/1‑150 (MQ=255)
                      tCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGtt                                                                                                                      <  1:1549065/150‑1 (MQ=255)
                      tCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGtt                                                                                                                      <  1:1915631/150‑1 (MQ=255)
                         cccGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTaa                                                                                                                   >  1:3542079/1‑150 (MQ=255)
                         cccGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTaa                                                                                                                   <  2:4507627/150‑1 (MQ=255)
                           cGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAAtt                                                                                                                 >  1:1823407/1‑150 (MQ=255)
                            gCATAAAAGTTCTCCATTCAATCTTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTa                                                                                                                <  2:2965150/150‑1 (MQ=255)
                            gCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTa                                                                                                                <  2:2965737/150‑1 (MQ=255)
                                aaaaGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTg                                                                                                            <  2:766133/150‑1 (MQ=255)
                                    gTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTAt                                                                                                        >  2:4727183/1‑150 (MQ=255)
                                      tctcCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTg                                                                                                      >  2:202537/1‑150 (MQ=255)
                                      tctcCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTg                                                                                                      <  1:74538/150‑1 (MQ=255)
                                      tctcCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTg                                                                                                      <  1:76728/150‑1 (MQ=255)
                                         ccATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGttt                                                                                                   <  2:4191935/150‑1 (MQ=255)
                                          cATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTc                                                                                                  <  2:502471/150‑1 (MQ=255)
                                          cATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTc                                                                                                  <  1:1674313/150‑1 (MQ=255)
                                          cATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTc                                                                                                  <  2:4161236/150‑1 (MQ=255)
                                           ataaaaTCGTATTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCa                                                                                                 <  2:4002327/146‑1 (MQ=255)
                                           aTTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCa                                                                                                 <  2:4001283/150‑1 (MQ=255)
                                              cAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCatat                                                                                              >  1:1226919/1‑150 (MQ=255)
                                              cAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCatat                                                                                              <  1:330474/150‑1 (MQ=255)
                                              cAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCatat                                                                                              <  1:330519/150‑1 (MQ=255)
                                               aaTCGTTTTAATGATTGAGTATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATAtt                                                                                             >  1:977473/1‑150 (MQ=255)
                                               aaTCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATAtt                                                                                             <  2:3766641/150‑1 (MQ=255)
                                                 tCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTc                                                                                           <  2:3692334/150‑1 (MQ=255)
                                                 tCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTc                                                                                           <  1:4320399/150‑1 (MQ=255)
                                                 acGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTc                                                                                           <  1:4320503/149‑1 (MQ=255)
                                                       tAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGtttt                                                                                     >  2:2721328/1‑150 (MQ=255)
                                                         aTGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTAc                                                                                   >  1:4430677/1‑150 (MQ=255)
                                                          tGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACt                                                                                  <  1:873362/150‑1 (MQ=255)
                                                             ttGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACtgtg                                                                               <  1:4587097/150‑1 (MQ=255)
                                                               gAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACtgtgtg                                                                             <  1:2753266/150‑1 (MQ=255)
                                                               gAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACtgtgtg                                                                             <  1:2525902/150‑1 (MQ=255)
                                                                aaTATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTgg                                                                            >  1:3788491/1‑150 (MQ=255)
                                                                aaTATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTgg                                                                            >  2:1938349/1‑150 (MQ=255)
                                                                         tttttATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTAt                                                                   >  1:1110989/1‑150 (MQ=255)
                                                                          ttttATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATg                                                                  >  1:1696849/1‑150 (MQ=255)
                                                                           tttATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTGTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGc                                                                 <  2:4738748/150‑1 (MQ=255)
                                                                            ttATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTTGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCt                                                                >  1:2917950/1‑150 (MQ=255)
                                                                            ttATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCt                                                                >  1:2917589/1‑150 (MQ=255)
                                                                             tataTCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCtt                                                               <  2:1740196/150‑1 (MQ=255)
                                                                             tataTCTAACTTAATGAGTCAATTACATAGTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCtt                                                               <  2:1740208/150‑1 (MQ=255)
                                                                              atatCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGtttt                                                              >  2:4569202/1‑150 (MQ=255)
                                                                              atatCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCttt                                                              >  2:2061900/1‑150 (MQ=255)
                                                                              atatCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCttt                                                              >  1:3780428/1‑150 (MQ=255)
                                                                              atatCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATAGTTCGGTTTTACTGTGTGGTTTTTTTATGCttt                                                              >  2:2062865/1‑150 (MQ=255)
                                                                                atCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTtg                                                            >  1:1418212/1‑150 (MQ=255)
                                                                                atCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTtg                                                            <  1:445482/150‑1 (MQ=255)
                                                                                atCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTtg                                                            >  2:3104859/1‑150 (MQ=255)
                                                                                atCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTGTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTtg                                                            >  2:552034/1‑150 (MQ=255)
                                                                                   tAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTgtgg                                                         >  1:694494/1‑150 (MQ=255)
                                                                                   tAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTgtgg                                                         >  1:4612521/1‑150 (MQ=255)
                                                                                   tAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTgtgg                                                         >  2:2826167/1‑150 (MQ=255)
                                                                                      cTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTTTGCTTTTGTGGTGc                                                      >  2:147737/1‑150 (MQ=255)
                                                                                       ttAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCt                                                     >  1:3003289/1‑150 (MQ=255)
                                                                                       ttAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCt                                                     >  1:3140411/1‑150 (MQ=255)
                                                                                         aaTGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCttt                                                   >  2:2321894/1‑150 (MQ=255)
                                                                                           tGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGtt                                                                          <  2:2873987/125‑1 (MQ=255)
                                                                                           tGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGtt                                                                          >  1:2873987/1‑125 (MQ=255)
                                                                                           tGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTa                                                 <  1:3156770/150‑1 (MQ=255)
                                                                                             aGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATc                                               <  1:2099438/150‑1 (MQ=255)
                                                                                               tCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTa                                             <  2:1016270/150‑1 (MQ=255)
                                                                                                  aTTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTAttt                                          <  1:1874073/150‑1 (MQ=255)
                                                                                                    tACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTaa                                        >  2:2159771/1‑150 (MQ=255)
                                                                                                        tatTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGc                                    >  2:2320424/1‑150 (MQ=255)
                                                                                                        tatTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGc                                    <  2:2643138/150‑1 (MQ=255)
                                                                                                           tGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCAt                                 <  2:2574515/150‑1 (MQ=255)
                                                                                                              tCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGcc                              >  1:3646154/1‑150 (MQ=255)
                                                                                                              tCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGcc                              <  1:1583435/150‑1 (MQ=255)
                                                                                                              tCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGcc                              <  1:3990075/150‑1 (MQ=255)
                                                                                                                    gTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTaga                        >  2:1558454/1‑150 (MQ=255)
                                                                                                                     tttATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTagag                       >  1:4488859/1‑150 (MQ=255)
                                                                                                                      ttATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAgg                      >  1:2679689/1‑150 (MQ=255)
                                                                                                                       tataTTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCCTTTATCTATTTAAGTGCCATGCCTTTAGAGGc                     >  1:1981764/1‑150 (MQ=255)
                                                                                                                           ttttGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCatat                 >  2:2822093/1‑150 (MQ=255)
                                                                                                                            tttGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCAtata                <  1:2117077/150‑1 (MQ=255)
                                                                                                                             ttGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATaa               >  1:2342499/1‑150 (MQ=255)
                                                                                                                              tGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAg              >  1:670295/1‑150 (MQ=255)
                                                                                                                                tttAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCg            <  2:4430677/150‑1 (MQ=255)
                                                                                                                                 tttgTATTGTATGTTTTTCACTATGCGCTATCTGTTTTTTGTTTAATTATCTGTTAGTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTTGAGGCATATTAGCgt           <  2:696137/147‑2 (MQ=255)
                                                                                                                                 ttAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGa           <  1:2721328/150‑1 (MQ=255)
                                                                                                                                 ttAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGATTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGa           <  1:2716333/150‑1 (MQ=255)
                                                                                                                                  tAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGaa          <  1:59191/150‑1 (MQ=255)
                                                                                                                                    gTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGaaaa        >  2:1340812/1‑150 (MQ=255)
                                                                                                                                    gTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGaaaa        >  2:3754266/1‑150 (MQ=255)
                                                                                                                                      aTTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa      <  1:202537/150‑1 (MQ=255)
                                                                                                                                         gAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCa   >  2:4632005/1‑150 (MQ=255)
                                                                                                                                         gAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCa   >  2:2865501/1‑150 (MQ=255)
                                                                                                                                          aaTGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATATCAt  >  1:4784032/1‑150 (MQ=255)
                                                                                                                                          aaTGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCAt  >  1:4784182/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTCAGGACTATTCCCAGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCAT  >  CP009273/581058‑581345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGT  >  CP009273/581073‑581350
                                                                                                                                  |                                                                                                                                                   
AGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTT                                                                                                                                  <  SRR19973078.3893791/150‑1 (MQ=60)
   AGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTT                                                                                                                               <  SRR19973078.4275391/150‑1 (MQ=60)
   AGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTT                                                                                                                               <  SRR19973078.4275369/150‑1 (MQ=60)
    GTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTG                                                                                                                              <  SRR19973078.3281124/150‑1 (MQ=60)
      TTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGT                                                                                                                            >  SRR19973078.1275429/1‑150 (MQ=60)
       TCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTT                                                                                                                           <  SRR19973078.1558309/150‑1 (MQ=60)
       TCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTT                                                                                                                           <  SRR19973078.447250/150‑1 (MQ=60)
       TCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTT                                                                                                                           <  SRR19973078.1926935/150‑1 (MQ=60)
          CCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAA                                                                                                                        >  SRR19973078.3562092/1‑150 (MQ=60)
            CGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATT                                                                                                                      >  SRR19973078.1834393/1‑150 (MQ=60)
                  AAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGT                                                                                                                >  SRR19973078.2714933/1‑150 (MQ=60)
                  AAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGT                                                                                                                >  SRR19973078.3637079/1‑150 (MQ=60)
                       TCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTG                                                                                                           <  SRR19973078.75317/150‑1 (MQ=60)
                       TCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTG                                                                                                           <  SRR19973078.77522/150‑1 (MQ=60)
                         TCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTT                                                                                                         >  SRR19973078.276471/1‑150 (MQ=60)
                           CATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTC                                                                                                       <  SRR19973078.1684497/150‑1 (MQ=60)
                              TCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATA                                                                                                    >  SRR19973078.1588385/1‑150 (MQ=60)
                               CAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATAT                                                                                                   <  SRR19973078.332678/150‑1 (MQ=60)
                               CAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATAT                                                                                                   <  SRR19973078.332633/150‑1 (MQ=60)
                               CAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATAT                                                                                                   >  SRR19973078.939913/1‑150 (MQ=60)
                               CAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATAT                                                                                                   >  SRR19973078.1234424/1‑150 (MQ=60)
                                AATCGTTTTAATGATTGAGTATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATT                                                                                                  >  SRR19973078.983581/1‑150 (MQ=60)
                                  TCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTC                                                                                                <  SRR19973078.4344787/150‑1 (MQ=60)
                                  ACGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTC                                                                                                <  SRR19973078.4344894/150‑1 (MQ=60)
                                          ATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTAC                                                                                        >  SRR19973078.4455653/1‑150 (MQ=60)
                                           TGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACT                                                                                       <  SRR19973078.878729/150‑1 (MQ=60)
                                              TTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTG                                                                                    <  SRR19973078.4613088/150‑1 (MQ=60)
                                                GAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTG                                                                                  <  SRR19973078.2540390/150‑1 (MQ=60)
                                                GAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTG                                                                                  <  SRR19973078.2768694/150‑1 (MQ=60)
                                                 AATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGG                                                                                 >  SRR19973078.3810236/1‑150 (MQ=60)
                                                     TGTATTTTTTATATCTAACTTAATGAGTCAATTACATTTTGCTCCACTGTTGATATTTTGTTTAGTATTGAATGATTTTCACAATGCGATATCTGTTTTTTGTTTTATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTT                                                                             <  SRR19973078.4563299/150‑1 (MQ=60)
                                                          TTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTAT                                                                        >  SRR19973078.1117728/1‑150 (MQ=60)
                                                           TTTTATATCTAACTTAATGATTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATG                                                                       <  SRR19973078.2101588/150‑1 (MQ=60)
                                                           TTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATG                                                                       >  SRR19973078.1707181/1‑150 (MQ=60)
                                                             TTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTTGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCT                                                                     >  SRR19973078.2934584/1‑150 (MQ=60)
                                                             TTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCT                                                                     >  SRR19973078.2934080/1‑150 (MQ=60)
                                                             TTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCT                                                                     >  SRR19973078.2934221/1‑150 (MQ=60)
                                                             TTATATCTAACTTAATGAGTCAATTACAGATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCT                                                                     <  SRR19973078.3092438/150‑1 (MQ=60)
                                                               ATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACCATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTT                                                                   >  SRR19973078.2333856/1‑150 (MQ=60)
                                                               ATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTT                                                                   >  SRR19973078.3803682/1‑150 (MQ=60)
                                                               ATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTT                                                                   >  SRR19973078.3802125/1‑150 (MQ=60)
                                                                 ATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTG                                                                 >  SRR19973078.1426982/1‑150 (MQ=60)
                                                                 ATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTG                                                                 <  SRR19973078.448304/150‑1 (MQ=60)
                                                                 ATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTG                                                                 >  SRR19973078.1710808/1‑150 (MQ=60)
                                                                    TAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGG                                                              >  SRR19973078.698806/1‑150 (MQ=60)
                                                                    TAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGGGG                                                              >  SRR19973078.4638648/1‑150 (MQ=60)
                                                                        TTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCT                                                          >  SRR19973078.3157996/1‑150 (MQ=60)
                                                                        TTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCT                                                          >  SRR19973078.3020406/1‑150 (MQ=60)
                                                                            TGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTagatcggaagagcacacgtctgaac                                                      >  SRR19973078.2890287/1‑125 (MQ=60)
                                                                            TGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTA                                                      <  SRR19973078.3174400/150‑1 (MQ=60)
                                                                             GAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTAT                                                     <  SRR19973078.607745/150‑1 (MQ=60)
                                                                              AGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATC                                                    <  SRR19973078.2111786/150‑1 (MQ=60)
                                                                                   ATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTT                                               <  SRR19973078.1885251/150‑1 (MQ=60)
                                                                                   ATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTT                                               <  SRR19973078.1357810/150‑1 (MQ=60)
                                                                                    TTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTA                                              <  SRR19973078.4308397/150‑1 (MQ=60)
                                                                                               TCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATTCTTTTGTGGTGCTTTTATCTATTTAATTGACATGCC                                   >  SRR19973078.3666915/1‑150 (MQ=60)
                                                                                               TCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCC                                   <  SRR19973078.1592802/150‑1 (MQ=60)
                                                                                               TCCACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCC                                   <  SRR19973078.4012558/150‑1 (MQ=60)
                                                                                                  ACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTT                                >  SRR19973078.1889800/1‑150 (MQ=60)
                                                                                                  ACTGTTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTT                                <  SRR19973078.728108/150‑1 (MQ=60)
                                                                                                      TTTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAATTGCCATGCCTTTAGAT                            >  SRR19973078.4514265/1‑150 (MQ=60)
                                                                                                       TTATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGG                           >  SRR19973078.2694917/1‑150 (MQ=60)
                                                                                                        TATATTTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCCTTTATCTATTTAAGTGCCATGCCTTTAGAGGC                          >  SRR19973078.1993238/1‑150 (MQ=60)
                                                                                                             TTTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATA                     <  SRR19973078.2129522/150‑1 (MQ=60)
                                                                                                              TTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAA                    >  SRR19973078.2355807/1‑150 (MQ=60)
                                                                                                              TTGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAA                    >  SRR19973078.4767244/1‑150 (MQ=60)
                                                                                                               TGTTTAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAG                   >  SRR19973078.674540/1‑150 (MQ=60)
                                                                                                                  TTATTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGATTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGA                <  SRR19973078.2731660/150‑1 (MQ=60)
                                                                                                                  TTAGTATTTAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGA                <  SRR19973078.2736664/150‑1 (MQ=60)
                                                                                                                   TAGTATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAA               <  SRR19973078.59825/150‑1 (MQ=60)
                                                                                                                      TATTGAATGATTTTCACAATGCGCTATCTGTTTTTGTTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGTTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAAT            <  SRR19973078.4123838/150‑1 (MQ=60)
                                                                                                                       ATTGAATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATA           <  SRR19973078.204195/150‑1 (MQ=60)
                                                                                                                           AATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATATCAT       >  SRR19973078.4810810/1‑150 (MQ=60)
                                                                                                                           AATGATTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCAT       >  SRR19973078.4810960/1‑150 (MQ=60)
                                                                                                                                TTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGT  <  SRR19973078.925580/150‑1 (MQ=60)
                                                                                                                                TTTTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCCTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGT  >  SRR19973078.385675/1‑150 (MQ=60)
                                                                                                                                  |                                                                                                                                                   
AGAAGTTTCGCCCGCATAAAAGTTCTCCATTCAATCGTTTTAATGATTGAATATGTATTTTTTATATCTAACTTAATGAGTCAATTACATATTGCTCCACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGT  >  CP009273/581073‑581350

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