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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,869,6370CA95.3% 188.9 / ‑2.3 64L87L (CTC→CTAcbrAcolicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (34/27);  total (36/28)
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.27e-01

GCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCAT  >  NC_000913/3869395‑3869883
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
gCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGCTGGGCTGACGCTTCCTGTCGATGGGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACCTCGCCGCATCGCTCACGCGTa                                                                                                                                                                                                                                                 >  2:521662/1‑250 (MQ=255)
 cGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAAc                                                                                                                                                                                                                                               <  1:403022/251‑1 (MQ=255)
    ttattGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGAAGCTTCCTGTCGATGTGATCGCTAATTCGCAGATTTTTAGCGTCAAAACCGTCGACGTCGCCGCATCGCTATCGCGTAACTAc                                                                                                                                                                                                                                            >  1:305547/1‑251 (MQ=255)
        tGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCTCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAg                                                                                                                                                                                                                                         <  1:335082/250‑1 (MQ=255)
        tGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGc                                                                                                                                                                                                                                        >  1:269888/1‑251 (MQ=255)
                   cGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATc                                                                                                                                                                                                                             <  1:356285/251‑1 (MQ=255)
                    ggCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATc                                                                                                                                                                                                                             >  2:475687/1‑250 (MQ=255)
                       aGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCCGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGCAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACACTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTAAGCGTCAAAACCGACGACGTCGCCGAATCGCTAACGCGTAACTACCAGCGAAGCTATATCAat                                                                                                                                                                                                                          >  1:368291/1‑250 (MQ=255)
                            cGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGCAGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCCATATTAAt                                                                                                                                                                                                                    >  2:211181/1‑251 (MQ=255)
                             ggCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATc                                                                                                                                                                                                                   >  1:215681/1‑251 (MQ=255)
                             ggCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATc                                                                                                                                                                                                                   <  2:573755/251‑1 (MQ=255)
                                    gCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGc                                                                                                                                                                                                            <  1:168356/251‑1 (MQ=255)
                                     cACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCt                                                                                                                                                                                                           >  1:321966/1‑251 (MQ=255)
                                     aaCGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGTCGCTTCCTGTCGATGTGATAGCTAATCCGCAGAATTTCAGCGTCAAAACCAGCGACGTCGCCGCATAGGTAACGCGTAACTACAAGCAAATCTTTATACATATTAATCGcactgct                                                                                                                                                                                                           >  1:321969/2‑245 (MQ=255)
                                        gAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTc                                                                                                                                                                                                        <  2:265601/251‑1 (MQ=255)
                                                      aaaaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTa                                                                                                                                                                                                                                                 <  2:103326/196‑1 (MQ=255)
                                                      aaaaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTa                                                                                                                                                                                                                                                 >  1:103326/1‑196 (MQ=255)
                                                       aaaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGaaa                                                                                                                                                                                         <  1:51309/251‑1 (MQ=255)
                                                             aGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTg                                                                                                                                                                                   <  1:206716/251‑1 (MQ=255)
                                                                ggATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATc                                                                                                                                                                                <  2:169047/250‑1 (MQ=255)
                                                                    cgcgCTGGATAAAAAGCCCCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCTATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCTCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATcccc                                                                                                                                                                             <  2:42317/250‑1 (MQ=255)
                                                                         tGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTATGGGTGAAATCGCttaaccccccccgc                                                                                                                                                                       >  2:386222/1‑238 (MQ=255)
                                                                         ccctccacgccgctctcccgatctcTGAAGGTTCCAGCAACCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTCCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCACCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTCCGACTTGTGGATGAAATCGCTGACCCCCGCCAGc                                                                                                                                                                       <  2:73150/227‑1 (MQ=255)
                                                                          ggATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGCTGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCACTATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCg                                                                                                                                                                      >  1:104438/1‑251 (MQ=255)
                                                                                   gCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGCCTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGttt                                                                                                                                                             >  2:171198/1‑251 (MQ=255)
                                                                                     aacAGTGTAGTACTGAAGGGTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGCCCCGCCGCGGACTTTTATTCGCGGTGGGCTGACGCGTTCTGTCGGTGAGCTCGCTAATCAGCAGATTTTCAGCATCAAAACCGTCGACGTCGCTGCATCGCTAACGCGTAACTACCAGCACAGATATATCAAAATTAATCGCCATGCTATCGACTTGTGGAGGAAATCGATGCACCCCGCCACCGTTGAGGGTTAt                                                                                                                                                           >  2:98044/3‑251 (MQ=255)
                                                                                      ccAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACACGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGAGATCGCTCATCCGCAGATTTTCCGAGTCAACACCGTCGACGTCGCCGCATCGCTAACGCGTAACTAACAGCGGAGCTATATCAATATTAATCGCAATGCTTTCGCCTTGTGGATGAAATCGCTGATCCCCGCCCGCGTTGAGGTTTATc                                                                                                                                                          >  2:24841/1‑251 (MQ=255)
                                                                                         gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACg                                                                                                                                                       >  1:435719/1‑251 (MQ=255)
                                                                                          ggtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGa                                                                                                                                                      <  1:401004/250‑1 (MQ=255)
                                                                                           gtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCATCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGAc                                                                                                                                                     <  1:318152/251‑1 (MQ=255)
                                                                                              gTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCGATATCAATATTAATCGCCATGCGTTCGACTTGTGAATAAAATCGCTGATCCCCGCCAGCGTTTAGGTTTATCACGACAGc                                                                                                                                                  >  2:25543/1‑251 (MQ=255)
                                                                                                 cTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGc                                                                                                                                                                       >  1:73150/1‑227 (MQ=255)
                                                                                                       gTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATTCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCg                                                                                                                                          >  2:167812/1‑250 (MQ=255)
                                                                                                             ccaaccccTGCGGCGGTCGGCTGGACCCGGCCCCGCAGCGATCTTTTATTCGCGAGGGGCGGACGCTTCCTGTCGATGTGACCGCTAATCCGCAGATTTCCAGCGTCAAAACCTTCGACGTCGCCGCACCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATc                                                                                                                                   <  2:323955/246‑1 (MQ=255)
                                                                                                               aaaCCCGGCGGTGGTCTGCTGTCACGGGACGCGCAGCGATCTTTTATTCGCGAGGGGCTGACGCTTCCTGTCGATGTGATCGATAATCCGCAGATTTTCAGCGTCAAAACCGGCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGTTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTg                                                                                                                                 <  2:187662/251‑1 (MQ=255)
                                                                                                                aaCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGAACTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCACTATTAATCGCCATGCTTTCGGCTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGTGGTTTATCACGACAGCCTGGGCCGAAAGATCTgg                                                                                                                                >  1:155135/1‑251 (MQ=255)
                                                                                                                        ggcggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAgg                                                                                                                         >  1:10104/1‑250 (MQ=255)
                                                                                                                          ggggtggcttgCACCGGACGGCCAGCGATTTTTTAGCCGCGATGGGCTGCCGCTTCCTGTCGATGTGTGCGCTAATCCGCAGATTTTCAGCGCAAAAACCGTCGACGTCGCCGCTCCGCTAACGGGTAACTACCGGCGAAGCTATACCAATATTAATCGCCATGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatg                                                                                                                      <  2:442386/241‑1 (MQ=255)
                                                                                                                          cggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatg                                                                                                                      >  2:315090/1‑251 (MQ=255)
                                                                                                                           ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatg                                                                                                                      >  2:318717/1‑250 (MQ=255)
                                                                                                                              cggctgGCCCCGGACGCGCAGCGATCTTTTATTCGCGATGGGCGGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATa                                                                                                                   <  2:156433/248‑1 (MQ=255)
                                                                                                                               tgctgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATaa                                                                                                                  >  2:149135/1‑250 (MQ=255)
                                                                                                                                        cGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTc                                                                                                        <  1:56555/251‑1 (MQ=255)
                                                                                                                                            cgcgCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAttt                                                                                                    >  1:125890/1‑251 (MQ=255)
                                                                                                                                            cgcgCAGCGATCTTTTATCCGCGATGGGCTGACGTTTCCTGTCGATGTGAGCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAACCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAttt                                                                                                    <  1:493511/251‑1 (MQ=255)
                                                                                                                                                            aTTCGCGATGGGCTGCCGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATaa                                                                                                                  >  1:438240/1‑221 (MQ=255)
                                                                                                                                                            aTTCGCGATGGGCTGACGCTGCCTGCCGAGGGTATCGCAACGCCGTAGAGTTTCAGCCGAAAAACCGCCGACGGCCCCGCATCGTTAACGCGTAACGACCAGCGAAGCTATATCAATATCAATCGCCAGGCTGCCTACTTGTGGATGAAACCGCGGAGCCCCGCCAGCGTTGAGGTTTATCCCGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATaa                                                                                                                  <  2:438240/221‑1 (MQ=255)
                                                                                                                                                                       gCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACat                                                                          >  1:384880/1‑250 (MQ=255)
                                                                                                                                                                           aCGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAc                                                                     >  2:549100/1‑251 (MQ=255)
                                                                                                                                                                                    gTCGATGTGCTCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATAGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGAATTGGGGCTGAAATCGCTGATCCCCGCCAGCGTTGCAGTTTATCACGACAGCCTGTGCCGCAAACTCTGGCGTGAGGTTGATAAAGGGCATGTACTTTTTCGGGCAGACGGCTGGGCGCAACATCTTACAGCCCGCt                                                             >  2:331266/1‑250 (MQ=255)
                                                                                                                                                                                        agtgaacctcTAACCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCAATATCAATATTAAGCGCCATGCTTTCGACTTGTGGATGAAATCGCGGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTg                                                        <  2:155135/242‑1 (MQ=255)
                                                                                                                                                                                          gtgAGCGCTACTCCACAGATTTGCAGCGTCAAACCCTTCAACGTCCCCGCAGCGCTAACGCGTAACTACCAGCGAAGCTATATAAATATTAATGGCCATGCTTTCGACTGGGGGATGAAATCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGt                                                      <  2:215681/251‑1 (MQ=255)
                                                                                                                                                                                                    aTCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAACCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGAt                                            <  1:235631/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGCAGATTTTCAGCGTCAAATCCGTCGACGTCGCCGCATCGCTAACGCGTATCTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATg                                           <  1:211181/250‑1 (MQ=255)
                                                                                                                                                                                                              tttCAGCGCCAAAACCGCCGACGCCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATACCAAAATTAATCGCCATGCTTTCGACTTGGGGATGAAACCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAAc                                   <  2:269888/250‑1 (MQ=255)
                                                                                                                                                                                                                  aGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGAt                              >  1:427141/1‑251 (MQ=255)
                                                                                                                                                                                                                     gTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGt                           <  1:167812/251‑1 (MQ=255)
                                                                                                                                                                                                                              gTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGAc                   >  1:469508/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCCATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACAtct               >  2:74592/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACAGCCCGCTATCTGGTCGGTGCCGCTGGCGCAAACTCGAGGGTGCGGCGACATctct             >  2:253439/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gTCGCCGCAGCGCAAACGCGTAACTACCAGCAAAGTGATAGCAATATTAACCGCCATGCTTTCGACTTGTGGAGGAAATCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTa            <  1:407117/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gcATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTTGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCgg       >  1:481151/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aTCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGatc    >  1:471250/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGatcat  >  2:36327/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCAT  >  NC_000913/3869395‑3869883

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