Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 22,373 T→C intergenic (+25/‑18) ribF → / → ileS bifunctional riboflavin kinase/FMN adenylyltransferase/isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091322,3730TC95.2% 173.5 / ‑0.7 63intergenic (+25/‑18)ribF/ileSbifunctional riboflavin kinase/FMN adenylyltransferase/isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (37/23);  total (38/25)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

CAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGAC  >  NC_000913/22128‑22618
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
cAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACgg                                                                                                                                                                                                                                                   >  1:360436/1‑250 (MQ=255)
   aCATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAAcc                                                                                                                                                                                                                                               <  1:635632/251‑1 (MQ=255)
    cATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCCGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGTCTGTTATGTAATCAAACCGAAACACGGAACCg                                                                                                                                                                                                                                              >  2:92727/1‑251 (MQ=255)
             cagCCCAACTGTTGCCGGTATTGGCCAGCCGCGGGAAGTGCTTTTGTTAGATGTGCCAATGGACCTTTCCGGTCGTCATATGCCAGTAGTGCTTCGTAAAAAAAAACGGAATGAGCGGCGATGTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAAAAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAGCCGAGAATCtga                                                                                                                                                                                                                                     <  2:167935/249‑1 (MQ=255)
             acGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCtg                                                                                                                                                                                                                                      >  1:181992/1‑250 (MQ=255)
                cccAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGCTGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCtgatg                                                                                                                                                                                                                                   >  1:530647/1‑250 (MQ=255)
                       gTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACt                                                                                                                                                                                                                            >  2:328465/1‑250 (MQ=255)
                          gCCGGTAGTCGCCAGCAGCTGGAAGTGCATTTGTTAGATTTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGTACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATaa                                                                                                                                                                                                                        <  2:45572/251‑1 (MQ=255)
                                   cGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAAccc                                                                                                                                                                                                                <  1:665989/250‑1 (MQ=255)
                                          agcTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACCAGTAGTGCTGCGTAAAAAAATACGCAATGGGCAGCGATTTGCGTCGCTGGACGAAATGAAAGAGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGATAACAAATCCGGCTTAAGCCTGTTATATAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTg                                                                                                                                                                                                        >  2:548111/1‑251 (MQ=255)
                                                aaGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGAGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCCAACCGAAACACGAAACCGCGAATCTGATGAGTGACTATAAATCAAACCTGAATTTGCCGGa                                                                                                                                                                                                   >  2:676752/1‑250 (MQ=255)
                                                              gctGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGCCGAACTGAAAGCGCAGATTGCGCGTGCTGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCACACCGAAACACGGACACGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCAGGAAACAGGGTTCCCGAt                                                                                                                                                                                    >  1:676612/3‑251 (MQ=255)
                                                              gATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGAt                                                                                                                                                                                    >  1:266758/1‑251 (MQ=255)
                                                               aTGTTGCAATGGGCCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATg                                                                                                                                                                                   <  2:266758/251‑1 (MQ=255)
                                                                   tGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAACAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAACCACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGcttg                                                                                                                                                                               >  1:516481/1‑248 (MQ=255)
                                                                               ttACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGTGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGcc                                                                                                                                                                    >  1:192234/1‑250 (MQ=255)
                                                                               ttACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGcc                                                                                                                                                                    <  1:328465/250‑1 (MQ=255)
                                                                                            ataCAAGTAGTGCTGAGTAAAAAAATACGAAATGAGCAGCGATTAGCGTCGCTGGACGACCTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGATAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGCATCTGATGAGAGACTATAAATCAACCCTGAATTTGCCGGAAAACGGGTTCCCGATGCGTGGCGATTCCGCCCAGCGCGCAccggg                                                                                                                                                      >  2:672895/1‑248 (MQ=255)
                                                                                               cAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGACTTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGAtct                                                                                                                                                                        >  2:597328/1‑230 (MQ=255)
                                                                                               cAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGAtct                                                                                                                                                                        <  1:597328/230‑1 (MQ=255)
                                                                                                 agtagaGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGc                                                                                                                                                 >  2:435486/1‑251 (MQ=255)
                                                                                                   tagtGCTGCGTAAAAAAATACGCCATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGAATATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTg                                                                                                                                               >  1:402199/1‑251 (MQ=255)
                                                                                                   tagtGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCACGCGCGAACCGGGAATGCTg                                                                                                                                               >  1:297583/1‑251 (MQ=255)
                                                                                                              aaaaaaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGa                                                                                                                                     <  2:181992/250‑1 (MQ=255)
                                                                                                                 aaaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtg                                                                                                                                  >  1:342798/1‑250 (MQ=255)
                                                                                                                  aaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgat                                                                                                                                <  2:676968/251‑1 (MQ=255)
                                                                                                                         cAATGCGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGACTTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCACGCGCGAACCGGGAATGCTGGCGCGTTGGACTGCTGTTGATc                                                                                                                         >  2:307024/1‑251 (MQ=255)
                                                                                                                                             cGCTGGACGAACTGAACGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGATTAAGCCTGTAATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGATAGACTATAAATCAACCCTGCATTTGCCGGAAACAGAGTTCCAGCTGCTATGCGATCACGCCAAGCGCGAACCGGGAATGCTGGCACGTTGGCATGATGATAATCTTTACGACATCATCCATGcc                                                                                                     >  1:273948/1‑250 (MQ=255)
                                                                                                                                               cTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGc                                                                                                   >  1:522496/1‑251 (MQ=255)
                                                                                                                                                tGGACGAACTGAACGCGAAGATTGCGCGTGATGATTTCACCGCCCGCGAATTTTTTGGGATAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGTAACACGGCACCGAGACTCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGACACAGGGTTCCCGCTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGCTGATCTGTACGGCATCATCCGAGCggat                                                                                                  >  2:276331/1‑249 (MQ=255)
                                                                                                                                                 ggACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGc                                                                                                   <  1:669245/249‑1 (MQ=255)
                                                                                                                                                       aCTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTGTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACGCGGAACCGAGAATCGGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTGCCCGAGGCGTGGCGATCGGGCCGAGGGCGAACCGGGAATGCTGGCGCGTGGTACGGATGATGATCTGTACGGCATCATCCGTGGGGCTAAAAAAg                                                                                           <  1:389366/251‑1 (MQ=255)
                                                                                                                                                          gAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAACCCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAACAAAGGc                                                                                         >  1:240636/1‑250 (MQ=255)
                                                                                                                                                          gAAAGCGCAGATTGCGCGTGATGAATTAACCGCCAGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGCTGCGTGGCGATCTCGCCAAGCGAGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGGACGGCATCATCCGTGCGGCTaaaaaggga                                                                                         >  2:116897/1‑246 (MQ=255)
                                                                                                                                                             aGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGTCGATCTCGCCAAGCGCGACCCGGGAATGCTGGCGCGTTGGGATGATGATGATCTGTACGGCATCATCCGGGCGGCTaacaacgccaaa                                                                                      >  2:459193/1‑240 (MQ=255)
                                                                                                                                                             aGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCaaaa                                                                                     <  2:342798/251‑1 (MQ=255)
                                                                                                                                                               cgcAGATTGCGCGTGATGAATTAACCGCCAGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGACTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGTTGATGATCTGTACGGCATCATACGAGCGGCTAAAAGAGGCaaaaa                                                                                    >  2:347187/1‑250 (MQ=255)
                                                                                                                                                                    aTTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCtt                                                                              >  1:572193/1‑251 (MQ=255)
                                                                                                                                                                       gcgcGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAAGCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCa                                                                            <  1:633930/250‑1 (MQ=255)
                                                                                                                                                                          cgtaatgaATTAACCGCCCGCGAATTTTTTGGGCTACCACCACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGACTTTGCCGGACCCAGGGTTCCCGCTCCGTGCCGATCTCGCCAAGCTCGAACCGGGAATGCTGTCGCGTTGGACTGATGCTGTTCTGTACGGCAGCATCCGTGCGGCTAAAAACGGCAAATAACCCTTCATTCt                                                                        >  2:25876/5‑251 (MQ=255)
                                                                                                                                                                           gTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATATCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCACACAAACCTTCATTCTg                                                                       >  1:311752/1‑251 (MQ=255)
                                                                                                                                                                                    tAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAAACGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGCCTGATGATGATCTGTACGACATCATCCGTGCGGCTAAAAAAGGCAAAAAACCCTTCATTCTGCATGATGGc                                                              >  2:666613/1‑251 (MQ=255)
                                                                                                                                                                                             gcgAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCACGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTAt                                                     >  2:120142/1‑251 (MQ=255)
                                                                                                                                                                                              cgAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGCACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCCTCATCCGTGCGGCTAAACAGGGCAAAAAAACCTTCATTCt                                                                        >  2:544895/1‑231 (MQ=255)
                                                                                                                                                                                              cgAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCt                                                                        <  1:544895/231‑1 (MQ=255)
                                                                                                                                                                                                     tttGGGCTTACAAAACCGGCTTAAGCCCTTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCGGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATg                                              <  2:673128/250‑1 (MQ=255)
                                                                                                                                                                                                         gactaaacaaccctgCTTAAGCCTGTTATGTAATCCAACCGAAACACGGAACCGAGAATCAGATGAGTGACTATAACTCAACCATGAATTTGCCGGAAACAGGCTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGTTGATCTGTACGGCTTCATCCGTACGGCTAAAAACGGCAAAAAAACCTTCATACTGCATGATGGCCCTCCTTATGCGACTGgaagc                                         >  2:65221/15‑251 (MQ=255)
                                                                                                                                                                                                          gCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAAACGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGCACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTACCTTCAGCATGATGGTCCTCCTTATGAGAATGgaagc                                         >  2:203341/1‑250 (MQ=255)
                                                                                                                                                                                                             aaCAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACAGAGAGTCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCCCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTACAAACGGCAGAAAAACCTTCATTCTGCATGATGACCCTCCTTATGCGACAGGCAGCAtt                                      >  2:338140/1‑250 (MQ=255)
                                                                                                                                                                                                              aCAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTc                                     >  1:105271/1‑250 (MQ=255)
                                                                                                                                                                                                              aCAAAACCGGCTTAAGCCTGTTATGTAATCAAACAGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                                                                                                                   <  1:172475/156‑1 (MQ=255)
                                                                                                                                                                                                              aCAAAACCGGCTTAAGCCTGTTATGTAATCAAACAGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                                                                                                                   >  2:172475/1‑156 (MQ=255)
                                                                                                                                                                                                                aaaaCCGGCTTAAGCCTGTTATGTAATCAAAACGAAACCCGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCat                                   >  1:183049/1‑250 (MQ=255)
                                                                                                                                                                                                                  aaCCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATAtt                                <  1:228810/251‑1 (MQ=255)
                                                                                                                                                                                                                   aCCGGCTTAAGCGTGTTATGAAATCAAACCGAAACACGGACCCGAGAATCTGCGGAGTGACTATAAATCAACCCTGAATTTGACAGAAACAGGGTTCCCGATGCGTGGCGATCTCGCAAAGCGCGAACCGGGAATGCTTGAGCATTGGACTGATGATGATCTGTACGGAATCATCCGTGCGGCTAAAAAAGGACACAAAACCTTCATTCTGCATCATGCACCACCTTAGACGAGTGGCAGAATTCaacttg                               >  2:73482/1‑246 (MQ=255)
                                                                                                                                                                                                                     cGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGt                             <  1:276331/251‑1 (MQ=255)
                                                                                                                                                                                                                          tAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACt                         <  1:363361/250‑1 (MQ=255)
                                                                                                                                                                                                                           aaGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTc                        <  1:335030/250‑1 (MQ=255)
                                                                                                                                                                                                                                   tATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGATTGACTATAAAGCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACa                <  2:294274/250‑1 (MQ=255)
                                                                                                                                                                                                                                       tAATCAACCCGAAACACGGAACCGTGAATCTGATGAGTGTCTATAAATCAACCCTGAATTTGCCGGAATCTGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGCATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGAt            <  2:145571/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGa       <  1:22285/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaaCCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGaa      <  2:102501/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTCAACAAGATTCTGAAAGAc  <  1:298835/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGAC  >  NC_000913/22128‑22618

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