Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 263,353 T→A 5.8% Q106L (CAA→CTA)  ykfI ← CP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913263,3530TA 5.8% 399.6 / 8.7 130Q106L (CAA→CTA) ykfICP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (60/62);  new base A (2/6);  total (62/68)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

TGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTG  >  NC_000913/263204‑263499
                                                                                                                                                     |                                                                                                                                                  
 ggTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTg                                                                                                                                                   <  2:799991/150‑1 (MQ=255)
    cgCATGATTGTTTGACCGTAAAACGATACCGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCg                                                                                                                                                <  2:309422/149‑1 (MQ=255)
     gCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGc                                                                                                                                               <  2:521355/150‑1 (MQ=255)
     gCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGc                                                                                                                                               <  2:703291/150‑1 (MQ=255)
      cATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGca                                                                                                                                              <  2:947505/150‑1 (MQ=255)
      cATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGca                                                                                                                                              >  2:1238183/1‑150 (MQ=255)
        tGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaac                                                                                                                                            >  2:199507/1‑150 (MQ=255)
          aTTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaacaa                                                                                                                                          >  2:1637518/1‑150 (MQ=255)
           ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGTGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaac                                                                                                                                         <  2:1310666/150‑2 (MQ=255)
           ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaac                                                                                                                                         <  1:1334701/150‑2 (MQ=255)
           ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaac                                                                                                                                         <  1:1404050/150‑2 (MQ=255)
           ttggtagtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTGCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaac                                                                                                                                         <  2:1660161/144‑2 (MQ=255)
             gtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGcc                                                                                                                                       <  2:563746/150‑1 (MQ=255)
               ttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAg                                                                                                                                     >  2:1826431/1‑150 (MQ=255)
                tgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaacgtagt                                                                                                                                    <  2:1151734/150‑7 (MQ=255)
                tgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGt                                                                                                                                    >  2:1426441/1‑150 (MQ=255)
                tgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGt                                                                                                                                    <  1:1351267/150‑1 (MQ=255)
                tgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGt                                                                                                                                    <  1:859212/150‑1 (MQ=255)
                 gtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGtt                                                                                                                                   >  2:670186/1‑150 (MQ=255)
                  tCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtg                                                                                                                                  <  1:1681475/150‑1 (MQ=255)
                  tCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtg                                                                                                                                  <  2:761556/150‑1 (MQ=255)
                   ccGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACTCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgc                                                                                                                                 <  2:1394530/150‑1 (MQ=255)
                   ccGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgc                                                                                                                                 >  2:42783/1‑150 (MQ=255)
                        aaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctgcc                                                                                                                            >  2:793290/1‑150 (MQ=255)
                        aaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctgcc                                                                                                                            >  1:1321073/1‑150 (MQ=255)
                        aaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATCTAAATCGGGAAGGTTGTACGCAATCTTCATCGTACTACGTTGTTACGGCTTTGCCGAAACAAGCCACTagactgcc                                                                                                                            >  2:1139464/1‑150 (MQ=255)
                           aGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCcgcg                                                                                                                         >  1:1271878/1‑150 (MQ=255)
                             aTAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCt                                                                                                                       <  1:121260/150‑1 (MQ=255)
                                aCGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGc                                                                                                                    >  2:683233/1‑150 (MQ=255)
                                     ccTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAAt                                                                                                               >  2:1749579/1‑150 (MQ=255)
                                     ccTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAAt                                                                                                               <  2:320026/150‑1 (MQ=255)
                                            ggTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCa                                                                                                        <  1:510203/150‑1 (MQ=255)
                                                  agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  2:899076/150‑1 (MQ=255)
                                                    cAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTATGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGgag                                                                                                >  1:1733183/1‑150 (MQ=255)
                                                     aGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGgaga                                                                                               <  2:1619797/150‑1 (MQ=255)
                                                     aGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGgaga                                                                                               <  2:1513946/150‑1 (MQ=255)
                                                     aGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGgaga                                                                                               <  2:594287/150‑1 (MQ=255)
                                                     aGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGGTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGgaga                                                                                               <  1:34902/150‑1 (MQ=255)
                                                      ggCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCCGCCCGGAGAt                                                                                              >  2:1731116/1‑150 (MQ=255)
                                                      ggCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGAt                                                                                              >  2:152957/1‑150 (MQ=255)
                                                      ggCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGAt                                                                                              >  2:1044091/1‑150 (MQ=255)
                                                         gcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAg                                                                                           <  1:1044091/150‑1 (MQ=255)
                                                          cATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAgg                                                                                          >  1:1419024/1‑150 (MQ=255)
                                                           aTTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGga                                                                                         <  1:1002784/150‑1 (MQ=255)
                                                            ttACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACACAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGAAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGgag                                                                                        >  2:455183/1‑150 (MQ=255)
                                                              aCGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAAGGTCTGCAGCCCGGAGATAAGGAGAt                                                                                      >  2:292885/1‑150 (MQ=255)
                                                              aCGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGAACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGAt                                                                                      >  1:914208/1‑150 (MQ=255)
                                                                gCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTg                                                                                    >  2:970684/1‑150 (MQ=255)
                                                                   aTAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTc                                                                                 <  1:33082/150‑1 (MQ=255)
                                                                      ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                                                              >  1:504407/1‑150 (MQ=255)
                                                                      ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                                                              >  2:198652/1‑150 (MQ=255)
                                                                          aacaacGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTGCATCGCACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAg                                                                          <  1:306219/150‑1 (MQ=255)
                                                                          aaCAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAg                                                                          >  2:542530/1‑150 (MQ=255)
                                                                            cAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCt                                                                        <  2:906960/150‑1 (MQ=255)
                                                                                   ggaAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATccc                                                                 >  1:569059/1‑150 (MQ=255)
                                                                                     aaGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCtt                                                               <  2:1746132/150‑1 (MQ=255)
                                                                                       gTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGGTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                             <  1:234712/150‑1 (MQ=255)
                                                                                         tCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTg                                                           <  1:122715/150‑1 (MQ=255)
                                                                                           aGAAATGTACATCGGGAAGGGTGTACGCAATGTTCAACGGACTACGTTGTTACGGATTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGACAGTGTTCCTGCCAGCTAAATCCCTTCCTGTc                                                         <  1:289982/150‑1 (MQ=255)
                                                                                           aGAAATGTAAATCGGTAATGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCATTGACGAAACAAGACAGTAGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGCAGATAAGGAGATTTTTCCTGCCAGCTAAATCCATACCTGTc                                                         >  1:1759815/1‑150 (MQ=255)
                                                                                             aaaTGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGa                                                       <  2:1635259/150‑1 (MQ=255)
                                                                                             aaaTGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGa                                                       <  1:1173926/150‑1 (MQ=255)
                                                                                              aaTGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGAt                                                      <  2:1429068/150‑1 (MQ=255)
                                                                                               aTGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATa                                                     >  1:1032292/1‑150 (MQ=255)
                                                                                                tGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATAc                                                    >  1:1485434/1‑150 (MQ=255)
                                                                                                tGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATAc                                                    <  2:1818544/150‑1 (MQ=255)
                                                                                                 gTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACg                                                   >  2:435799/1‑150 (MQ=255)
                                                                                                 gTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACg                                                   <  1:800253/150‑1 (MQ=255)
                                                                                                  tAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGACTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGa                                                  >  1:108266/1‑150 (MQ=255)
                                                                                                    aaTCGGGAAGGTTGTACGCAAAGTTCATCGTACTACGTTGTAACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATAGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAAc                                                >  1:1527454/1‑150 (MQ=255)
                                                                                                     aTCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGATCCTGCCAGCTAAATCCCTTCCTGTCGATACGAAcc                                               <  1:454913/150‑1 (MQ=255)
                                                                                                      tCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCa                                              >  1:691346/1‑150 (MQ=255)
                                                                                                         ggAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCt                                           >  1:1554974/1‑150 (MQ=255)
                                                                                                            aGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGt                                        <  2:1321073/150‑1 (MQ=255)
                                                                                                                  tACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAttttt                                  >  1:430175/1‑150 (MQ=255)
                                                                                                                   aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTc                                 <  1:1637518/150‑1 (MQ=255)
                                                                                                                    cGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCt                                <  1:983794/150‑1 (MQ=255)
                                                                                                                     gCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTa                               <  2:1518198/150‑1 (MQ=255)
                                                                                                                       aaTGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTAcc                             <  1:42783/150‑1 (MQ=255)
                                                                                                                        aTGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCa                            <  1:964911/150‑1 (MQ=255)
                                                                                                                         tGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAg                           >  1:462600/1‑150 (MQ=255)
                                                                                                                         tGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGCAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAg                           >  1:122099/1‑150 (MQ=255)
                                                                                                                          gTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGACCCAGCTCGTATTTTTCTACCAGa                          >  1:1002609/1‑150 (MQ=255)
                                                                                                                             cATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaaaa                       >  2:877727/1‑150 (MQ=255)
                                                                                                                                 gTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCa                   <  1:164236/150‑1 (MQ=255)
                                                                                                                                  tactacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCAc                  >  1:1764723/1‑150 (MQ=255)
                                                                                                                                    ctggcttgttgCGGCATAGCCGTAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGACCCAGCTCGTTTTTTTCTACCAGAAAATTCACgg                >  1:1090947/6‑150 (MQ=255)
                                                                                                                                    ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                <  1:1238183/150‑1 (MQ=255)
                                                                                                                                     tacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGc               <  2:1259454/150‑1 (MQ=255)
                                                                                                                                     tacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGc               <  2:850786/150‑1 (MQ=255)
                                                                                                                                       cGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCag             >  2:1690995/1‑149 (MQ=255)
                                                                                                                                       cGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCAt             >  2:1697879/1‑150 (MQ=255)
                                                                                                                                           gttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGc         <  2:7479/150‑1 (MQ=255)
                                                                                                                                             tACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCCGCTCGTATTTTTCTACCAGAAAATTCCCGGCCTCGGCTa       >  1:1150823/1‑150 (MQ=255)
                                                                                                                                             tACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTa       >  2:976007/1‑150 (MQ=255)
                                                                                                                                               cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCCCGGCATCGGCTAgg     >  1:1660674/1‑150 (MQ=255)
                                                                                                                                               cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAgg     >  1:1558186/1‑150 (MQ=255)
                                                                                                                                                ggCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGCACCAGCTCGTATTTTTCTACCAGCAAATTCACGGCGTCGGCTAggg    >  1:859650/1‑150 (MQ=255)
                                                                                                                                                  cTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTg  <  1:1311939/150‑1 (MQ=255)
                                                                                                                                                  cTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGGCTGCAGTCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTg  <  2:96613/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTG  >  NC_000913/263204‑263499

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