Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,351 A→T N20Y (AAT→TAT)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,3510AT99.3% 514.1 / NA 137N20Y (AAT→TAT) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (66/69);  total (66/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

GCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGG  >  NC_000913/3937076‑3937625
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                  
gtgtttcctcctattaccagatgatcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATc                                                                                                                                                                                                                                                                          <  2:911067‑M2/123‑1 (MQ=255)
gtggtgcctcctattaacagatgatcgtcaaatctgtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATc                                                                                                                                                                                                                                                                          <  1:692429‑M2/123‑1 (MQ=255)
       ctcctattaccagatgatcgtcaaagcggtggtgattttgcaggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTcc                                                                                                                                                                                                                                                                    <  1:87740‑M2/129‑1 (MQ=255)
          ctagtaccaggtgtacggaaaagcgggggtgattttgctggcggggctggtaaacaacaaaaagcagtaataacgcctccaggccatcttgaatatgaacatgaaacaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGAAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACt                                                                                                                                                                                                                                                                <  2:20579‑M2/133‑1 (MQ=255)
            ctgtacagattatcgtcaaagctgtggtgattttgatggcggggctggtagacaacaaaaagcagtaataacgactacaggacatcgtgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTcc                                                                                                                                                                                                                                                              <  2:831753‑M2/135‑1 (MQ=255)
             ttaccagatgatcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgctaccttcagAAAGCAGTGCACGCGCCACTCTAACACGGTGGCGCATTTTATGGACATCCCGAATCTGCAAAACGCAGGCAGCTTCGTTGTTTTTGGCATACATAATGCTGACCACATTCTTTATCCTCAATCTTTTCCTActcac                                                                                                                                                                                                                                                             >  1:388614‑M2/151‑284 (MQ=255)
                 cagatgatcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCg                                                                                                                                                                                                                                                         <  2:164286‑M2/140‑1 (MQ=255)
                    atgatcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACCCGGCGGCGCCTTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACAACATTCTTTATCTTCAATCTTTTCCTACTCAAGGCgaac                                                                                                                                                                                                                                                      >  2:898182‑M2/144‑285 (MQ=255)
                    atgatcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGaaa                                                                                                                                                                                                                                                      <  2:89731‑M2/143‑1 (MQ=255)
                       atcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGACTATGCAAAACGCAGGCCGCCTCGTTGTTCTTGTCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTCCTCCAGGCCATAcc                                                                                                                                                                                                                                                    >  1:112316‑M2/141‑285 (MQ=255)
                       atcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAAcc                                                                                                                                                                                                                                                    >  1:701753‑M2/141‑285 (MQ=255)
                       atcgtcaaagcggtggtgattttgctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCg                                                                                                                                                                                                                                                   >  2:530415‑M2/141‑286 (MQ=255)
                               agcggtggggattttgatggcggtgctggtagaaaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagagccaccgtaagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCgg                                                                                                                                                                                                                                            <  2:847604‑M2/153‑1 (MQ=255)
                               agaggtggtgattttgctggcgtggctggtagagaagaaaaagcaggaataacgactaaaggacatgttgaatatgaacatgaaaagactgactaccctggtttccgctgttgagctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAAAATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACAACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCgg                                                                                                                                                                                                                                            <  1:713865‑M2/153‑1 (MQ=255)
                                  gtgggtgattttgctggcggttatggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAAc                                                                                                                                                                                                                                        <  2:753042‑M2/157‑1 (MQ=255)
                                        gattttgctgacgatactcgtagacaacaaaaatcagtaataacgcctacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccactgttgcgctaagcgccaccgtcagAAAGCAGGGCACACGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTa                                                                                                                                                                                                                                   <  2:388614‑M2/162‑1 (MQ=255)
                                              gctggcggtgctggtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGtt                                                                                                                                                                                                                            >  2:62254‑M2/118‑286 (MQ=255)
                                                           gtagacaacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCa                                                                                                                                                                                                               >  1:63207‑M2/105‑286 (MQ=255)
                                                                caacaaaaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGgc                                                                                                                                                                                                          <  1:530415‑M2/187‑1 (MQ=255)
                                                                      aaagcagtaataacgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAAt                                                                                                                                                                                                    <  2:645468‑M2/193‑1 (MQ=255)
                                                                           agtaataacgactaccggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGc                                                                                                                                                                                               <  2:63207‑M2/198‑1 (MQ=255)
                                                                                  acgactacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTCATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTACAGGCGAAACCGTAACCGGTAACCACTACCAGGTTGCCTTTGGCGGCAAAGGCGCGCATCAGGCtgtg                                                                                                                                                                                           >  2:347949‑M2/82‑283 (MQ=255)
                                                                                   cgactacaggacatcttgaatatgaacatgaaaaaactggctcccctggtttcctctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgctg                                                                                                                                                                                       <  1:790988‑M2/206‑1 (MQ=255)
                                                                                     actacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTaa                                                                                                                                                                                                                                         <  2:590460‑M2/156‑1 (MQ=255)
                                                                                     actacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTaa                                                                                                                                                                                                                                         >  1:590460‑M2/79‑234 (MQ=255)
                                                                                     actacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgctgct                                                                                                                                                                                     <  2:332355‑M2/208‑1 (MQ=255)
                                                                                       tacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATTTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTAACAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTgg                                                                                                                                                                                   >  1:467022‑M2/77‑286 (MQ=255)
                                                                                       tacaggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCCCTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTgg                                                                                                                                                                                   >  1:111819‑M2/77‑286 (MQ=255)
                                                                                         caggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGCCGCGAATCAGTCTGTGTCTGCtggga                                                                                                                                                                                 >  2:612711‑M2/75‑285 (MQ=255)
                                                                                         caggacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACAACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGc                                                                                                                                                                                 >  1:497208‑M2/75‑286 (MQ=255)
                                                                                          aggacatcttgaatatgaacctgaacaaactggctcccctgttttccgctgtgccgctaagcgccaccgtcagAAAGCAGGGCCCGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCg                                                                                                                                                                                <  1:73502‑M2/213‑1 (MQ=255)
                                                                                            gacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTa                                                                                                                                                                              >  1:381093‑M2/72‑286 (MQ=255)
                                                                                            gacatcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTa                                                                                                                                                                              >  1:337147‑M2/72‑286 (MQ=255)
                                                                                               atcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTa                                                                                                                                                                                                                                   <  2:349064‑M2/162‑1 (MQ=255)
                                                                                               atcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTa                                                                                                                                                                                                                                   >  1:349064‑M2/69‑230 (MQ=255)
                                                                                                tcttgaatattaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                            <  1:898182‑M2/217‑1 (MQ=255)
                                                                                                tcttgaatatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                            <  1:62254‑M2/217‑1 (MQ=255)
                                                                                                      atatgaacatgaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCAAGGCAAAACCGTAACCGGTATCCACTATCAGGTTGCATTTGGCGGAGAAGGCGCGAATCAGGCTGTGGCTGCTGGGTGTAGCGGTGAGaa                                                                                                                                                                    >  1:347173‑M2/62‑286 (MQ=255)
                                                                                                           aacatgaaaaaaccggctccccgggtttccgtcgtggcgataagcgacaacgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGAATTTTATGGAAAGCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTAGCTTCAATCTTTTCCTAATCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATcg                                                                                                                                                               <  2:165031‑M2/230‑1 (MQ=255)
                                                                                                                gaaaaaactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGttt                                                                                                                                                          <  1:57045‑M2/235‑1 (MQ=255)
                                                                                                                     aactggctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAGTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGc                                                                                                                                                     <  1:430923‑M2/240‑1 (MQ=255)
                                                                                                                          gctaccctggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTCGCTGCTGGGCGTTGCGGTCGGACTTTCGCCTTTTTTCCCTTTa                                                                                                                                                >  2:419296‑M2/42‑286 (MQ=255)
                                                                                                                               cctggtttccgcggttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGt                                                                                                                                           <  1:450438‑M2/250‑1 (MQ=255)
                                                                                                                               cctggtttcccctgttgcgctaagcgacacagtcagAAAGCAGGGCACGCGCCACCCTAGCACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGt                                                                                                                                           <  2:381474‑M2/250‑1 (MQ=255)
                                                                                                                                  ggtttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACCCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCTTTTTGCGGCAAAGCCGCGAATCAGGCTGTGCCTGCTGGGCGTAGCGGTGAGAATTACGTGTTTACTGCCTGTACGGGtga                                                                                                                                         >  2:201492‑M2/34‑285 (MQ=255)
                                                                                                                                    tttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGgc                                                                                                                                                                                                          >  2:672286‑M2/32‑218 (MQ=255)
                                                                                                                                    tttccgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGgc                                                                                                                                                                                                          <  1:672286‑M2/187‑1 (MQ=255)
                                                                                                                                        cgctgttgcgctaagcgccaccgtcagAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGc                                                                                                                                  <  1:40601‑M2/259‑1 (MQ=255)
                                                                                                                                                                        aGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCAc                                                                                                  <  1:214651/286‑1 (MQ=255)
                                                                                                                                                                          ggCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTAGTGCCTGTACGGGTGATGACAGCATGGGTGAGAGCGGTGGCAAGGAGCCGGCCACTg                                                                                                >  2:214651/1‑286 (MQ=255)
                                                                                                                                                                           gCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                            >  1:216561/1‑209 (MQ=255)
                                                                                                                                                                             ctgggcgaCCCTAACACGGTGTGGAATGTAATGGAATGCGCGAAAATGCAAAGCGCAGGAAGACTCGTTGTTCTTGGAAGCATTAATGCTGACCACATTCTTTAGCTTCAATATTTTCCTACACCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                              <  2:53573/278‑1 (MQ=255)
                                                                                                                                                                             ccgcgCCACCCTAGCGTGTTGGCGCTTTTTATGTACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                              <  2:30785/284‑1 (MQ=255)
                                                                                                                                                                             aCGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                              <  2:381093/285‑1 (MQ=255)
                                                                                                                                                                              cgcgCCACCCTAACACGGTGGCGCATTTTATGGACAGCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGa                                                                                               <  1:917331/283‑1 (MQ=255)
                                                                                                                                                                                cgcCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGTTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTTCGAATATCGCGTTTATTTCCTTTACGGGTTCTGGCCGCATTGGTGATAGCGTTCTCCAGCAGCTGCCCACTGATAACa                                                                                          >  1:554095/1‑286 (MQ=255)
                                                                                                                                                                                 ggcACCCTAACACTGTGGCGCATTTTATGGACATCCCGAATAGGCAAAACGCAGGCAGACTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAAGCTTTTCCTACGCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTGATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                         <  1:508202/284‑1 (MQ=255)
                                                                                                                                                                                  ccgcccTACCCCGGTGGCGCATTTTTTGGACTTCCCGAATATGCAAACCGCAGGCAGCCTCGTTTTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                         <  2:449894/282‑1 (MQ=255)
                                                                                                                                                                                    aCCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                            <  2:839618/200‑1 (MQ=255)
                                                                                                                                                                                    aCCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                            >  1:839618/1‑200 (MQ=255)
                                                                                                                                                                                    aCCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCAtt                                                                                                                               >  1:215924/1‑245 (MQ=255)
                                                                                                                                                                                    aCCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGTCaggaagttcgcgag                                                                                      >  2:261211/1‑273 (MQ=255)
                                                                                                                                                                                         aaCACGGTGGCTCATTTTATTGGCATCCCGAATATGCAAAACGCAGGCCGCCTCGTTTTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCt                                                                                                        <  2:916204/263‑1 (MQ=255)
                                                                                                                                                                                                tGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTTGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGTGGTGCGAATATCGCGTTTATTGCCTTTACGGGTGATGACAGCATTGGTGAGAGCGTTAGCCAGTATCGCGGCACTGATCACCTTGATGTTCCTCCgg                                                                           >  2:79131/1‑285 (MQ=255)
                                                                                                                                                                                                tGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATAtt                                                                                  >  1:449894/1‑278 (MQ=255)
                                                                                                                                                                                                tGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACt                                                                                                 >  1:916204/1‑263 (MQ=255)
                                                                                                                                                                                                tGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCgg                                                                           >  1:586260/1‑285 (MQ=255)
                                                                                                                                                                                                tGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCgg                                                                           >  1:851066/1‑285 (MQ=255)
                                                                                                                                                                                                            tGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGCCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTa                                                                                                                                                         >  1:359540/1‑195 (MQ=255)
                                                                                                                                                                                                                 aTCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGc                                                          >  1:845098/1‑285 (MQ=255)
                                                                                                                                                                                                                  tCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTTCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGa                                                        >  1:908675/1‑286 (MQ=255)
                                                                                                                                                                                                                  tCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCt                                                                                                        <  1:256921/238‑1 (MQ=255)
                                                                                                                                                                                                                  tCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATCTTACTCCGGTCAGCGTGATCAAAGgaga                                                        >  2:852694/1‑286 (MQ=255)
                                                                                                                                                                                                                  tCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATCGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCt                                                                                                        >  2:256921/1‑238 (MQ=255)
                                                                                                                                                                                                                   cccGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGc                                                                                                                                  >  1:11690/1‑211 (MQ=255)
                                                                                                                                                                                                                   cccGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGc                                                                                                                                  <  2:11690/211‑1 (MQ=255)
                                                                                                                                                                                                                   cccGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                         <  1:181219/252‑1 (MQ=255)
                                                                                                                                                                                                                           tGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                            >  1:61805/1‑161 (MQ=255)
                                                                                                                                                                                                                              aaaaCTCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACt                                                                               <  2:492318/251‑1 (MQ=255)
                                                                                                                                                                                                                              aaaaCGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACt                                                                               >  1:492318/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTg                                     >  1:144453/1‑286 (MQ=255)
                                                                                                                                                                                                                                             tCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAAACTTGCTATTACGCCGGTCAGCGTGATCAACGGCGCATCAACAGCTGTGGCGCTGAttttttt                             >  2:51694/1‑284 (MQ=255)
                                                                                                                                                                                                                                                   ttCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGata                                                                                    <  2:492637/225‑1 (MQ=255)
                                                                                                                                                                                                                                                   aagagacaggcaGTAATGATGACCAAATTCGTTAGCTTCAAGCTTTGCCGACGCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATg                        <  2:897473/276‑1 (MQ=255)
                                                                                                                                                                                                                                                       tGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACTGGTGTGGCGCTGATTTTTGTTAATGGCGaa                   >  1:285091/1‑286 (MQ=255)
                                                                                                                                                                                                                                                        ggCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgtg                                             >  1:62502/1‑259 (MQ=255)
                                                                                                                                                                                                                                                        ggCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgtg                                             <  2:62502/259‑1 (MQ=255)
                                                                                                                                                                                                                                                             caTTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGc                                                                                                                                                                                                                                                          >  2:706399/1‑49 (MQ=255)
                                                                                                                                                                                                                                                                       gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCt                                                                                                        >  1:679500/1‑185 (MQ=255)
                                                                                                                                                                                                                                                                       gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGa                                                                                               >  1:81728/1‑194 (MQ=255)
                                                                                                                                                                                                                                                                       gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGa                                                                                               <  2:81728/194‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAACGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAgg                 >  1:897473/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                       gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGc                                                          >  2:235595/1‑231 (MQ=255)
                                                                                                                                                                                                                                                                       gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGc                                                          <  1:235595/231‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCg   >  1:553149/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                        aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGgc                                                                                                                                                                                                          <  1:153612/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTa                                                                                                                                                <  1:465334/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTa                                                                                                                                                >  2:465334/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                        aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCgg  >  2:359231/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                             aTTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTAtt                                                                                                                                                       >  2:55827/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                  
GCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGG  >  NC_000913/3937076‑3937625

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