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,6370CA100.0% 170.0 / NA 48L87L (CTC→CTAcbrAcolicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/21);  total (27/21)

TAGTCTTTTGCCATCGTGATGAATGTGAGTGAAAAGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACAT  >  NC_000913/3869341‑3869931
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      
tAGTCTTTTGCCATCGTGATGAATGTGAGTGAAAAGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAAcg                                                                                                                                                                                                                                                                                                     <  2:7171/300‑1 (MQ=255)
            aTCGTGATGAATGTGAGTGAAAAGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCa                                                                                                                                                                                                                                                                                          >  2:346855/1‑299 (MQ=255)
              cGTGATGAATGTGAGTGAAAAGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGc                                                                                                                                                                                                                                                                                        >  1:476957/1‑299 (MQ=255)
               gTGATGAATGTGATTGAAAAGAGTGAACATTTCAACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCGAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCg                                                                                                                                                                                                                                                                                       <  2:253756/299‑1 (MQ=255)
                         tgagtgaAAAGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGNAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCa                                                                                                                                                                                                                                                                            <  1:446743/300‑1 (MQ=255)
                                  aGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATc                                                                                                                                                                                                                                                                   <  1:181172/300‑1 (MQ=255)
                                  aGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACAAGCGAAGCTATATCAATATTAATc                                                                                                                                                                                                                                                                   <  1:403921/300‑1 (MQ=255)
                                   gATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATc                                                                                                                                                                                                                                                                   >  2:549208/1‑299 (MQ=255)
                                                                                            aCGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCa                                                                                                                                                                                                         <  1:465932/300‑1 (MQ=255)
                                                                                                     aGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCt                                                                                                                                                                                                <  1:533772/300‑1 (MQ=255)
                                                                                                               aTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGaaaa                                                                                                                                                                                      >  1:352833/1‑300 (MQ=255)
                                                                                                                                ggATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatga                                                                                                                                                                     >  1:411452/1‑300 (MQ=255)
                                                                                                                                       aaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATgg                                                                                                                                                              <  1:387228/300‑1 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGTGATCTTTTATTCGCGATGGGCTGACGCTTACTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACNAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGAAAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAt                                                                                                                                                      >  2:184539/1‑300 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAt                                                                                                                                                      >  1:42876/1‑300 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAt                                                                                                                                                      >  2:130639/1‑300 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCCGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAt                                                                                                                                                      >  1:459831/1‑300 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCAATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATCAATGGCATGTCAt                                                                                                                                                      >  1:186238/1‑300 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCAGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAt                                                                                                                                                      >  1:187591/1‑300 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCGTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCAGCATCGCTAACGCGTAACTACCAGCGAAGCTATATTAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGAATGTCAt                                                                                                                                                      >  1:479922/1‑300 (MQ=255)
                                                                                                                                               gtgtGGTACTGAAGGTNTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCNAAAAATCTGGCGTGAGGANGATAAATGGCATGTc                                                                                                                                                        >  1:509505/1‑298 (MQ=255)
                                                                                                                                                                  aGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCNGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCt                                                                                                                                     >  2:485731/1‑298 (MQ=255)
                                                                                                                                                                      aaCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGc                                                                                                                               <  1:68568/300‑1 (MQ=255)
                                                                                                                                                                        cccTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCaa                                                                                                                             >  1:109081/1‑300 (MQ=255)
                                                                                                                                                                                 ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGGCGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGNTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTTAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAc                                                                                                                     <  2:442223/299‑1 (MQ=255)
                                                                                                                                                                                 ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCNCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAc                                                                                                                     <  1:132371/299‑1 (MQ=255)
                                                                                                                                                                                 ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGNGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAc                                                                                                                     <  2:476957/299‑1 (MQ=255)
                                                                                                                                                                                 ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGNTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAc                                                                                                                     <  2:442230/299‑1 (MQ=255)
                                                                                                                                                                                    ctgctgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGCCGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGCCAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTAGTGCAGACGGCTGGGAGCAACATATTAacgcc                                                                                                                 >  1:410212/1‑300 (MQ=255)
                                                                                                                                                                                    ctgctgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgcc                                                                                                                 >  1:87407/1‑300 (MQ=255)
                                                                                                                                                                                      gctgGCACCGGACGCGCGGCGATCTTTTATTCGCGATGTGCTTACGCTTCCTGTCGATTTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCATCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAc                                                                                                                     <  2:243038/294‑1 (MQ=255)
                                                                                                                                                                                      gctgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTNGCCGCATCGCTAACGCGTAACTACCAGCGAAGNTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgccc                                                                                                                >  2:457110/1‑299 (MQ=255)
                                                                                                                                                                                               ggACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGt                                                                                                      <  2:491138/300‑1 (MQ=255)
                                                                                                                                                                                                gACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTc                                                                                                     <  2:427220/300‑1 (MQ=255)
                                                                                                                                                                                                          gATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATg                                                                                           <  1:345876/300‑1 (MQ=255)
                                                                                                                                                                                                            tCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgc                                                                                         >  2:502606/1‑300 (MQ=255)
                                                                                                                                                                                                            tCTTTTATTCGCGATGGGCTGACGCTCCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGTGAAGCTATATCAATATTAATCGACATGCTCTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAAGTTTATCACGACAGCCTGTGCAGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTAGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgc                                                                                         >  2:502603/1‑300 (MQ=255)
                                                                                                                                                                                                                           gggCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGTTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGNTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGt                                                                           >  2:197549/1‑299 (MQ=255)
                                                                                                                                                                                                                           gggCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTNCTTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGt                                                                           >  1:446140/1‑299 (MQ=255)
                                                                                                                                                                                                                           gggCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGt                                                                           >  1:192709/1‑299 (MQ=255)
                                                                                                                                                                                                                                 aCGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGa                                                                    <  1:5102/300‑1 (MQ=255)
                                                                                                                                                                                                                                    cTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACAt                                                                 >  1:150929/1‑300 (MQ=255)
                                                                                                                                                                                                                                                         aaTCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCAGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCAATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGCGGTTTATCACGCCAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATCAATGGCATGTCATTTTTCGTGCCGACGGCTGGGAGCAACATATTACCGCCCGCTATTTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATc                                            >  2:502935/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                 aCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCTAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCcagcag                    <  2:388003/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 aCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCcagcag                    <  1:540714/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    gTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGt                   >  1:28979/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                          gTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCTACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGTGGTTcgc            <  2:523589/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   tCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACAt  >  1:142108/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      
TAGTCTTTTGCCATCGTGATGAATGTGAGTGAAAAGATGGAACATTTCGACGTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCATCAAATCCGTAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACAT  >  NC_000913/3869341‑3869931

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