Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,869,637 C→A L87L (CTC→CTA)  cbrA → colicin M resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,869,6370CA96.2% 307.6 / ‑1.4 106L87L (CTC→CTA) cbrAcolicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (46/55);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 106 total aligned reads displayed.

GTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATC  >  NC_000913/3869392‑3869881
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAAcgcg                                                                                                                                                                                                                                                   <  1:847601/251‑1 (MQ=255)
 gtgCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGt                                                                                                                                                                                                                                                  <  1:489491/249‑1 (MQ=255)
    cgattaattgcatCGGCCCGGCAGGATCGGCGTTGGCACGAATGGTAGCAGGCAAAATGCAGGGGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTATCAGCAAACCCTTCTGCGGTCTGCTGGCGCCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACTCTTCCCTTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAACCCGTCGACGTCGCCGCATCGCTAACGCGTA‑Ac                                                                                                                                                                                                                                              <  2:81542/239‑1 (MQ=255)
     gATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACt                                                                                                                                                                                                                                             <  1:421673/251‑1 (MQ=255)
           tGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGc                                                                                                                                                                                                                                       >  1:853050/1‑251 (MQ=255)
            ggCCTCGGCCCGGCAGGAGCGGCGTTGGCACGAAAGTGAGCAGGCAAAATGCAGGTGATCGCGCTGGAGAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATGCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTCCCAGCg                                                                                                                                                                                                                                      <  2:676789/251‑1 (MQ=255)
                   gCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCtat                                                                                                                                                                                                                               <  1:83748/251‑1 (MQ=255)
                       ggCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATc                                                                                                                                                                                                                            >  1:906916/1‑250 (MQ=255)
                             aTCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTa                                                                                                                                                                                                                    <  1:264503/251‑1 (MQ=255)
                               cGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAAt                                                                                                                                                                                                                  <  1:967329/251‑1 (MQ=255)
                                ggCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCCA‑TATTAATc                                                                                                                                                                                                                 >  1:701660/1‑251 (MQ=255)
                                 gCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCGAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCg                                                                                                                                                                                                                <  2:646147/251‑1 (MQ=255)
                                   gTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCCCCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTTACGCCTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCCGCGTCAAAACCTTCGACGTCCCCGCCTCGCTTACGCGAC‑ACTACCAGCGAAGCTATATCCA‑TATTAATCGcc                                                                                                                                                                                                              >  2:24233/1‑251 (MQ=255)
                                      ggCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATg                                                                                                                                                                                                           <  1:50087/251‑1 (MQ=255)
                                      ggCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCGGGATAAAAAGCACCAGTGTGGTACGGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGAAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATg                                                                                                                                                                                                           <  2:588232/251‑1 (MQ=255)
                                          cgcaaGTTAGCAGGCAAAATGCAGGAGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTAAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGAAGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGAACGCTAATCCGCAGATTTTCAGCGTCCAAACCGTAGACGTCGCCGCATCGCTAACGCGTA‑ACTACCATCGACGCTATATCAA‑TATTAATCCCCagtcttt                                                                                                                                                                                                       >  2:785493/4‑245 (MQ=255)
                                           gAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTc                                                                                                                                                                                                      <  1:458783/251‑1 (MQ=255)
                                             aaGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCg                                                                                                                                                                                                     >  1:509293/1‑250 (MQ=255)
                                                    cAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGAAGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACCCGTA‑ACTACCCGCAAAGCCATATCAA‑TATTAATCGCCATGCTTTCGACTTGTgg                                                                                                                                                                                             >  2:612601/1‑251 (MQ=255)
                                                      ggCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATAGCCATGCTTTCGACTTGTggct                                                                                                                                                                                           >  1:546674/1‑249 (MQ=255)
                                                        cAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATg                                                                                                                                                                                          >  1:428545/1‑250 (MQ=255)
                                                          aaaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGaa                                                                                                                                                                                        >  1:300648/1‑250 (MQ=255)
                                                           aaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGGAGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCATCGACGTCGCCGCATCGCTAACGCGTACAC‑ACCAGCGAAGCTATATCAC‑TATTAATCGCCATGCTTTCGACTTGTgggtgcaa                                                                                                                                                                                       >  2:762356/1‑244 (MQ=255)
                                                               cAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCt                                                                                                                                                                                  >  1:214545/1‑251 (MQ=255)
                                                                    gATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGCTTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATcc                                                                                                                                                                             >  1:966448/1‑251 (MQ=255)
                                                                    gATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATcc                                                                                                                                                                             <  1:494463/251‑1 (MQ=255)
                                                                    gATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCGTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATcc                                                                                                                                                                             <  1:212704/251‑1 (MQ=255)
                                                                     aTCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATccc                                                                                                                                                                            >  1:420816/1‑251 (MQ=255)
                                                                      tCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATcccc                                                                                                                                                                           >  1:926602/1‑251 (MQ=255)
                                                                       cgcgCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGACGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCAAGCTTTCGACTTGTGGATGAAATCGCTGATccccc                                                                                                                                                                          >  2:966163/1‑250 (MQ=255)
                                                                         cgcTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGCAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                        >  2:193931/1‑251 (MQ=255)
                                                                         cgcTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                        >  1:197141/1‑251 (MQ=255)
                                                                                 aaaaaGCACCAGTGTGGTACTGACGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGCCGTCGCCGCCTCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TCTTAATCGCCATGCTTTCGTCTTGTGGATGAAATCGCTGATCCCCGCCAGCATTgg                                                                                                                                                                >  2:377772/1‑250 (MQ=255)
                                                                                  aaaaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCGGATCCCCGCCAGCGTTGAg                                                                                                                                                               <  1:221150/251‑1 (MQ=255)
                                                                                         ccAGTGTGGTACTGAAGGTTTCAGCAACCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATc                                                                                                                                                        <  2:150040/251‑1 (MQ=255)
                                                                                         ccAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATc                                                                                                                                                        <  1:795826/251‑1 (MQ=255)
                                                                                           aGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCAc                                                                                                                                                      >  1:88398/1‑251 (MQ=255)
                                                                                           aGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAAAGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCAc                                                                                                                                                      >  1:278769/1‑251 (MQ=255)
                                                                                              gtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGGGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGa                                                                                                                                                    >  1:753568/1‑250 (MQ=255)
                                                                                                   aCTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCt                                                                                                                                              <  1:895792/251‑1 (MQ=255)
                                                                                                      gAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCTATTTTTTATTCGCGAGGGGCTGACGCTTCCTGTCGATGTGATCGCTAACCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATAGCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCtgt                                                                                                                                            <  1:599061/250‑1 (MQ=255)
                                                                                                      gAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGGGGGTTATCCCGACAGCCtgtg                                                                                                                                           >  2:466317/1‑251 (MQ=255)
                                                                                                             tcatcaaccccTGCGGCGGTCTGCTGGCACCGGTCGGGCCGCGATTTTTTTTTCGGGATGGGCTGACGCTTCCGGGCGAGGTGATCGCGAATCCCCAGATTTTCAGCTTCAAACCCGTCGACGTCGCCTCACCGCTCACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCGTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGaaa                                                                                                                                     <  2:509293/242‑1 (MQ=255)
                                                                                                              cAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCTATCTTTTAGGCGCGATGGGCTGACGCTTCCTGGCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCATCGAAGCTATATCAA‑TATTAAGCGCCATGCTTCCGACTTGTGGATGAAATCTCTGAGCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGaaaaa                                                                                                                                   <  1:68232/251‑1 (MQ=255)
                                                                                                                     cccgcggggggTCTGCTGGCACCGGACGGGCAGCGTTCTTTTATTCGGAATGGGCTGACGCTTCCTGTCGAGGTGATCGCTAATCCGCAGATTTGCAGCGACAAAACCGTCGACGTCGCCGCACCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGAGGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCg                                                                                                                            <  2:830912/242‑1 (MQ=255)
                                                                                                                           ggcgtTCTGCTGGCACCGGACTGGCAGCTATTTTTGATGCGCGTTGGGCTGACGCTTCCCGTCGATGAGAACGCTAATCCGCAGATTTTGAGCGTCAAACCCTTCGACGTCGCCGCAACGCTAACGCGTA‑ACTCCCAGCGAATCTA‑CTCAAATATTATCCGCCATGCTTTCGACTTTGGGATGAAATCGCTGATCCCCTCCATCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAATTTTGGCGTGAGga                                                                                                                      <  2:554555/251‑1 (MQ=255)
                                                                                                                           ggcggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCGGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGga                                                                                                                      <  1:445873/251‑1 (MQ=255)
                                                                                                                                 ctgctgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATa                                                                                                                 >  1:90222/1‑250 (MQ=255)
                                                                                                                                   gctgGCACCGGACGCGCAGAGATCTTTTATTCGCGATGGGCTGACGGTTCCTGCCGAGGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGGCGACGTCGCCGCATCGCTACCGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAAGCGCCATGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATa                                                                                                                 <  1:235397/248‑1 (MQ=255)
                                                                                                                                     tgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATgg                                                                                                            >  1:934632/1‑251 (MQ=255)
                                                                                                                                     gggCACCGGACGGGCAGCGATCTTTTATTCGCGATGGGCGGACGCTTCCTGTCGATGGGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAAGCGCCATGCTTTCGACTTGGGGAGGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATgg                                                                                                            <  2:482618/250‑1 (MQ=255)
                                                                                                                                     gggCACCGGACGCGCAGCGATTTTTTTGTCGCGAGGGGCTGACGCTTCCTGTCGATGGGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATgg                                                                                                            <  2:853050/250‑1 (MQ=255)
                                                                                                                                       gCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCTTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCa                                                                                                          <  1:137988/251‑1 (MQ=255)
                                                                                                                                           cGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCCA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAACTCTGGCGTGAGGATGAAAAATGGCATGTc                                                                                                      >  2:524005/1‑251 (MQ=255)
                                                                                                                                             gACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGAGGTTAAATGGCATGTCAt                                                                                                    >  2:809595/1‑251 (MQ=255)
                                                                                                                                                  gcAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTc                                                                                               <  1:698265/251‑1 (MQ=255)
                                                                                                                                                   cAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGGTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCg                                                                                              <  1:822126/251‑1 (MQ=255)
                                                                                                                                                     gCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGACCGCTAATCCGCAGATTTTCAGCGTCAAACCCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATGAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTTAGTATGATAAATGGCATGTCATTTTTCGTg                                                                                            <  1:755063/251‑1 (MQ=255)
                                                                                                                                                         tCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTAAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGc                                                                                           >  1:645464/1‑248 (MQ=255)
                                                                                                                                                           ttttatcggggtggggttGACGCTTCCTTGCGATGGGAGCGTTAATCCGCAGATTTCAAGCGTAAAAACCGTCGACTTCCCCGAATCGCTACCGGGGA‑ACTCCCAGCGAAGCT‑TAACAAATATTAATCGCCAGGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGAc                                                                                       <  2:56916/233‑1 (MQ=255)
                                                                                                                                                                  cgcgATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTggg                                                                                <  1:24233/250‑1 (MQ=255)
                                                                                                                                                                   gcgATGGGTTGACGCTCCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAg                                                                              <  1:466317/251‑1 (MQ=255)
                                                                                                                                                                    cgATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGCAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTgggggc                                                                             >  2:154711/1‑248 (MQ=255)
                                                                                                                                                                    cgATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGc                                                                             >  1:124297/1‑251 (MQ=255)
                                                                                                                                                                    cgATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGc                                                                             >  2:163701/1‑251 (MQ=255)
                                                                                                                                                                     gATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATACCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGc                                                                             >  1:484373/1‑250 (MQ=255)
                                                                                                                                                                      aTGGGCGGACGCTTCCTGTCTGTGTGATCGCTAATCCGCAGATGTTCAGCTTCAAGTCCTTGGACGGCGCCGCATCGCTAACGCGTA‑ACTACCAGAGAAGCTATATCAA‑TATTAATCGCCATGTTTTCGACTTGTGGATGAAATCGCTGATCCACGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAGTTTTCGTGCAGACGGCTGGGAGCat                                                                           <  2:331298/251‑2 (MQ=255)
                                                                                                                                                                       tGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAC‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCACCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGCTAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCaa                                                                           >  2:874913/1‑250 (MQ=255)
                                                                                                                                                                       tGGGCTGACGCTCCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCaa                                                                           <  2:529255/250‑1 (MQ=255)
                                                                                                                                                                        gggCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGCAAAATCTGGCGTGAGGATGATAACTGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAAc                                                                          >  2:318239/1‑250 (MQ=255)
                                                                                                                                                                         ggCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGGCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACat                                                                        <  1:706762/251‑1 (MQ=255)
                                                                                                                                                                           cTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCAATGCTTTCGAATTGTGGATGAAATCGCTGATCCACGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGTTGATAAATGGAATGTCATTTTTCGTGCAGACGACTGGGGGCAACatat                                                                      >  2:422114/1‑251 (MQ=255)
                                                                                                                                                                             gACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATAtt                                                                     <  1:498040/250‑1 (MQ=255)
                                                                                                                                                                                 cTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGGCGTCGACGCATCGATAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGAGTTGAGGTTTATCAAGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCCTTTTTCGTGCAGACGGCTGGGAGCACAATATTAccgc                                                                >  1:95132/1‑251 (MQ=255)
                                                                                                                                                                                 cTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgc                                                                >  1:75828/1‑251 (MQ=255)
                                                                                                                                                                                  ttCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgcc                                                               >  1:4238/1‑251 (MQ=255)
                                                                                                                                                                                    cctggcgatgagtctgtaaacccGCATATTTTCATCGTAAAAGCCTTCGCCGGCGCCGCAGCGCGACCGCGTA‑CCTGCTAGCGAGGCTATATCAA‑TATAAATCGCCATGCGTGCGTTTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgccc                                                              <  2:841129/229‑1 (MQ=255)
                                                                                                                                                                                    ccTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCTTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgccc                                                              <  1:749029/250‑1 (MQ=255)
                                                                                                                                                                                     cTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGCGGTTTATCACGACAGCCTGTGCCGAAAAATCCGGCGTGGGGATGATAAGTGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgcccgc                                                            >  2:739120/1‑251 (MQ=255)
                                                                                                                                                                                       gTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTCCCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCt                                                           <  1:612601/250‑1 (MQ=255)
                                                                                                                                                                                       gTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTa                                                          <  1:318239/251‑1 (MQ=255)
                                                                                                                                                                                           atgttatcgcgcaccctcAGATCTTCAACGTCAAAACCGTCGCCGTCGCCCCATCGCTACCGCTTA‑TCTACCAGCGAACTTATATCAA‑TATTAACCGCCGTGCTGCCGACTTGTGTATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCTCCCGCTATCTg                                                      <  1:250409/234‑1 (MQ=255)
                                                                                                                                                                                              tgATCGCTAATCCGCAGATTTTCAGCGTCAAAAACGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATACA‑TATTAATCGCCATGCTTTCGACATGTGGATGAAATCGCTGTTCCCCGCCAGCGTTGAGGTTTCTCACTACAGCCTGTCCCGCACAATCTGGCGTGAGGATGATACATGGCATGTCATTTTTCGTGCAGACGGCGGGGAGAAAAATATTAACGCCCGCTATCAGGTc                                                   >  2:574309/1‑251 (MQ=255)
                                                                                                                                                                                                aTCGCTAATCCGCCGATTTTCAGCGTCAAAACCGTCGCCGTCGCCGCATCGCTCACGCGTA‑ACTACCAGCGAAGCTATCTCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGCTCCCCGCCAGCGTTGAGGATTATCACGACATCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATAGCATGTCATTTTTCGGGCAGACTGCTGGGGGCACAATATTACCGCCCGCGCTCTGGtagg                                                 >  2:618100/1‑248 (MQ=255)
                                                                                                                                                                                                  cGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTg                                               >  1:766045/1‑251 (MQ=255)
                                                                                                                                                                                                   gctaacccgaagtttttaagactCAAAACCGTGGCCGCCGCCGAATCGCTAACGCTAA‑ACTACCAGCGAAGCTATATCAA‑TATTAACCGCCATGCTTGCGACTGTTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGTAGCAACATATTACCGCCCGCTATCTGGTCGGTGc                                              <  2:936288/229‑1 (MQ=255)
                                                                                                                                                                                                    cTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGCGGTTTATCACGACAGCCTGTGCCGACAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGAAGACGGCTGGGAGCAACATATTCCCGCCCGCTATCTGGTCGGTGcc                                             >  2:790443/1‑251 (MQ=255)
                                                                                                                                                                                                               aTTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTTGTCGGTGCCGATGGCGCaa                                   >  1:358284/1‑250 (MQ=255)
                                                                                                                                                                                                                    cAGCGTCAAAACCTTCGACGTCGCCGAATCGCTAACGCGTA‑ACTACCAGCGAAGTTATATCAT‑TATTAATCGCCTTGCTTTCGACTTGTGGATGAATTCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTGCCGCCCGCGATCTGGTCGGTGCCGATGGCGCAAACTCg                              <  1:762356/250‑1 (MQ=255)
                                                                                                                                                                                                                      gCGTCAAACCCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATg                           <  1:203618/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGACAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCATACGGCTGGGAGCAACATATTACCGCCCGCAATCTGGTCGGTGCCCATGGCGCAAACTCGATgg                          >  2:809042/1‑251 (MQ=255)
                                                                                                                                                                                                                         acAAAACCGGCGCCGTCGCCGAAACGATAACGGGTA‑ACTACAAGCGAAGCTATATCAA‑TATTAACCGCCATTCTTGCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGACTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTg                        <  2:273085/250‑1 (MQ=255)
                                                                                                                                                                                                                          cAAACCCGTCGCCGTCGCCACATCGCTACCGCGTA‑ACTTCCTGCGAAGCTATATCAA‑TATTAATCGCCTTGCTTACGACTTGTGGATGACATCGTTGATCGCCGCCATCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTg                        <  1:436363/250‑1 (MQ=255)
                                                                                                                                                                                                                             aaCCGTCGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCACGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTgcggc                    >  1:94637/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTCGACGTCCCCGCATCGCTAACGCTTA‑ACTACCAGCGAAGCTATAACAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACa                <  2:906916/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gcGACGTCGCCGCATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACAt               <  1:481636/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cgccgcATCGCTAACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGATTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTAc         <  1:352052/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cgcATCGCTAACGCGTA‑ACCACCAGCGAAGCTATATCAA‑TATTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTCCCCgg     <  1:785493/251‑1 (MQ=255)
                                                                                                                                                                                                                                              caacGCTAACGCGAA‑ACTCCCAGCGAAGCTCTTACAA‑TATTACTCTCCATGCTTCCGACTTGTGGAGGTAATCGCTGAGCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGat   <  1:520342/248‑1 (MQ=255)
                                                                                                                                                                                                                                               agcGCTAACGCGTA‑ACTACCAGCGAAGCTATACCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGCGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGatc  <  1:32687/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTA‑ACTACCAGCGAAGCTATATCAA‑TATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATC  >  NC_000913/3869392‑3869881

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑