Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,066,228 G→T 6.5% R99L (CGG→CTG)  phoR → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,066,2280GT 6.5% 233.4 / 2.8 73R99L (CGG→CTG) phoRtwo‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base T (3/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTGCGCTTGCACGACTGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGACGGCGGCCAGCAAGTGACCAATCT  >  NC_002947/6066080‑6066369
                                                                                                                                                    |                                                                                                                                             
ctgcGCTTGCACGACTGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCgg                                                                                                                                              <  1:1329220/150‑1 (MQ=255)
      ttGCACGACTGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAg                                                                                                                                        <  1:665120/150‑1 (MQ=255)
        gCACGACTGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGa                                                                                                                                      >  1:236650/1‑150 (MQ=255)
         cACGACTGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCTGGTGCAGGAg                                                                                                                                     >  2:576378/1‑150 (MQ=255)
         cACGACTGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCTGGTGCAGGAg                                                                                                                                     >  2:576374/1‑150 (MQ=255)
               tGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCAcc                                                                                                                               >  2:376240/1‑150 (MQ=255)
                 gCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACcgc                                                                                                                             >  2:867267/1‑150 (MQ=255)
                   tGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCTGGTGCAGGAGTCCACcgctg                                                                                                                           <  1:576378/150‑1 (MQ=255)
                   tGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCTGGTGCAGGAGTCCACcgctg                                                                                                                           <  1:576374/150‑1 (MQ=255)
                      gcAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTCTTCGATAGCATCTACCACCGCCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCCGTGATCGACCGCGTGCAGGAGTCCACcgctgcgc                                                                                                                        >  1:1018957/1‑150 (MQ=255)
                          tCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCATGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTcc                                                                                                                                  >  1:206691/1‑136 (MQ=255)
                          tCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCATGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTcc                                                                                                                                  <  2:206691/136‑1 (MQ=255)
                           cACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGt                                                                                                                   <  1:782629/150‑1 (MQ=255)
                           cACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGt                                                                                                                   <  2:550683/150‑1 (MQ=255)
                            aCCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCgg                                                                                                                                              <  1:30770/122‑1 (MQ=255)
                            aCCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCgg                                                                                                                                              >  2:30770/1‑122 (MQ=255)
                            aCCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTg                                                                                                                  >  1:1043834/1‑150 (MQ=255)
                              cAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGAc                                                                                                                >  2:157330/1‑150 (MQ=255)
                                aCCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGAcgc                                                                                                              >  1:591929/1‑150 (MQ=255)
                                 cccGACGCAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCAACAACAACGCCGCGACTAAAGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGAcgcg                                                                                                             <  1:115436/150‑1 (MQ=255)
                                 cccGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGATGTATTCGATAGCATCTACCACCTGCCACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGAcgcg                                                                                                             >  1:528790/1‑150 (MQ=255)
                                 cccGACGAAGAACCGACCGATGGCTACTGCCGGGGGGGCGAGTTATTCGATAGCATCGACCACCTGCAACGCCGCGACCAGCGCGTCCGGGGCCGACTGCAGGCAGCGATCGACCTGGTGCAGGAGGCCACCGCTGCGCTGCGTGAcgcg                                                                                                             >  2:1177493/1‑150 (MQ=255)
                                  ccGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGGGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCgg                                                                                                            <  2:224439/150‑1 (MQ=255)
                                      cgaCGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGAt                                                                                                        <  1:70321/150‑1 (MQ=255)
                                      cgaAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCggg                                                                                                                                             <  1:754775/113‑1 (MQ=255)
                                      cgaAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCggg                                                                                                                                             >  2:754775/1‑113 (MQ=255)
                                           cACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGc                                                                                                   <  2:1043834/150‑1 (MQ=255)
                                            aCCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCt                                                                                                  >  2:620868/1‑150 (MQ=255)
                                            aCCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCt                                                                                                  >  2:927971/1‑150 (MQ=255)
                                             ccGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTc                                                                                                 <  2:1406173/150‑1 (MQ=255)
                                                cccGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGGGACGCGGTGATCATGGTCggc                                                                                              >  2:761250/1‑148 (MQ=255)
                                                   gATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGc                                                                                           >  2:706041/1‑150 (MQ=255)
                                                        ctatggCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACgg                                                                                      <  1:1381325/146‑1 (MQ=255)
                                                         tACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGc                                                                                     >  2:133527/1‑150 (MQ=255)
                                                           cgacgtgtgGGGCGAGGGATTCGATATCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGTCCAGCTGGAGGCAGGGATGCACCGGGGGCAGGAGTGCACCGCTGGGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCaa                                                                                   <  2:264954/145‑1 (MQ=255)
                                                            ggCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAAc                                                                                  <  1:96497/150‑1 (MQ=255)
                                                               cTGTGGGGCGATGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTg                                                                               <  2:528790/150‑1 (MQ=255)
                                                               cTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTg                                                                               <  2:591929/150‑1 (MQ=255)
                                                               cTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTg                                                                               >  1:1554525/1‑150 (MQ=255)
                                                                      gCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAAtggt                                                                        >  1:1633531/1‑150 (MQ=255)
                                                                      gCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAAtggt                                                                        <  1:157330/150‑1 (MQ=255)
                                                                      gCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCACCTGCAGGCAGTGATCGACCGGGCGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTTATCATGCTTGACAGCGACGGCAACCTGGAAtggt                                                                        <  2:598010/150‑1 (MQ=255)
                                                                         aGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGa                                                                     <  2:1633531/150‑1 (MQ=255)
                                                                          ggTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGaa                                                                    >  2:821591/1‑150 (MQ=255)
                                                                          ggTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGaa                                                                    <  2:1018957/150‑1 (MQ=255)
                                                                           gTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAAc                                                                   <  2:1343641/150‑1 (MQ=255)
                                                                              ttCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGGTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTAGAATGGTGGAACcgc                                                                >  1:1252291/1‑150 (MQ=255)
                                                                              ttCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCTGTGATCATGCTCGACAGCGACTGCAACCTGGAATGGTGGAACcgc                                                                <  2:236650/150‑1 (MQ=255)
                                                                              ttCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACcgc                                                                <  1:71489/150‑1 (MQ=255)
                                                                                     gCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAAc                                                                   >  1:234858/1‑140 (MQ=255)
                                                                                     gCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAAc                                                                   <  2:234858/140‑1 (MQ=255)
                                                                                      cATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCga                                                        <  1:867267/150‑1 (MQ=255)
                                                                                             caccTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTctg                                                 >  1:1284798/1‑150 (MQ=255)
                                                                                                  gCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACcgc                                                                >  2:1454613/1‑130 (MQ=255)
                                                                                                  gCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACcgc                                                                <  1:1454613/130‑1 (MQ=255)
                                                                                                  gCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTggg                                            <  1:133527/150‑1 (MQ=255)
                                                                                                       gccgcGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCa                                       >  2:487896/1‑150 (MQ=255)
                                                                                                          gcgACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGa                                    >  2:1049720/1‑150 (MQ=255)
                                                                                                             aCCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGTCGCCGAGACTCTGCTGGGCTTCAAGAccc                                 <  2:1554525/150‑1 (MQ=255)
                                                                                                             aCCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGAccc                                 <  2:1284798/150‑1 (MQ=255)
                                                                                                               cAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCca                               >  1:209224/1‑150 (MQ=255)
                                                                                                                 gcgcGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCcaca                             <  1:761250/150‑1 (MQ=255)
                                                                                                                        cGGGGCCGCCTGCAGGCAGTGATCGACCGGGCGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGAGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGAcggc                      >  2:747889/1‑150 (MQ=255)
                                                                                                                         ggggCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGAcggcg                     <  2:1124230/150‑1 (MQ=255)
                                                                                                                         ggggCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGAcggcg                     <  1:620868/150‑1 (MQ=255)
                                                                                                                            gccgccTGCAGGCAGTGATCGACCGGGTGCAGGAGTCAACCGCTGAGCTGCGAGACCGGGTGATCATGCTCGTCAACGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGACGGCGGcc                  >  1:199858/1‑150 (MQ=255)
                                                                                                                               gccTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGAc                                                                                        <  2:55279/77‑1 (MQ=255)
                                                                                                                               gccTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGAc                                                                                        >  1:55279/1‑77 (MQ=255)
                                                                                                                               gccTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGACGGCGGCCAGc               <  1:1049720/150‑1 (MQ=255)
                                                                                                                                 cTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGAccc                                 <  1:1448015/130‑1 (MQ=255)
                                                                                                                                 cTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGAccc                                 >  2:1448015/1‑130 (MQ=255)
                                                                                                                                      atctgtaaTCGACCGGGTGCAGGAGTCCACCGCTTCGCTGCGTAACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGCTGGAATCACGCCGCGGCGACTCCGCTGGGCTTCAAGACCCCACAGGACGTCGCCCAGCAAGTGAc        <  2:83922/143‑1 (MQ=255)
                                                                                                                                           tGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGAc                                   <  1:561375/118‑1 (MQ=255)
                                                                                                                                           tGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGAc                                   >  2:561375/1‑118 (MQ=255)
                                                                                                                                           tGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGAACCAACAGGAGCGTGGCCAGCAAGTGACCagat   >  1:83922/1‑147 (MQ=255)
                                                                                                                                            gATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGACGGCGGCCAGCAAGTGACCAATCt  <  1:476118/150‑1 (MQ=255)
                                                                                                                                                 aCCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGACGGCGGCCAGCaa             >  1:988346/1‑134 (MQ=255)
                                                                                                                                                 aCCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGACGGCGGCCAGCaa             <  2:988346/134‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CTGCGCTTGCACGACTGGCTGCGCAATCACCAACCCGACGAAGCACCGCCCGATGGCTACGGCCTGTGGGGCGAGGTATTCGATAGCATCTACCACCTGCAACGCCGCGACCAGCGCGTACGGGGCCGCCTGCAGGCAGTGATCGACCGGGTGCAGGAGTCCACCGCTGCGCTGCGTGACGCGGTGATCATGCTCGACAGCGACGGCAACCTGGAATGGTGGAACCGCGCCGCCGAGACTCTGCTGGGCTTCAAGACCCCACAGGACGGCGGCCAGCAAGTGACCAATCT  >  NC_002947/6066080‑6066369

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