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 3,263,409 C→T 100% intergenic (+448/+166) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,4090CT92.6% 236.6 / 7.8 103intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base T (29/62);  total (35/68)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 106 total aligned reads displayed.

TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGA  >  CP009273/3263264‑3263554
                                                                                                                                                 |                                                                                                                                                 
ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAtttt                                                                                                                                                 <  1:1466804/148‑1 (MQ=255)
 tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                                <  2:2300457/148‑1 (MQ=37)
 tttGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATACGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                                <  2:1545477/148‑1 (MQ=255)
  ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                               <  2:765609/148‑1 (MQ=37)
  ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                               <  2:765646/148‑1 (MQ=37)
  ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                               <  2:766537/148‑1 (MQ=37)
       tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGccc                                                                                                                                          >  2:590048/1‑148 (MQ=37)
       tataGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCCAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTATTTTGGGTATAGTGTCGTGGACAGTGAATAAATTTTCTGccc                                                                                                                                          >  1:253135/1‑148 (MQ=11)
         taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCt                                                                                                                                        <  1:1135861/148‑1 (MQ=37)
          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                       <  2:972112/148‑1 (MQ=37)
          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                       <  2:848846/148‑1 (MQ=37)
          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                       <  2:804238/148‑1 (MQ=37)
          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                       >  2:3115914/1‑148 (MQ=255)
          aGTGTCGTGGACAGTCAATCATCTTTATGCCCCTCCAAAAGTAAAAACCCCCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCGTTTGGGTATAGTGTGGTGGACAGTCCTTCATTTTTGGGCCCCTc                                                                                                                                       >  1:3478105/1‑148 (MQ=11)
                 tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGc                                                                                                                                <  1:619683/148‑1 (MQ=37)
                 tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGc                                                                                                                                <  1:619686/148‑1 (MQ=37)
                     cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                            <  2:114001/148‑1 (MQ=37)
                     cAGTCATTCATCTTTCTGCCACTCCACAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                            <  1:2557736/148‑1 (MQ=255)
                      aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGGCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaaa                                                                                                                           <  1:1914541/148‑1 (MQ=255)
                       gTCATTCATCTTTCGGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAc                                                                                                                          <  1:2418047/148‑1 (MQ=255)
                        tcattcatCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                         <  1:3058756/148‑1 (MQ=255)
                        tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                         <  1:3306762/148‑1 (MQ=37)
                        tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                         <  1:3058739/148‑1 (MQ=37)
                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACccg                                                                                                                       <  2:3176292/148‑1 (MQ=37)
                           ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGACGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACccgc                                                                                                                      <  1:377246/148‑1 (MQ=255)
                            tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                                                                                                     <  1:1930861/148‑1 (MQ=37)
                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                                                                                   <  2:148421/148‑1 (MQ=37)
                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                                                                                  >  1:2565004/1‑148 (MQ=37)
                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCCTTCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                           >  2:3517331/1‑148 (MQ=255)
                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                        <  1:674287/148‑1 (MQ=37)
                                            ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTGCTTCATTTTTATGCCGGTCCAAAAGCTAATACACGCCGAAGCGGGTTGTTACg                                                                                                     >  2:399777/1‑148 (MQ=11)
                                            ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTCTGCCCCTCCAAAAGCAAATACCCGCCGCCGCGGGTTTTTACg                                                                                                     >  2:2269427/1‑148 (MQ=11)
                                            ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGTGGACAGTCATTCATGTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGGTGTTaag                                                                                                     >  2:2269194/1‑146 (MQ=11)
                                                aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                 <  2:481939/148‑1 (MQ=37)
                                                   tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCCGGCATTCATTTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGCAGCGGGTGTTTACGTAAACCa                                                                                              >  1:2460773/1‑148 (MQ=11)
                                                     aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                                                                            <  1:524383/148‑1 (MQ=37)
                                                        aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCATTCATCTTTATGCCCCACCAAAAGCAAAAAGCCGCCGACGCGGGTTTTTACGTGAACCAGGTGa                                                                                         >  2:3280503/1‑148 (MQ=11)
                                                        aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                                                                                                         >  2:3280540/1‑100 (MQ=37)
                                                        aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                                                                                                         >  2:1202538/1‑100 (MQ=37)
                                                              cgAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                   <  2:1482743/148‑1 (MQ=255)
                                                              cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                   <  2:930158/148‑1 (MQ=37)
                                                              cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                   <  2:930102/148‑1 (MQ=37)
                                                              cgAAGCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                   <  1:1834247/148‑1 (MQ=255)
                                                               gAAGCGGGTTTTGACGTAAATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGTGTCGTGGACAGTCATTCAGTTTTCTGCCCCTCCAAACGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                                                                  <  2:725919/148‑1 (MQ=11)
                                                                aaGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                                                 <  2:1756969/148‑1 (MQ=37)
                                                                 agggggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCg                                                                                <  2:1266827/145‑1 (MQ=12)
                                                                     ggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATaa                                                                            >  2:2371375/1‑148 (MQ=17)
                                                                      gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                                           <  1:590048/148‑1 (MQ=37)
                                                                      gTTTTTACGTAAAACAGGTGAAACTGACCGATAACCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCCCAGCGGGGTTTTACGTAAACCAGGTGAAAATTACCCATAAg                                                                           >  2:2819822/1‑148 (MQ=11)
                                                                        ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGTAACTGACCGATGAGcc                                                                         >  1:525818/1‑148 (MQ=255)
                                                                        ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                                         >  2:579293/1‑148 (MQ=37)
                                                                        ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAACCAAAAACCCGCCCAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                                         >  1:2942869/1‑148 (MQ=255)
                                                                         tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                                        <  1:2014665/148‑1 (MQ=37)
                                                                         tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                                        <  1:1171860/148‑1 (MQ=37)
                                                                            aCGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                                     <  2:83436/148‑1 (MQ=255)
                                                                             cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGcttt                                                                    <  1:3292042/148‑1 (MQ=255)
                                                                              gTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                   <  1:2951480/148‑1 (MQ=255)
                                                                               tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                  >  2:285663/1‑148 (MQ=25)
                                                                                aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAAATTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                 >  1:203831/1‑148 (MQ=255)
                                                                                 aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                                <  2:2176627/148‑1 (MQ=255)
                                                                                      ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTCTGCCCCTCCAAAAGCCAAAACCCGCCGCAGCGGGTTTTTACGTTAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTgtgt                                                           >  2:1822567/1‑148 (MQ=11)
                                                                                      ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                                           >  2:1089717/1‑148 (MQ=12)
                                                                                      ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTGTTgggg                                                           >  2:1822676/1‑147 (MQ=255)
                                                                                        taaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                         <  1:2054165/146‑1 (MQ=17)
                                                                                        tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATCGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                         <  1:3119443/148‑1 (MQ=25)
                                                                                          aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAg                                                       >  2:76659/1‑148 (MQ=255)
                                                                                            aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcg                                                     <  1:3120963/148‑1 (MQ=25)
                                                                                            aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcg                                                     <  1:3119828/148‑1 (MQ=25)
                                                                                            aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcg                                                     <  1:2850327/148‑1 (MQ=25)
                                                                                            aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATGTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcc                                                     <  1:3119824/148‑2 (MQ=25)
                                                                                            aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAAATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcg                                                     <  1:1347763/148‑1 (MQ=25)
                                                                                             cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                    <  1:1126505/148‑1 (MQ=37)
                                                                                             cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                    <  2:1566133/148‑1 (MQ=37)
                                                                                              tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                   <  1:837112/148‑1 (MQ=37)
                                                                                              tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                   <  1:1722056/148‑1 (MQ=37)
                                                                                                  cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                               <  1:1161977/148‑1 (MQ=255)
                                                                                                  cGATAAGCCGCATTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                               <  1:1161969/148‑1 (MQ=255)
                                                                                                  cGACAAGCCGCTTTCTTTTGGCTATAGTGACGTGGACAGTCATTCCTTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAACCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                               <  1:1161967/148‑1 (MQ=11)
                                                                                                    aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTTCCCCTCCAAAATCAAAAACCCCCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACAGAACGATAATCCGCTTTCTTTTTGGTATAGCGTCGTGGAc                                             >  2:3052499/1‑148 (MQ=11)
                                                                                                    aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                             <  2:2333357/148‑1 (MQ=37)
                                                                                                      aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCGTGTACAg                                           >  2:1377298/1‑148 (MQ=255)
                                                                                                        gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAAATTTTCTGCCCCACCAAAAGCAAAAACCCCCCGAATCGGGTTTTTTAGTAAACCAGGTGAAACAGACCGATAATCCGCTTTATTTTGGGTATATCGTCGTGGGCAAtc                                         >  2:906443/1‑148 (MQ=11)
                                                                                                         ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtca                                        >  2:3143729/1‑148 (MQ=37)
                                                                                                          cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                                                                                                          <  2:3240270/65‑1 (MQ=35)
                                                                                                             tttctttTGGGTATAGTGTCGTGGACAGTCTTTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca                                    <  1:1107413/148‑1 (MQ=255)
                                                                                                               tctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATc                                  <  1:573030/148‑1 (MQ=37)
                                                                                                               tctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATc                                  <  2:2565004/148‑1 (MQ=37)
                                                                                                                  tttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                               <  1:3115914/148‑1 (MQ=255)
                                                                                                                  tttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                               >  2:802377/1‑148 (MQ=37)
                                                                                                                  cttGGGTATAGTGTCGTGGACAGTCAATCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGATTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                               <  2:1911221/147‑1 (MQ=255)
                                                                                                                         ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGcccc                        <  1:3280503/148‑1 (MQ=37)
                                                                                                                           cgtgtCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                      <  2:924924/147‑1 (MQ=255)
                                                                                                                           aGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAACTCAATCATCTTTCTGCCCCTc                      >  2:3425469/1‑148 (MQ=255)
                                                                                                                           aGTGTCGTGGACAGTCAATCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACACGCAATCATCATTCTAccccaa                      >  2:3425651/1‑146 (MQ=11)
                                                                                                                                  ttgACAGTCATTCTTTTGTCTGACCCTCCAAAAGCAAATATCCTCCGATGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGc               <  2:2087179/146‑1 (MQ=11)
                                                                                                                                     aCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGGTTTCTGTTGGGTATAGCGACGTGGACAGTCCATCATCTTTGATCCCCTCCAAAAGCaaa            >  2:3258283/1‑148 (MQ=11)
                                                                                                                                       aGTCATTCATTTTTCGTCCCCTCCACTAGCATAAACCCGCCGTATCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaaa          <  2:1311716/148‑1 (MQ=11)
                                                                                                                                         tcattgatTTTACGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTATTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc        <  2:1240062/148‑1 (MQ=11)
                                                                                                                                              catTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATCTTTATTCCCCACCAAAAACAAAAACccgccg   >  2:2758400/1‑148 (MQ=11)
                                                                                                                                               atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTGTGGGGTATAGCGTCGTGGACAGTCATTCATGATAGAGCCCCTCCAACACCAAAAACccgccgc  >  2:1771158/1‑147 (MQ=11)
                                                                                                                                                 |                                                                                                                                                 
TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGA  >  CP009273/3263264‑3263554

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
CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGC  >  CP009273/3263263‑3263557
                                                                                                                                                  |                                                                                                                                                     
CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTT                                                                                                                                                    <  SRR19973067.68664/150‑1 (MQ=32)
CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTT                                                                                                                                                    <  SRR19973067.1879864/150‑1 (MQ=32)
 TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTC                                                                                                                                                   <  SRR19973067.1476588/150‑1 (MQ=31)
         ATAGTGTCGTGGACAGTCAATCATCTTTATGCCCCTCCAAAAGTAAAAACCCCCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCGTTTGGGTATAGTGTGGTGGACAGTCCTTCATTTTTGGGCCCCTC                                                                                                                                           >  SRR19973067.3500185/1‑150 (MQ=24)
          TAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCC                                                                                                                                          <  SRR19973067.1143512/150‑1 (MQ=32)
                 GTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTTTTATTCCCCACCAAAAGCA                                                                                                                                   >  SRR19973067.1109868/1‑150 (MQ=28)
                  TGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAA                                                                                                                                  <  SRR19973067.623762/150‑1 (MQ=32)
                  TGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAA                                                                                                                                  <  SRR19973067.623759/150‑1 (MQ=32)
                    GACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGACCGGCCTTAATTTTTATGCCCCGCCAAAAGCAAAA                                                                                                                                >  SRR19973067.372467/1‑150 (MQ=26)
                      CAGTCATTCATCTTTCTGCCACTCCACAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAC                                                                                                                              <  SRR19973067.2574456/150‑1 (MQ=29)
                       AGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGGCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACC                                                                                                                             <  SRR19973067.1927611/150‑1 (MQ=31)
                        GTCATTCATCTTTCGGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCC                                                                                                                            <  SRR19973067.2433885/150‑1 (MQ=31)
                        GTCATTCATCGTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCC                                                                                                                            <  SRR19973067.2434120/150‑1 (MQ=31)
                         TCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCG                                                                                                                           <  SRR19973067.3078231/150‑1 (MQ=31)
                         TCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCG                                                                                                                           <  SRR19973067.3327625/150‑1 (MQ=32)
                         TCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCG                                                                                                                           <  SRR19973067.3078214/150‑1 (MQ=32)
                            TTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGACGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCG                                                                                                                        <  SRR19973067.379868/150‑1 (MQ=31)
                             TCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCAGTCATTTTTCTGCCCCTACAAAAGCAAAAACCCGCCGA                                                                                                                       >  SRR19973067.2305764/1‑150 (MQ=29)
                             TCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGA                                                                                                                       <  SRR19973067.1944033/150‑1 (MQ=32)
                              CATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAA                                                                                                                      >  SRR19973067.2581757/1‑150 (MQ=31)
                                        CCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTT                                                                                                            <  SRR19973067.3474847/150‑1 (MQ=32)
                                          CCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAC                                                                                                          <  SRR19973067.678919/150‑1 (MQ=32)
                                          CCCCCAAAAGCAAAAACACGCCGAAGCGGTTATTTACGTATAACATGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTAC                                                                                                          <  SRR19973067.337609/150‑1 (MQ=0)
                                               AAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTGTTGGGTATAGAGGCGAGGACAGTAATTCATTATGCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTCAA                                                                                                    >  SRR19973067.1057843/1‑150 (MQ=16)
                                                AAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCCGTCATTCATATTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAA                                                                                                   >  SRR19973067.98615/1‑150 (MQ=29)
                                                AAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGTGGACAGTGCGTAATTTTTATGCCCCACCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTTAAA                                                                                                    >  SRR19973067.2109568/1‑150 (MQ=25)
                                                  AGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCCGGCATTCATTTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGCAGCGGGTGTTTACGTAAACCA                                                                                                 >  SRR19973067.2476954/1‑150 (MQ=29)
                                                   gacaatacgtgactggtgctcagtcttgggctcaaccgatatACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAACCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAG                                                                                                <  SRR19973067.1217253/108‑1 (MQ=0)
                                                      AAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTG                                                                                             <  SRR19973067.527916/150‑1 (MQ=32)
                                                            cgccgtataggtTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGGTTTCTTTTGGGTTTTGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCGGAAGCGGGTTTTTACGTAAACCAGGTGAAACTG                                                                                       <  SRR19973067.339534/138‑1 (MQ=28)
                                                               CGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACC                                                                                    <  SRR19973067.1846859/150‑1 (MQ=32)
                                                               CGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACC                                                                                    <  SRR19973067.2511925/150‑1 (MQ=32)
                                                               CGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACC                                                                                    <  SRR19973067.2511910/150‑1 (MQ=32)
                                                               CGAAGCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACC                                                                                    <  SRR19973067.1846731/150‑1 (MQ=31)
                                                                GAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCG                                                                                   <  SRR19973067.1048649/150‑1 (MQ=32)
                                                                GAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCG                                                                                   <  SRR19973067.1048129/150‑1 (MQ=32)
                                                                       GTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTAATTTTTGGGCCCCTACAAAAACAAAAAACCCCCGCAGCGGGTTTTTACGTAAACCAGGGGAAACagtcccagaaacc                                                                            >  SRR19973067.2408642/1‑137 (MQ=46)
                                                                       GTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGTAACTGACCGATGAGCC                                                                            >  SRR19973067.529359/1‑150 (MQ=29)
                                                                       GTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCC                                                                            <  SRR19973067.593807/150‑1 (MQ=32)
                                                                       GTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAACCAAAAACCCGCCCAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCC                                                                            >  SRR19973067.2961516/1‑150 (MQ=58)
                                                                          TTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCT                                                                         <  SRR19973067.2028212/150‑1 (MQ=36)
                                                                          TTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCT                                                                         <  SRR19973067.1179629/150‑1 (MQ=36)
                                                                              CGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCT                                                                     <  SRR19973067.3312744/150‑1 (MQ=44)
                                                                               GTAAATCAGGTGAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTT                                                                    <  SRR19973067.3467521/150‑1 (MQ=28)
                                                                               GTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTT                                                                    <  SRR19973067.2970152/150‑1 (MQ=44)
                                                                               GTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAAATTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTT                                                                    >  SRR19973067.205375/1‑150 (MQ=21)
                                                                                TAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCCCCCAAGCGGGTTTTTACGTAAAACAGGTGTAACAGACCGATAATCCCCGTTCTTT                                                                   >  SRR19973067.473188/1‑150 (MQ=17)
                                                                                TAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTT                                                                   <  SRR19973067.3287240/150‑1 (MQ=42)
                                                                                  AAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGCAGCGGGGGTTTTCGTAAACCCGGGGAAACTGACCGATAAGCCGCgtgggtttg                                                                 >  SRR19973067.1602512/1‑141 (MQ=25)
                                                                                      AGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTA                                                             >  SRR19973067.1488572/1‑150 (MQ=32)
                                                                                      AGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCACCAAAATCAAAAACCCGCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGACCGATAATCCGCGTTCTTTTGGGgt                                                             >  SRR19973067.3359151/1‑148 (MQ=28)
                                                                                       GGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCCTTGATTTTTATGCCCCGCCAAAAACAAAAACCCCCCGCAGCGGGTTTTTAAGTAAAACAGGTGAAAaatacccctaagccccgttgttttggggtt                                                            >  SRR19973067.2731534/1‑120 (MQ=1)
                                                                                         TGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATCGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAG                                                          <  SRR19973067.3139275/150‑1 (MQ=25)
                                                                                         TGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAG                                                          <  SRR19973067.2067869/150‑1 (MQ=25)
                                                                                             ACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTC                                                      <  SRR19973067.3140801/150‑1 (MQ=29)
                                                                                             ACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTC                                                      <  SRR19973067.3139662/150‑1 (MQ=29)
                                                                                             ACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTC                                                      <  SRR19973067.2868611/150‑1 (MQ=29)
                                                                                             ACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATGTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCCTC                                                      <  SRR19973067.3139658/150‑1 (MQ=22)
                                                                                             ACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAAATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTC                                                      <  SRR19973067.1356981/150‑1 (MQ=29)
                                                                                              CTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCG                                                     <  SRR19973067.1134112/150‑1 (MQ=32)
                                                                                               TGGCCGATAAGCCGCTTTGTTTTGGGTATAGTGTCGTGGGCAGTCAATAATTTTTCTTCCCCCCCAAAAGCAAAAAaccccccaagccgggtttttagtaaaacacgtggaacagaacgagaaacagcattatttttggtttttgggagt                                                     >  SRR19973067.625692/1‑76 (MQ=6)
                                                                                               TGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGT                                                    <  SRR19973067.842802/150‑1 (MQ=32)
                                                                                               TGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGT                                                    <  SRR19973067.1733854/150‑1 (MQ=32)
                                                                                                   CGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAC                                                <  SRR19973067.1169709/150‑1 (MQ=30)
                                                                                                   CGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAC                                                <  SRR19973067.3157065/150‑1 (MQ=32)
                                                                                                   CGATAAGCCGCATTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAC                                                <  SRR19973067.1169701/150‑1 (MQ=30)
                                                                                                   CGACAAGCCGCTTTCTTTTGGCTATAGTGACGTGGACAGTCATTCCTTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAACCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAC                                                <  SRR19973067.1169699/150‑1 (MQ=27)
                                                                                                         GCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCATAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAT                                          <  SRR19973067.420053/150‑1 (MQ=30)
                                                                                                         GCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAACGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAT                                          <  SRR19973067.1002490/150‑1 (MQ=31)
                                                                                                         GCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAT                                          <  SRR19973067.2222008/150‑1 (MQ=31)
                                                                                                          CCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATT                                         <  SRR19973067.2901228/150‑1 (MQ=32)
                                                                                                              TTTCTTTTGGGTATAGTGTCGTGGACAGTCTTTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATC                                     <  SRR19973067.1114935/150‑1 (MQ=31)
                                                                                                                TCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTT                                   <  SRR19973067.576734/150‑1 (MQ=32)
                                                                                                                   TTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                <  SRR19973067.3135724/150‑1 (MQ=31)
                                                                                                                   TTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                <  SRR19973067.1821414/150‑1 (MQ=32)
                                                                                                                   TTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                <  SRR19973067.1821424/150‑1 (MQ=32)
                                                                                                                          ATAGTGTCGTGGTCAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTC                         <  SRR19973067.1210399/150‑1 (MQ=31)
                                                                                                                          ATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTTGGTATAGCGTCGTGGTCAGTCCTTCATATTTGTGCCCCGC                         >  SRR19973067.2591835/1‑150 (MQ=60)
                                                                                                                          ATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTC                         <  SRR19973067.3301104/150‑1 (MQ=32)
                                                                                                                          ATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTC                         <  SRR19973067.3301067/150‑1 (MQ=32)
                                                                                                                               GTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCCATCATCTTTATGCCCCTCCAAAA                    >  SRR19973067.2039280/1‑150 (MQ=29)
                                                                                                                                            ATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCC       <  SRR19973067.39232/150‑1 (MQ=32)
                                                                                                                                                 TTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGGCAGGCATTCATagttgtggcccgccaaaatcaaaaaacccccgcagc  >  SRR19973067.1519394/1‑114 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGC  >  CP009273/3263263‑3263557

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: