Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,953 T→C E449E (GAA→GAG)  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,9530TC93.9% 190.1 / 2.9 66E449E (GAA→GAG) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (30/32);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

CAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGC  >  NC_000913/75707‑76195
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
cAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGa                                                                                                                                                                                                                                                <  2:431276/251‑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                                                                                                                        <  1:170534/244‑1 (MQ=255)
                                                                                                                               ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                        >  2:170534/2‑245 (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)
                                                                                                                                                                                                                                     gttaa‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACt            >  2:595019/5‑251 (MQ=255)
                                                                                                                                                                                                                                        aa‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGt          <  1:560336/249‑1 (MQ=255)
                                                                                                                                                                                                                                        aa‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGt          <  1:513627/249‑1 (MQ=255)
                                                                                                                                                                                                                                               gATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAg   <  1:550940/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCATTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTATTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCGGGTACTGTTCCta     >  1:583831/1‑248 (MQ=255)
                                                                                                                                                                                                                                                atatCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGc  <  1:333883/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
CAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGC  >  NC_000913/75707‑76195

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: ‑