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 F368 I0 R1 1472 214.4 7116750 95.9% 6824963 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 724,813 A→T 100% intergenic (+17/‑226) ybfA → / → rhsC DUF2517 family protein/Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273724,8130AT99.5% 623.8 / ‑6.4 205intergenic (+17/‑226)ybfA/rhsCDUF2517 family protein/Rhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base T (82/119);  total (82/123)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 205 total aligned reads displayed.

TGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGA  >  CP009273/724669‑724956
                                                                                                                                                |                                                                                                                                               
 gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCAtt                                                                                                                                             >  1:3482514/1‑148 (MQ=255)
 gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCAtt                                                                                                                                             <  1:1332150/148‑1 (MQ=255)
 gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCAtt                                                                                                                                             >  2:584443/1‑148 (MQ=255)
  ccTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTg                                                                                                                                            >  1:1748644/1‑148 (MQ=255)
     ttCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGaaa                                                                                                                                         >  2:3030820/1‑148 (MQ=255)
     ttCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGaaa                                                                                                                                         >  1:2557769/1‑148 (MQ=255)
       cATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAAc                                                                                                                                       <  2:2478202/148‑1 (MQ=255)
          gCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGc                                                                                                                                    >  2:2915189/1‑148 (MQ=255)
          gCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGATACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGc                                                                                                                                    >  2:632505/1‑148 (MQ=255)
           cTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGcc                                                                                                                                   >  2:3259189/1‑148 (MQ=255)
             tttCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAg                                                                                                                                 <  1:512822/148‑1 (MQ=255)
             tttCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAg                                                                                                                                 >  2:3392294/1‑148 (MQ=255)
                cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCgg                                                                                                                              <  2:1511856/148‑1 (MQ=255)
                  ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaa                                                                                                                            <  1:780864/148‑1 (MQ=255)
                  ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaa                                                                                                                            <  1:937730/148‑1 (MQ=255)
                   gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                           >  2:1509046/1‑148 (MQ=255)
                   gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                           <  2:3068639/148‑1 (MQ=255)
                      aGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGc                                                                                                                        <  1:309747/148‑1 (MQ=255)
                      aGACCGCGCCCGCTTTTACAGCTACCTGAATCGCGTCTGGTCGAAAACCATCGATAAACCGGTGTGGATGTATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGc                                                                                                                        <  1:3363125/148‑1 (MQ=255)
                       gACCGCGCCCGCATTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAACAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCt                                                                                                                       <  2:2744363/148‑1 (MQ=255)
                        aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTg                                                                                                                      <  2:2723595/148‑1 (MQ=255)
                           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCg                                                                                                                   >  1:2018537/1‑148 (MQ=255)
                              cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                <  1:3460671/148‑1 (MQ=255)
                               ccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGACTTCATTGAAAACCGCCAGCGGAAACGCGGGCGGttt                                                                                                               >  2:1964761/1‑148 (MQ=255)
                                cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtttt                                                                                                              <  2:597205/148‑1 (MQ=255)
                                cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttt                                                                                                              <  2:1711697/148‑1 (MQ=255)
                                 gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTGTTGAATTCATTGAAATCCGCCAGCGGAAACGCTGGCGGtgttt                                                                                                             >  2:47859/1‑148 (MQ=255)
                                   ttttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTg                                                                                                           >  1:650553/1‑148 (MQ=255)
                                     ttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTa                                                                                                         <  2:3480226/148‑1 (MQ=255)
                                          gCTACCTGCATCGCGTGTGGTCGAAAACCAGCGATATACCTTTGTGTATGTTTCAGGCCGAAAAGGCTACTTGTGATTTTTATTTTTTTCTTAGTTTTGAATTCATTGAAAACCTCCAGCGGAAACGCTGGCGGTTTTTTGTATCAgg                                                                                                    <  2:1412371/148‑1 (MQ=255)
                                           cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCATGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTAATTGAAAACCGCCATCGGAAACGCTGGCGGTTTTTTTTATAAGGa                                                                                                   >  1:169171/1‑148 (MQ=255)
                                             aCCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGTATGGATCAGGCCTAAAAGGAAACTGTTGATTTTTATTGATTGCTGAGTTTTGATTTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAg                                                                                                 <  1:3448424/148‑1 (MQ=255)
                                             aCCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAg                                                                                                 <  2:609317/148‑1 (MQ=255)
                                              ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGa                                                                                                <  1:1398080/148‑1 (MQ=255)
                                                 gCATCGCGTCTGGTCTAAAACCAGCTATAAACCGGTGTGTATGGATCCGGCCGAAAAGGCAAGTGGTGATTTTTATTGATTGCTGAGTTTTTAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGt                                                                                             <  1:717901/148‑1 (MQ=255)
                                                 gCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGt                                                                                             <  2:3226555/148‑1 (MQ=255)
                                                    tCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                          >  1:1109649/1‑148 (MQ=255)
                                                      gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                        <  2:2962698/148‑1 (MQ=255)
                                                       cgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagag                                                                                       <  1:2331654/148‑1 (MQ=255)
                                                          cTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGaa                                                                                    <  2:3150067/148‑1 (MQ=255)
                                                             gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAatat                                                                                 <  1:2743481/148‑1 (MQ=255)
                                                              tCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                >  2:1832942/1‑148 (MQ=255)
                                                               cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                               >  1:207748/1‑148 (MQ=255)
                                                                 aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                             <  1:515501/148‑1 (MQ=255)
                                                                 aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                             >  2:2361371/1‑148 (MQ=255)
                                                                 aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                             <  2:210758/148‑1 (MQ=255)
                                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCg                                                                       <  1:1949444/148‑1 (MQ=255)
                                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCg                                                                       <  2:3233362/148‑1 (MQ=255)
                                                                        gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                      <  2:1227314/148‑1 (MQ=255)
                                                                        gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                      <  1:427329/148‑1 (MQ=255)
                                                                           aTAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtat                                                                   <  2:1839639/148‑1 (MQ=255)
                                                                            tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtata                                                                  >  2:875489/1‑148 (MQ=255)
                                                                            tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtata                                                                  <  2:2661083/148‑1 (MQ=255)
                                                                             aaaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGATGTTAAGAGGAATATTTAGACATAGTtataa                                                                 >  1:2443143/1‑147 (MQ=255)
                                                                              aaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                >  2:1129186/1‑148 (MQ=255)
                                                                                ccGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                              <  2:2627080/148‑1 (MQ=255)
                                                                                 cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAAt                                                             <  1:2915189/148‑1 (MQ=255)
                                                                                      tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                        <  1:1233789/148‑1 (MQ=255)
                                                                                       ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                       <  1:103451/148‑1 (MQ=255)
                                                                                         atggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                     >  1:1450960/1‑148 (MQ=255)
                                                                                           ggatCAGGCCGAAACGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatat                                                   <  1:983217/148‑1 (MQ=255)
                                                                                           ggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatat                                                   <  2:2297228/148‑1 (MQ=255)
                                                                                            gatCAGGCCGAAAAGTCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTTAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCATGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                  <  2:2048355/148‑1 (MQ=255)
                                                                                            gatCAGGCCGAAAAGGCAACTGTTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                  <  1:538084/148‑1 (MQ=255)
                                                                                            gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                  <  2:1189154/148‑1 (MQ=255)
                                                                                            gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                  <  2:3302448/148‑1 (MQ=255)
                                                                                             atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatatat                                                 <  2:1177731/148‑1 (MQ=255)
                                                                                             atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatatat                                                 >  2:358518/1‑148 (MQ=255)
                                                                                              tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                <  1:96620/148‑1 (MQ=255)
                                                                                               cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTc                                               <  2:739042/148‑1 (MQ=255)
                                                                                                 ggCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTAAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTAAt                                             >  2:727894/1‑148 (MQ=255)
                                                                                                    cGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaa                                          <  2:1344577/148‑1 (MQ=255)
                                                                                                        aaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatat                                      >  2:167524/1‑148 (MQ=255)
                                                                                                          ggCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatat                                    <  1:2700755/148‑1 (MQ=255)
                                                                                                          ggCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatat                                    >  2:861230/1‑148 (MQ=255)
                                                                                                           gCAACTTGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatata                                   <  2:1887526/148‑1 (MQ=255)
                                                                                                           gCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatata                                   <  2:2002116/148‑1 (MQ=255)
                                                                                                             aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAAtatatata                                 >  2:2736833/1‑148 (MQ=255)
                                                                                                               cTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaaa                               >  2:2957242/1‑148 (MQ=255)
                                                                                                               cTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATAGTTATACCAATCGAAACATATATTCATGAAATATATATaaa                               >  2:1246560/1‑148 (MQ=255)
                                                                                                                tGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGAATTTCAGAGGAAAATTTATACAAAGTTATACCAATCCAAACAAAAAATAAATAAAAATATATaaaa                              >  1:1048930/1‑147 (MQ=255)
                                                                                                                   tGATTTTTATGGCTTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAtt                           <  1:2607537/148‑1 (MQ=255)
                                                                                                                    gATTTTTTTGGATTTCTGTGTTTTGAATTCTCTGAAAACCGCCAGCGGAAAGTCTGGGGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACAAATCGAAACATTTATTCATGAAATATATATAAATAttt                          <  1:2592433/148‑1 (MQ=255)
                                                                                                                    gATTTTTATTGTTTGCTGAGTTTTGAATACATAGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAttt                          <  1:1849695/148‑1 (MQ=255)
                                                                                                                    gATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAttt                          >  2:528005/1‑148 (MQ=255)
                                                                                                                      tttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                        <  1:3259189/148‑1 (MQ=255)
                                                                                                                       ttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTcc                       <  1:873462/148‑1 (MQ=255)
                                                                                                                           attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATTAAATATATATAAAAATTTTaataat                   >  2:2451143/1‑148 (MQ=255)
                                                                                                                           attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAt                   <  1:3255303/148‑1 (MQ=255)
                                                                                                                            ttgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATAAGGAATATTTTCAGTGGAAAAATTATACAAAGTTATAAAAATAGAAAAAAAAAATAAATAAATAAATATAAAAAttttaaaaaat                  >  1:410650/1‑141 (MQ=255)
                                                                                                                               attgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGtt               <  2:1126692/148‑1 (MQ=255)
                                                                                                                                ttgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGtta              >  2:2832629/1‑147 (MQ=255)
                                                                                                                                  gCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCtt            <  1:2962605/148‑1 (MQ=255)
                                                                                                                                   cTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTa           >  2:1445005/1‑148 (MQ=255)
                                                                                                                                       gTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGTAGATGTTAAGAGGAATATTTTGAAATAGTTATAACCAACGAAACAAATATTCATGAAATATATATAAATATTTTACTAATTGTTATTATATg       >  2:1209044/1‑148 (MQ=255)
                                                                                                                                        ttttGAATTCATTCAAAACCGCCAGCTGAATCGGTGGCGGTTTTTTGTATCATGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTAATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGa      <  1:903660/148‑1 (MQ=255)
                                                                                                                                           tGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg   <  2:1109649/148‑1 (MQ=255)
                                                                                                                                           tGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg   <  2:160457/148‑1 (MQ=255)
                                                                                                                                           tGAATTAATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATAGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg   <  1:584493/148‑1 (MQ=255)
                                                                                                                                            gAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTCATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGa  <  1:3391407/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGA  >  CP009273/724669‑724956

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATA  >  CP009273/724668‑724958
                                                                                                                                                 |                                                                                                                                                 
CTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATT                                                                                                                                               >  SRR19973067.3504612/1‑150 (MQ=60)
 TGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTG                                                                                                                                              <  SRR19973067.371910/150‑1 (MQ=60)
 TGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTG                                                                                                                                              >  SRR19973067.781291/1‑150 (MQ=60)
 TGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTG                                                                                                                                              >  SRR19973067.1760567/1‑150 (MQ=60)
  GCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCTCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGA                                                                                                                                             <  SRR19973067.1341246/150‑1 (MQ=60)
    CTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAA                                                                                                                                           >  SRR19973067.2574489/1‑150 (MQ=60)
    CTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGAATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAA                                                                                                                                           >  SRR19973067.3275569/1‑150 (MQ=60)
         ATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGC                                                                                                                                      >  SRR19973067.2326206/1‑150 (MQ=60)
              TTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCG                                                                                                                                 <  SRR19973067.516206/150‑1 (MQ=60)
              TTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCG                                                                                                                                 <  SRR19973067.516312/150‑1 (MQ=60)
              TTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCG                                                                                                                                 <  SRR19973067.516769/150‑1 (MQ=60)
               TTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGG                                                                                                                                <  SRR19973067.172673/150‑1 (MQ=60)
                   GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAC                                                                                                                            <  SRR19973067.786268/150‑1 (MQ=60)
                   GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAC                                                                                                                            <  SRR19973067.944103/150‑1 (MQ=60)
                   GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAC                                                                                                                            <  SRR19973067.2968281/150‑1 (MQ=60)
                   GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTAGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAC                                                                                                                            <  SRR19973067.84436/150‑1 (MQ=60)
                       AGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTG                                                                                                                        <  SRR19973067.311769/150‑1 (MQ=60)
                       AGACCGCGCCCGCTTTTACAGCTACCTGAATCGCGTCTGGTCGAAAACCATCGATAAACCGGTGTGGATGTATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTG                                                                                                                        <  SRR19973067.3384473/150‑1 (MQ=60)
                          CCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCACCTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCG                                                                                                                     >  SRR19973067.2032101/1‑150 (MQ=60)
                          CCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCG                                                                                                                     >  SRR19973067.2032100/1‑150 (MQ=60)
                               CCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT                                                                                                                <  SRR19973067.3482671/150‑1 (MQ=60)
                                CCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTT                                                                                                               >  SRR19973067.1106016/1‑150 (MQ=60)
                                 CGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTT                                                                                                              <  SRR19973067.539349/150‑1 (MQ=60)
                                 CGATTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCTATAAACCGGTGTGGATGGATCAGGCCTAAAAGGAAACTGGTGATTTTTATTGATTGCTGATTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTT                                                                                                              <  SRR19973067.2422437/150‑1 (MQ=60)
                                  GCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCGGGCGGTTTTTTT                                                                                                             >  SRR19973067.654966/1‑150 (MQ=60)
                                    TTTTACCGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTA                                                                                                           <  SRR19973067.2850556/150‑1 (MQ=60)
                                      TTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTtttta                                                                                                         >  SRR19973067.935508/1‑145 (MQ=60)
                                         CAGCTACCTGCATCGCGTCTGTTCTAAAACCAGCTATAAACGGGTGTGGATTTATCAGGCCTAAAAGGAAAGTGGTGATTTTTATTGATTGCTGATTTTTTAATTAATTTAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGG                                                                                                      <  SRR19973067.699891/150‑1 (MQ=60)
                                         CAGCTACCTGCATCGCGTCTGGTCCAAAACCAGCTTTAACCCGTTGTGGATGGATCAGGCCGAAAAGTCAACTGGTTATTTTTATTGATTGCTGCGTTTTGATTTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGG                                                                                                      <  SRR19973067.699780/150‑1 (MQ=60)
                                          AGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCATGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTAATTGAAAACCGCCATCGGAAACGCTGGCGGTTTTTTTTATAAGGA                                                                                                     >  SRR19973067.170549/1‑150 (MQ=60)
                                              ACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGTATGGATCAGGCCTAAAAGGAAACTGTTGATTTTTATTGATTGCTGAGTTTTGATTTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAT                                                                                                 <  SRR19973067.3470363/150‑1 (MQ=60)
                                               CCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATG                                                                                                <  SRR19973067.1407570/150‑1 (MQ=60)
                                                  GCATCGCGTCTGGTCTAAAACCAGCTATAAACCGGTGTGTATGGATCCGGCCGAAAAGGCAAGTGGTGATTTTTATTGATTGCTGAGTTTTTAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTC                                                                                             <  SRR19973067.722878/150‑1 (MQ=60)
                                                   CATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGATGTTCA                                                                                            >  SRR19973067.1117177/1‑150 (MQ=60)
                                                    ATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAG                                                                                           <  SRR19973067.1144583/150‑1 (MQ=60)
                                                        CGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTTTGAGGAAGATtgtcataggc                                                                                       >  SRR19973067.194061/1‑140 (MQ=60)
                                                        CGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGA                                                                                       <  SRR19973067.2346948/150‑1 (MQ=60)
                                                             GGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATT                                                                                  <  SRR19973067.2620768/150‑1 (MQ=60)
                                                              GTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTT                                                                                 <  SRR19973067.2761205/150‑1 (MQ=60)
                                                              GTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTT                                                                                 >  SRR19973067.209309/1‑150 (MQ=60)
                                                                  AAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAC                                                                             <  SRR19973067.518999/150‑1 (MQ=60)
                                                                        AGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTTTAAGGAAGATttttagagtaatatatataaaaacgt                                                                       >  SRR19973067.2729117/1‑124 (MQ=60)
                                                                        AGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTT                                                                       <  SRR19973067.1962722/150‑1 (MQ=60)
                                                                         GCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTTTTTTTATTGATTGCTTAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTA                                                                      <  SRR19973067.430226/150‑1 (MQ=60)
                                                                         GCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTA                                                                      <  SRR19973067.430296/150‑1 (MQ=60)
                                                                          CTATAACCCGGTGTGGATGGATCAGGCCGAAAAGCCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTAT                                                                     <  SRR19973067.3018083/150‑1 (MQ=60)
                                                                            ATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGATGTTAAGAGGAATATTTAGACATAGTTATAA                                                                   >  SRR19973067.2459188/1‑150 (MQ=60)
                                                                              AAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATGAGGAAGATGTTCAGAGGAATATTTATACATAGTTATAACA                                                                 >  SRR19973067.2760083/1‑150 (MQ=60)
                                                                                  CGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCG                                                             <  SRR19973067.2933743/150‑1 (MQ=60)
                                                                                       TGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACA                                                        <  SRR19973067.1242063/150‑1 (MQ=60)
                                                                                        GGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACAT                                                       <  SRR19973067.104485/150‑1 (MQ=60)
                                                                                        GGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACAT                                                       >  SRR19973067.1460680/1‑150 (MQ=60)
                                                                                         GATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATA                                                      <  SRR19973067.3413874/150‑1 (MQ=60)
                                                                                            GTATCAGGCAGAAAAGCAAACTGTTGATTTTTATTGATTGCTGTTTTTTGAATTCATTCAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAT                                                   <  SRR19973067.315678/150‑1 (MQ=60)
                                                                                            GGATCAGGCCGAAACGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAT                                                   <  SRR19973067.989972/150‑1 (MQ=60)
                                                                                             GATCAGGCCGAAAAGGCAACTGTTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATT                                                  <  SRR19973067.541676/150‑1 (MQ=60)
                                                                                             GATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATT                                                  <  SRR19973067.1788738/150‑1 (MQ=60)
                                                                                               TCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCA                                                <  SRR19973067.97621/150‑1 (MQ=60)
                                                                                                     CGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTttataatgaatatgttaaaaggaatatttataaaacgttataacaaaagaaaaaaatattaaataaaa                                          >  SRR19973067.438807/1‑82 (MQ=60)
                                                                                                     CGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTtttatgaagatttttacaggaaaatttttaaaaaggtataaaaaaacaaaaaaaaaataaataaat                                          >  SRR19973067.936914/1‑84 (MQ=60)
                                                                                                      GAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATA                                         <  SRR19973067.1258598/150‑1 (MQ=60)
                                                                                                        AAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTTTCAGGAAGATGTTCAGATGAATATTTAGAAATCGTTATAACAATCCAAACAAAAATTCATTAAAAATA                                       >  SRR19973067.3476296/1‑150 (MQ=60)
                                                                                                          AGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCCGCGGAAACGCTGGCGGTTTTTTtttttatgaatatgtttatagttaaaattttacaaagttataacaaaagaaaaaaaaaataaataaaaatata                                     >  SRR19973067.1023/1‑77 (MQ=60)
                                                                                                           GGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATAT                                    <  SRR19973067.2718129/150‑1 (MQ=60)
                                                                                                           GGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATAT                                    <  SRR19973067.2594026/150‑1 (MQ=60)
                                                                                                            GCATCTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTTAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATA                                   <  SRR19973067.981741/150‑1 (MQ=60)
                                                                                                            GCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATA                                   <  SRR19973067.2772082/150‑1 (MQ=60)
                                                                                                               ACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGAATTTCAGAGGAAAATTTATACAAAGTTATACCAATCCAAACAaaaaataaataaaaatatataaaa                                >  SRR19973067.1056130/1‑126 (MQ=60)
                                                                                                               ACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATaatggagatgttgagaggcataattagtaaacggtataacaaacgaaaaatatattaaagaaaaatatataaaa                                >  SRR19973067.1056096/1‑76 (MQ=60)
                                                                                                                  GGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATT                             <  SRR19973067.1408044/150‑1 (MQ=60)
                                                                                                                    TGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTT                           <  SRR19973067.2628098/150‑1 (MQ=60)
                                                                                                                    TGATTTTTATGGCTTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTT                           <  SRR19973067.2624476/150‑1 (MQ=60)
                                                                                                                     GATTTTTTTGGATTTCTGTGTTTTGAATTCTCTGAAAACCGCCAGCGGAAAGTCTGGGGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACAAATCGAAACATTTATTCATGAAATATATATAAATATTTTC                          <  SRR19973067.2609309/150‑1 (MQ=60)
                                                                                                                     GATTTTTATTGTTTGCTGAGTTTTGAATACATAGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTC                          <  SRR19973067.1862306/150‑1 (MQ=60)
                                                                                                                     GATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTC                          <  SRR19973067.1619784/150‑1 (MQ=60)
                                                                                                                       TTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTGGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCT                        <  SRR19973067.3355160/150‑1 (MQ=60)
                                                                                                                       TTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCT                        <  SRR19973067.3279552/150‑1 (MQ=60)
                                                                                                                       TTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCT                        <  SRR19973067.3408116/150‑1 (MQ=60)
                                                                                                                        TTTTTTTGATTGCTGTGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTA                       <  SRR19973067.879294/150‑1 (MQ=60)
                                                                                                                        TTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTA                       <  SRR19973067.881332/150‑1 (MQ=60)
                                                                                                                           TATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTttataaggaatattttcagtggaaaaattatacaaagttataaaaatagaaaaaaaaaataaataaataaatataaaaattttaaaaaat                    >  SRR19973067.413509/1‑60 (MQ=60)
                                                                                                                            ATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTG                   <  SRR19973067.3275683/150‑1 (MQ=60)
                                                                                                                            ATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTG                   <  SRR19973067.3275652/150‑1 (MQ=60)
                                                                                                                              TGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATAAGGAATATGTTAAGAGGAAAATTTAGAAAAAGTTATAACAATCGAAAAAAATATTAATTAAAAATATATAAAAATTTTacaaattttt                 >  SRR19973067.2534949/1‑140 (MQ=56)
                                                                                                                                 TTGCTGTGTTTTTAATTCATTGAAAACCGCCAGCGTAACGGCGTGGGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTT              <  SRR19973067.636740/150‑1 (MQ=60)
                                                                                                                                 TTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTT              <  SRR19973067.3487446/150‑1 (MQ=60)
                                                                                                                                   GCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTAT            <  SRR19973067.2981317/150‑1 (MQ=60)
                                                                                                                                    CTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATAAGGAAGATGTTAAGAGGAAaaattagaaaaacttataacaaaacaaacaaaaattaattaaatatatataaaaaatttaataaattttattata           >  SRR19973067.1348839/1‑75 (MQ=60)
                                                                                                                                       AGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGA        <  SRR19973067.195469/150‑1 (MQ=60)
                                                                                                                                         TTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGCTCTTATCTGACA      <  SRR19973067.2032029/150‑1 (MQ=60)
                                                                                                                                         TTTTGAATTCATTCAAAACCGCCAGCTGAATCGGTGGCGGTTTTTTGTATCATGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTAATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACA      <  SRR19973067.909627/150‑1 (MQ=60)
                                                                                                                                         TTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGGAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACA      <  SRR19973067.909541/150‑1 (MQ=60)
                                                                                                                                            TGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGAT   <  SRR19973067.965792/150‑1 (MQ=60)
                                                                                                                                            TGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGAT   <  SRR19973067.588184/150‑1 (MQ=60)
                                                                                                                                            TGAATTAATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATAGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGAT   <  SRR19973067.588234/150‑1 (MQ=60)
                                                                                                                                             GAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTCATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATA  <  SRR19973067.3412983/150‑1 (MQ=60)
                                                                                                                                                 |                                                                                                                                                 
CTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATA  >  CP009273/724668‑724958

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