Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 263,316 T→C 5.2% intergenic (‑104/+12) ykfN ← / ← ykfI pseudogene, YdiA family/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,3160TC 5.2% 251.3 / 2.8 77intergenic (‑104/+12)ykfN/ykfIpseudogene, YdiA family/CP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base C (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTA  >  NC_000913/263168‑263462
                                                                                                                                                    |                                                                                                                                                  
aaTTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGtt                                                                                                                                                    <  1:443158/149‑1 (MQ=255)
 aTTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTg                                                                                                                                                   >  3:543723/1‑149 (MQ=255)
 aTTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTg                                                                                                                                                   >  1:295270/1‑149 (MQ=255)
  ttAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGt                                                                                                                                                  <  1:233493/149‑1 (MQ=255)
   tAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTa                                                                                                                                                 <  3:110477/149‑1 (MQ=255)
    aGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTAc                                                                                                                                                <  2:128032/149‑1 (MQ=255)
    aGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTAc                                                                                                                                                <  1:440541/149‑1 (MQ=255)
         aTTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAAt                                                                                                                                           >  2:444782/1‑149 (MQ=255)
             tGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTc                                                                                                                                       <  3:94289/149‑1 (MQ=255)
              gATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCa                                                                                                                                      <  1:184591/149‑1 (MQ=255)
                  ggTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGt                                                                                                                                  >  4:418804/1‑149 (MQ=255)
                   gTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGta                                                                                                                                 <  4:543723/149‑1 (MQ=255)
                     tAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGtact                                                                                                                               <  1:60880/149‑1 (MQ=255)
                      aaTACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGtacta                                                                                                                              >  1:189056/1‑149 (MQ=255)
                           ttCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACgttg                                                                                                                         >  4:238711/1‑149 (MQ=255)
                           ttCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACgttg                                                                                                                         >  2:4138/1‑149 (MQ=255)
                           ttCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACgttg                                                                                                                         >  3:505624/1‑149 (MQ=255)
                              aaCTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTa                                                                                                                      >  2:122092/1‑149 (MQ=255)
                               aCTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTAc                                                                                                                     >  1:140160/1‑149 (MQ=255)
                                 ttCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACgg                                                                                                                   >  1:184835/1‑149 (MQ=255)
                                 ttCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGAGCATTACGCAAGAGGTAAACAAGGGAGGACGTTCTGAGATGTAAATGGGGAAGGTTGTACGCGATGTTCATCGTACTACGTTGTTACgg                                                                                                                   >  3:497810/1‑149 (MQ=255)
                                       ttGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGc                                                                                                             <  1:139923/149‑1 (MQ=255)
                                           aTGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGca                                                                                                         >  3:28193/1‑149 (MQ=255)
                                            tGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaa                                                                                                        <  1:21343/149‑1 (MQ=255)
                                            tGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGCGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaa                                                                                                        >  1:218928/1‑149 (MQ=255)
                                               ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaacaa                                                                                                     <  3:380303/149‑1 (MQ=255)
                                               ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaacaa                                                                                                     <  1:4138/149‑1 (MQ=255)
                                               ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaacaa                                                                                                     <  2:393492/149‑1 (MQ=255)
                                                   ttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCa                                                                                                 <  2:415490/149‑1 (MQ=255)
                                                        cGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgc                                                                                            <  3:238711/149‑1 (MQ=255)
                                                        cGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgc                                                                                            >  2:473147/1‑149 (MQ=255)
                                                           aaaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTTAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctg                                                                                         >  3:301194/1‑149 (MQ=255)
                                                           aaaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctg                                                                                         >  4:234734/1‑149 (MQ=255)
                                                           aaaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctg                                                                                         >  4:170844/1‑149 (MQ=255)
                                                           aaaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctg                                                                                         >  1:407062/1‑149 (MQ=255)
                                                               aGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCcgc                                                                                     >  2:179183/1‑149 (MQ=255)
                                                               aGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCcgc                                                                                     >  3:303473/1‑149 (MQ=255)
                                                                      gcgcCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCa                                                                              >  4:217326/1‑149 (MQ=255)
                                                                              cGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCt                                                                      <  4:526615/149‑1 (MQ=255)
                                                                              cGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCt                                                                      <  1:122092/149‑1 (MQ=255)
                                                                                 gtagtaGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCa                                                                   >  3:250563/1‑149 (MQ=255)
                                                                                      agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCg                                                              >  1:145884/1‑149 (MQ=255)
                                                                                        cAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGga                                                            <  2:235788/149‑1 (MQ=255)
                                                                                        cAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGga                                                            <  2:134221/149‑1 (MQ=255)
                                                                                        cAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGga                                                            <  1:47665/149‑1 (MQ=255)
                                                                                         aGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGgag                                                           >  3:1651/1‑149 (MQ=255)
                                                                                           gcgcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGCTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCTGGGAGAt                                                         <  2:537609/149‑1 (MQ=255)
                                                                                             gcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                       >  2:25533/1‑149 (MQ=255)
                                                                                             gcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                       <  4:244198/149‑1 (MQ=255)
                                                                                              cATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGt                                                                                               >  3:201168/1‑108 (MQ=255)
                                                                                              cATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGt                                                                                               <  4:201168/108‑1 (MQ=255)
                                                                                                ttACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTCGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                       <  1:549906/146‑1 (MQ=255)
                                                                                                ttACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTCGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                       >  2:549906/1‑146 (MQ=255)
                                                                                                ttACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTCGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                       <  1:189647/146‑1 (MQ=255)
                                                                                                ttACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTCGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                       >  2:189647/1‑146 (MQ=255)
                                                                                                 tACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGgag                                                   <  4:279641/149‑1 (MQ=255)
                                                                                                  aCGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGgaga                                                  <  1:191579/149‑1 (MQ=255)
                                                                                                          ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCt                                          <  2:189056/149‑1 (MQ=255)
                                                                                                           gTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                         <  1:250639/149‑1 (MQ=255)
                                                                                                              aaCAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCa                                      >  1:425901/1‑149 (MQ=255)
                                                                                                                   gggAGGAAGTTCAGAAATGTAAGTCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTGGTTACGGCTTTGCCGGAACAAGCCAGTTGCCTGCCGCGCTGGCAGAATGTCTGCAGCCCGGAGATAGGGGGATTGTTCCTGCCAGCTaa                                 >  2:271957/1‑149 (MQ=255)
                                                                                                                   gggAGGAAGTTCACAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTaa                                 >  4:400939/1‑149 (MQ=255)
                                                                                                                     gaggaAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGTCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAAt                               >  4:335754/1‑149 (MQ=255)
                                                                                                                        gaAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATccc                            >  4:100409/1‑149 (MQ=255)
                                                                                                                          aGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCtt                          <  4:303473/149‑1 (MQ=255)
                                                                                                                            ttCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCt                           <  2:52699/146‑1 (MQ=255)
                                                                                                                            ttCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCt                           >  1:52699/1‑146 (MQ=255)
                                                                                                                            ttCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                        >  2:153797/1‑149 (MQ=255)
                                                                                                                               aGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGt                     <  4:1651/149‑1 (MQ=255)
                                                                                                                               aGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGt                     >  2:561281/1‑149 (MQ=255)
                                                                                                                                gAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTc                    <  4:329066/149‑1 (MQ=255)
                                                                                                                                  aaTGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGa                  <  4:151920/149‑1 (MQ=255)
                                                                                                                                  aaTGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGa                  <  4:497810/149‑1 (MQ=255)
                                                                                                                                  aaTGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGa                  <  4:505624/149‑1 (MQ=255)
                                                                                                                                        aaTCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGaa            <  2:371962/149‑1 (MQ=255)
                                                                                                                                             ggAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGTAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGc       >  1:118329/1‑149 (MQ=255)
                                                                                                                                             ggAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGc       >  1:388411/1‑149 (MQ=255)
                                                                                                                                                  gTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTa  >  4:193361/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTA  >  NC_000913/263168‑263462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: