Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,950 A→G S450S (AGT→AGC 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,9500AG92.4% 183.4 / 0.5 66S450S (AGT→AGCsgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (29/32);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

CAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGC  >  NC_000913/75707‑76195
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
cAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGa                                                                                                                                                                                                                                                <  2:431276/251‑1 (MQ=255)
          aTACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCCGATATCGCTCTCGATCcctccttc                                                                                                                                                                                                                                      >  2:551284/1‑249 (MQ=255)
                    gCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGa                                                                                                                                                                                                                            <  1:162673/251‑1 (MQ=255)
                    gCCCCGGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTg                                                                                                                                                                                                                             <  1:709384/250‑1 (MQ=255)
                          gAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTTATTTCGCCAGTGCGCCAGCGCGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTa                                                                                                                                                                                                                       <  1:113786/250‑1 (MQ=255)
                              gATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCCCTCCTTCATGCCACTGAGCGTAgcgg                                                                                                                                                                                                                   <  2:152535/250‑4 (MQ=255)
                               aTCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGGGTTGGAGCAGCGGCCCCTCGCACACGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGGTGGCGCTGTTAAGCCAGAAATCGCTCTAGATCTCTCCTTCCTGCCACTGATCGTAgctgat                                                                                                                                                                                                                 >  2:519518/1‑246 (MQ=255)
                                           gtggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGAttt                                                                                                                                                                                                      <  1:17773/241‑1 (MQ=255)
                                            tggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCTTTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTCAGCCAGATATCGCTCTCGATCTCTCCTTAATGCCACTGATCGTAGCTGTTCTCTTTGAtttc                                                                                                                                                                                                     >  1:125343/10‑250 (MQ=255)
                                                     ggCACCATTGCTTGGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGCTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGAGCGTAGCTGATCTCTTTGATTTCCAGCGTGAc                                                                                                                                                                                            <  2:121576/251‑1 (MQ=255)
                                                      gCACCAACGCTTTGCTGGCGACCAGTTGCTGGCCCCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGGGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTGAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGAcc                                                                                                                                                                                           <  2:249012/251‑1 (MQ=255)
                                                       cACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCt                                                                                                                                                                                          <  1:217257/251‑1 (MQ=255)
                                                        aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGCTCGGAATCCAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGAtctc                                                                                                                                                                                                                                            >  2:406102/1‑200 (MQ=255)
                                                        aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCGCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGAtctc                                                                                                                                                                                                                                            <  1:406102/200‑1 (MQ=255)
                                                           attgCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATAGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTg                                                                                                                                                                                      >  2:513627/4‑251 (MQ=255)
                                                                    cTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGAGGCACCTCGCACAAGTGCGCGAACAGCGCAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATAACGCTCTCGATCCCTCCTTCATGCCACTGATCGTAGCTAATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGcccg                                                                                                                                                                             >  2:541533/1‑249 (MQ=255)
                                                                            cAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATctgct                                                                                                                                                                     <  1:519518/251‑1 (MQ=255)
                                                                                tGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGCTTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                                                                                                                                                 >  2:560336/1‑251 (MQ=255)
                                                                                ggcgggCACCAGTTCGCCAGGTTCATTTCGCCAGTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCCCCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                                                                                                                                                 <  2:635438/247‑1 (MQ=255)
                                                                                             tcgccagggtcattgcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCTCCTCGCACAAGCGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAAt                                                                                                                                                     <  2:246722/235‑1 (MQ=255)
                                                                                                 cgggttcattccGCCATTGCGCCAGCGGGCGGCGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAAGCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCAccc                                                                                                                                                <  2:277345/240‑1 (MQ=255)
                                                                                                  aggttCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCg                                                                                                                                               >  1:350258/4‑251 (MQ=255)
                                                                                                        atttcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTc                                                                                                                                         >  1:703806/4‑251 (MQ=255)
                                                                                                            cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTccatg                                                                                                                                     >  2:143041/1‑247 (MQ=255)
                                                                                                            cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACtg                                                                                                                                     >  1:58436/1‑251 (MQ=255)
                                                                                                                aTTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGAt                                                                                                                                 <  1:317737/251‑1 (MQ=255)
                                                                                                                  tGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATc                                                                                                                                <  1:455182/250‑1 (MQ=255)
                                                                                                                  ggcgcCAGCGGGCGGGGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAATTTGGCGCTGGTAAGCCAGATATCGCTCTCTATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCTCCCGATGCTCACTGTGATcc                                                                                                                               <  2:13816/250‑1 (MQ=255)
                                                                                                                    cgcCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCTGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTg                                                                                                                             <  1:210514/251‑1 (MQ=255)
                                                                                                                           gggcggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                        >  1:36657/6‑249 (MQ=255)
                                                                                                                           gggcggcGTCGGCTTGCCAGTCGATCGGAATACAGGGTGGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCACTCTCGATCTCTCCTTCAGGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGCTaaagg                                                                                                                      >  1:181635/6‑249 (MQ=255)
                                                                                                                            ggcggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGGGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAgg                                                                                                                     <  2:667825/247‑1 (MQ=255)
                                                                                                                               ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                        <  1:170534/244‑1 (MQ=255)
                                                                                                                               ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                        >  2:170534/2‑245 (MQ=255)
                                                                                                                               ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGtc                                                                                                                   >  1:301539/2‑249 (MQ=255)
                                                                                                                                  ggcGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGGTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACGGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATTCTCACTGTGATCCTGATAAAAGGTCAgga                                                                                                               <  2:620400/249‑2 (MQ=255)
                                                                                                                                  gTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAgga                                                                                                               <  1:551284/251‑2 (MQ=255)
                                                                                                                                           ccAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCttt                                                                                                      >  1:702439/1‑251 (MQ=255)
                                                                                                                                           ccAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCttt                                                                                                      >  1:268914/1‑251 (MQ=255)
                                                                                                                                            cagtcgatcggcatacagtgttggAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTACAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTACATCGTGACCTGGTGACTTGCCAGCATCTGCTGCATGATCCCGGCAATCACCCGCTGCTCACTGTGATCCTGATAAAAGGTCCGGATGAGGCTTTc                                                                                                     >  2:518493/19‑251 (MQ=255)
                                                                                                                                             agtcgatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTcc                                                                                                    <  1:12260/243‑1 (MQ=255)
                                                                                                                                                cgatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAgg                                                                                                 >  1:494151/6‑251 (MQ=255)
                                                                                                                                                 gatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGCGGCTTTcccggg                                                                                                >  2:541974/5‑247 (MQ=255)
                                                                                                                                                    cggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCGTTCCAGGCCAg                                                                                             <  2:242798/250‑1 (MQ=255)
                                                                                                                                                             gtgttggAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGCTAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGttttt                                                                                    >  1:397623/2‑251 (MQ=255)
                                                                                                                                                                                tCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGATGTTTAGCCAGATATCATTTTCGATCGCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCgg                                                                   <  1:249174/249‑1 (MQ=255)
                                                                                                                                                                                         tgcgcgaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGcc                                                         >  2:249124/19‑250 (MQ=255)
                                                                                                                                                                                            gcgaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGATGCCAACg                                                     >  1:610071/16‑251 (MQ=255)
                                                                                                                                                                                              gaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTACGCCAGATATCGCTCTCGATATCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGCTGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCACTTTATGGGGCGGGCCTGGTGCAAACggg                                                   >  2:41020/14‑251 (MQ=255)
                                                                                                                                                                                                  agcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAg                                               <  1:9836/242‑1 (MQ=255)
                                                                                                                                                                                                                cGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCg                                 <  1:541119/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCggg                               <  1:41020/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCAGCGTAAAGTTGGCGCGTTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTTACCTGGTGACTTGCCAGAATCCGCTGCATGAGCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCggg                               <  2:449298/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCATTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGCCCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTACAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAAcc                           >  2:609615/2‑251 (MQ=255)
                                                                                                                                                                                                                        tAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACcac                         >  1:121123/1‑251 (MQ=255)
                                                                                                                                                                                                                        tAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGAATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGACCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGACCATAAGCCGGGAACcac                         >  1:625164/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGACAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGc                   >  1:26759/1‑251 (MQ=255)
                                                                                                                                                                                                                                gcgcTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGATCCACAGTTGCTg                 <  1:577362/251‑1 (MQ=255)
                                                                                                                                                                                                                                gcgcTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTg                 <  1:167223/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGACAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACTGTTGCTGGTAc             >  1:701329/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gttaaGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACt            >  2:595019/5‑251 (MQ=255)
                                                                                                                                                                                                                                        aaGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGt          <  1:560336/249‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGt          <  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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGC  >  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: