Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,938 T→A L454L (CTA→CTT)  sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,9380TA93.5% 199.4 / 5.0 64L454L (CTA→CTT) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base A (30/28);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTAC  >  NC_000913/75695‑76184
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
ggTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCGCCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑Gcc                                                                                                                                                                                                                                                    >  1:205906/1‑248 (MQ=255)
          cgtaGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGACGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGatagcgctctc                                                                                                                                                                                                                                       >  2:536089/4‑243 (MQ=255)
            cAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGa                                                                                                                                                                                                                                     <  2:431276/251‑1 (MQ=255)
             aGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCa                                                                                                                                                                                                                                                   >  1:143648/1‑236 (MQ=255)
             aGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCa                                                                                                                                                                                                                                                   <  2:143648/236‑1 (MQ=255)
                      aTACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCCGATATCGCTCTCGATCcctccttc                                                                                                                                                                                                                           >  2:551284/1‑249 (MQ=255)
                                gCCCCTGAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGa                                                                                                                                                                                                                 <  1:162673/251‑1 (MQ=255)
                                gCCCCGGAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTg                                                                                                                                                                                                                  <  1:709384/250‑1 (MQ=255)
                                      gAATGATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTTATTTCGCCAGTGCGCCAGCGCGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTa                                                                                                                                                                                                            <  1:113786/250‑1 (MQ=255)
                                          gATCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCCCTCCTTCATGCCACTGAGCGTAgcgg                                                                                                                                                                                                        <  2:152535/250‑4 (MQ=255)
                                           aTCAGCCAGTGGTG‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGGGTTGGAGCAGCGGCCCCTCGCACACGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGGTGGCGCTGTTAA‑GCCAGAAATCGCTCTAGATCTCTCCTTCCTGCCACTGATCGTAgctgat                                                                                                                                                                                                      >  2:519518/1‑246 (MQ=255)
                                                       gtggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGAttt                                                                                                                                                                                           <  1:17773/241‑1 (MQ=255)
                                                        tggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCTTTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTCA‑GCCAGATATCGCTCTCGATCTCTCCTTAATGCCACTGATCGTAGCTGTTCTCTTTGAtttc                                                                                                                                                                                          >  1:125343/10‑250 (MQ=255)
                                                                 ggCACCATTGCTTGGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGCTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGAGCGTAGCTGATCTCTTTGATTTCCAGCGTGAc                                                                                                                                                                                 <  2:121576/251‑1 (MQ=255)
                                                                  gCACCAACGCTTTGCTGGCGACCAGTTGCTGGCCCCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGGGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTGA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGAcc                                                                                                                                                                                <  2:249012/251‑1 (MQ=255)
                                                                   cACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCt                                                                                                                                                                               <  1:217257/251‑1 (MQ=255)
                                                                    aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGCTCGGAATCCAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGAtctc                                                                                                                                                                                                                                 >  2:406102/1‑200 (MQ=255)
                                                                    aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCGCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGAtctc                                                                                                                                                                                                                                 <  1:406102/200‑1 (MQ=255)
                                                                       attgCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATAGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTg                                                                                                                                                                           >  2:513627/4‑251 (MQ=255)
                                                                                cTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGAGGCACCTCGCACAAGTGCGCGAACAGCGCAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATAACGCTCTCGATCCCTCCTTCATGCCACTGATCGTAGCTAATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGcccg                                                                                                                                                                  >  2:541533/1‑249 (MQ=255)
                                                                                        cAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATctgct                                                                                                                                                          <  1:519518/251‑1 (MQ=255)
                                                                                            tGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGCTTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                                                                                                                                      >  2:560336/1‑251 (MQ=255)
                                                                                            ggcgggCACCAGTTCGCCAGGTTCATTTCGCCAGTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCCCCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                                                                                                                                      <  2:635438/247‑1 (MQ=255)
                                                                                                         tcgccagggtcattgcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCTCCTCGCACAAGCGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAAt                                                                                                                                          <  2:246722/235‑1 (MQ=255)
                                                                                                             cgggttcattccGCCATTGCGCCAGCGGGCGGCGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAAGCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCAccc                                                                                                                                     <  2:277345/240‑1 (MQ=255)
                                                                                                              aggttCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCg                                                                                                                                    >  1:350258/4‑251 (MQ=255)
                                                                                                                    atttcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTc                                                                                                                              >  1:703806/4‑251 (MQ=255)
                                                                                                                        cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTccatg                                                                                                                          >  2:143041/1‑247 (MQ=255)
                                                                                                                        cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACtg                                                                                                                          >  1:58436/1‑251 (MQ=255)
                                                                                                                            aTTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGAt                                                                                                                      <  1:317737/251‑1 (MQ=255)
                                                                                                                              tGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATc                                                                                                                     <  1:455182/250‑1 (MQ=255)
                                                                                                                              ggcgcCAGCGGGCGGGGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAATTTGGCGCTGGTAA‑GCCAGATATCGCTCTCTATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCTCCCGATGCTCACTGTGATcc                                                                                                                    <  2:13816/250‑1 (MQ=255)
                                                                                                                                cgcCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCTGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTg                                                                                                                  <  1:210514/251‑1 (MQ=255)
                                                                                                                                       gggcggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                             >  1:36657/6‑249 (MQ=255)
                                                                                                                                       gggcggcGTCGGCTTGCCAGTCGATCGGAATACAGGGTGGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCACTCTCGATCTCTCCTTCAGGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGCTaaagg                                                                                                           >  1:181635/6‑249 (MQ=255)
                                                                                                                                        ggcggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGGGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAgg                                                                                                          <  2:667825/247‑1 (MQ=255)
                                                                                                                                           ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                             >  2:170534/2‑245 (MQ=255)
                                                                                                                                           ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                             <  1:170534/244‑1 (MQ=255)
                                                                                                                                           ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGtc                                                                                                        >  1:301539/2‑249 (MQ=255)
                                                                                                                                              ggcGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGGTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACGGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATTCTCACTGTGATCCTGATAAAAGGTCAgga                                                                                                    <  2:620400/249‑2 (MQ=255)
                                                                                                                                              gTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAgga                                                                                                    <  1:551284/251‑2 (MQ=255)
                                                                                                                                                       ccAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCttt                                                                                           >  1:702439/1‑251 (MQ=255)
                                                                                                                                                       ccAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCttt                                                                                           >  1:268914/1‑251 (MQ=255)
                                                                                                                                                        cagtcgatcggcatacagtgttggAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTACAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTACATCGTGACCTGGTGACTTGCCAGCATCTGCTGCATGATCCCGGCAATCACCCGCTGCTCACTGTGATCCTGATAAAAGGTCCGGATGAGGCTTTc                                                                                          >  2:518493/19‑251 (MQ=255)
                                                                                                                                                         agtcgatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTcc                                                                                         <  1:12260/243‑1 (MQ=255)
                                                                                                                                                            cgatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAgg                                                                                      >  1:494151/6‑251 (MQ=255)
                                                                                                                                                             gatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGCGGCTTTcccggg                                                                                     >  2:541974/5‑247 (MQ=255)
                                                                                                                                                                cggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCGTTCCAGGCCAg                                                                                  <  2:242798/250‑1 (MQ=255)
                                                                                                                                                                         gtgttggAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGCTAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGttttt                                                                         >  1:397623/2‑251 (MQ=255)
                                                                                                                                                                                            tCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGATGTTTA‑GCCAGATATCATTTTCGATCGCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCgg                                                        <  1:249174/249‑1 (MQ=255)
                                                                                                                                                                                                     tgcgcgaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGcc                                              >  2:249124/19‑250 (MQ=255)
                                                                                                                                                                                                        gcgaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGATGCCAACg                                          >  1:610071/16‑251 (MQ=255)
                                                                                                                                                                                                          gaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAC‑GCCAGATATCGCTCTCGATATCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGCTGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCACTTTATGGGGCGGGCCTGGTGCAAACggg                                        >  2:41020/14‑251 (MQ=255)
                                                                                                                                                                                                              agcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAg                                    <  1:9836/242‑1 (MQ=255)
                                                                                                                                                                                                                            cGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCg                      <  1:541119/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCggg                    <  1:41020/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCAGCGTAAAGTTGGCGC‑GTTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTTACCTGGTGACTTGCCAGAATCCGCTGCATGAGCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCggg                    <  2:449298/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cgTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCATTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGCCCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTACAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAAcc                >  2:609615/2‑251 (MQ=255)
                                                                                                                                                                                                                                    tAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACcac              >  1:121123/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGAATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGACCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGACCATAAGCCGGGAACcac              >  1:625164/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGACAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGc        >  1:26759/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gcgcTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGATCCACAGTTGCTg      <  1:577362/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gcgcTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTg      <  1:167223/251‑1 (MQ=255)
                                                                                                                                                                                                                                                tGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGACAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACTGTTGCTGGTAc  >  1:701329/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
GGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTAC  >  NC_000913/75695‑76184

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