Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,569,046 Δ1 bp 5.2% coding (847/879 nt) PP_2250 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,569,0460G. 5.2% 450.3 / 18.6 134coding (847/879 nt)PP_2250LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (65/62);  new base . (2/5);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 134 total aligned reads displayed.

GAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCC  >  NC_002947/2568898‑2569193
                                                                                                                                                     |                                                                                                                                                    
 agggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTg                                                                                                                                                    >  1:757785/1‑150 (MQ=255)
  gggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGc                                                                                                                                                   <  1:2814690/150‑1 (MQ=255)
  gggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGc                                                                                                                                                   <  1:2336317/150‑1 (MQ=255)
   cccgcACCGCGTCATGAGCGCGCTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCa                                                                                                                                                  >  1:3188330/3‑150 (MQ=255)
       cAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCct                                                                                                                                             <  1:4276987/150‑1 (MQ=255)
           gcgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGGCAGCTGCTCGATCAGCAGGCGTGCAGCctgctg                                                                                                                                         <  2:2541600/150‑1 (MQ=255)
           gcgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGGCAGCTGCTCGATCAGCAGGCGTGCAGCctgctg                                                                                                                                         <  2:2541700/150‑1 (MQ=255)
             gcTCATGACGCGGTGCGGGGCCTATCGACCGAGGTATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctgat                                                                                                                                       <  1:1571134/150‑3 (MQ=255)
             gcTCATGACGCGGTGCGGGGCCTATCGACCGAGGTATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctgat                                                                                                                                       <  1:1570867/150‑3 (MQ=255)
             gcTCATGACGCGGTGCGGGGCCTATCGACCGAGGTATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctgat                                                                                                                                       <  1:1570845/150‑3 (MQ=255)
             gcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctgat                                                                                                                                       <  1:2647407/150‑3 (MQ=255)
              cTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAg                                                                                                                                      >  1:2022111/1‑150 (MQ=255)
               tCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGc                                                                                                                                     <  1:3018093/150‑1 (MQ=255)
               tCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGc                                                                                                                                     <  1:2374626/150‑1 (MQ=255)
                 aTGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCgg                                                                                                                                   >  2:127862/1‑150 (MQ=255)
                  tGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGc                                                                                                                                  <  2:299626/150‑1 (MQ=255)
                  tGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGc                                                                                                                                  <  2:299622/150‑1 (MQ=255)
                  tGACGCGGTGCGGGGCCTATCGACCGAGAGATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGc                                                                                                                                  <  2:4758099/150‑1 (MQ=255)
                   gACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCt                                                                                                                                 <  1:1350037/150‑1 (MQ=255)
                       gcgGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCtgtt                                                                                                                              >  2:255407/1‑150 (MQ=255)
                        gggTGCGGGGCCTATCGTCCGAGATATCGGCTTGATGGTTAAAAGATTGAAGTGCTTAGGGTAAATTCTTACAGACTTTTCGAAGGACGTCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCTGCtgttg                                                                                                                             <  2:3017094/149‑1 (MQ=255)
                         ggTGCGTGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCTGCCTGCTGCAGCGGCtgttgt                                                                                                                            <  2:2843455/150‑1 (MQ=255)
                         ggTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCtgttgt                                                                                                                            <  2:5018192/150‑1 (MQ=255)
                             cGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTgc                                                                                                                        <  1:3114680/150‑1 (MQ=255)
                                 gCCTATCGACAGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCca                                                                                                                    <  2:5085381/150‑1 (MQ=255)
                                  ccTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcac                                                                                                                   >  2:814136/1‑150 (MQ=255)
                                  ccTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcac                                                                                                                   >  2:1022295/1‑150 (MQ=255)
                                   cTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                  <  1:201631/150‑1 (MQ=255)
                                   cTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                  <  1:2859652/150‑1 (MQ=255)
                                   cTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                  <  1:2034584/150‑1 (MQ=255)
                                   cTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                  <  1:3262211/150‑1 (MQ=255)
                                   cTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                  <  1:367759/150‑1 (MQ=255)
                                   cTATCGACCGAGATATCGGCGTGAGGGGTCAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                  >  2:3142006/1‑150 (MQ=255)
                                    tATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAg                                                                                                                 >  2:5077742/1‑150 (MQ=255)
                                      tCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGa                                                                                                               <  2:1621158/150‑1 (MQ=255)
                                      tCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGa                                                                                                               >  2:3137022/1‑150 (MQ=255)
                                      tCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGa                                                                                                               <  1:4073782/150‑1 (MQ=255)
                                         aCCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGatat                                                                                                            >  2:578737/1‑150 (MQ=255)
                                         aCCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCAGCTGCAGCGGCTGTTGTTTGCGCCACAGGatat                                                                                                            >  2:435372/1‑150 (MQ=255)
                                                tatCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACggg                                                                                                     >  1:666766/1‑150 (MQ=255)
                                                  tCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCa                                                                                                   <  2:4335475/150‑1 (MQ=255)
                                                   cGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAg                                                                                                  <  1:69545/150‑1 (MQ=255)
                                                    ggCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGt                                                                                                 >  1:402456/1‑150 (MQ=255)
                                                     gCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGtt                                                                                                >  1:2562669/1‑150 (MQ=255)
                                                         gAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCaa                                                                                            >  2:1274648/1‑150 (MQ=255)
                                                           ggggTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAAtt                                                                                          <  1:3978216/150‑1 (MQ=255)
                                                            gggTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTc                                                                                         >  2:1705110/1‑150 (MQ=255)
                                                             ggTAACAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCg                                                                                        >  2:2773937/1‑150 (MQ=255)
                                                             ggTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCg                                                                                        >  2:285837/1‑150 (MQ=255)
                                                               tAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      >  2:2682687/1‑150 (MQ=255)
                                                               tAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      >  1:1215690/1‑150 (MQ=255)
                                                               tAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      >  2:4639113/1‑150 (MQ=255)
                                                               tAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      <  2:5006505/150‑1 (MQ=255)
                                                                   aGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCgt                                                                                  <  1:1022295/150‑1 (MQ=255)
                                                                   aGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCgt                                                                                  <  2:4623653/150‑1 (MQ=255)
                                                                     aTTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGt                                                                                                                                                     >  2:3318037/1‑82 (MQ=255)
                                                                      ttGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTa                                                                               <  1:255407/150‑1 (MQ=255)
                                                                        gAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTAtt                                                                             <  2:1636736/150‑1 (MQ=255)
                                                                         aaGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTa                                                                            >  2:185036/1‑150 (MQ=255)
                                                                         aaGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTa                                                                            >  1:1703370/1‑150 (MQ=255)
                                                                            tGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACgg                                                                         <  1:1873685/150‑1 (MQ=255)
                                                                                 aGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTc                                                                    <  1:3128834/150‑1 (MQ=255)
                                                                                  gggAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTACAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTcc                                                                   >  1:1865768/1‑150 (MQ=255)
                                                                                   ggAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCa                                                                  >  1:5041091/1‑150 (MQ=255)
                                                                                     aaaaTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAgg                                                                >  1:3552505/1‑150 (MQ=255)
                                                                                      aaaTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGgc                                                               <  2:1022176/150‑1 (MQ=255)
                                                                                       aaTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGgcg                                                              >  1:1116839/1‑150 (MQ=255)
                                                                                        aTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCtgttgt                                                                                                                            >  2:604810/1‑87 (MQ=255)
                                                                                             tACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGGTGTTTGCGCCACAGGATCTGCACGGGCAGTTGCCATTCGTTGCGTGTATTACGGAATTCCAGGCGCAcccg                                                        >  2:3710678/1‑148 (MQ=255)
                                                                                                aGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTCCGGAATTCCAGGCGCACCAGACg                                                     >  2:4374042/1‑150 (MQ=255)
                                                                                                aGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCACACCAGACg                                                     >  1:4748228/1‑150 (MQ=255)
                                                                                                aGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCACACCAGACg                                                     >  1:4748174/1‑150 (MQ=255)
                                                                                                    gTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCg                                                 >  2:3673694/1‑150 (MQ=255)
                                                                                                             gacggCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCcagc                                        <  1:957541/150‑1 (MQ=255)
                                                                                                                   aGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCggg                                  <  1:2284624/150‑1 (MQ=255)
                                                                                                                    gTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCgggg                                 <  1:1034672/150‑1 (MQ=255)
                                                                                                                     tCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGGTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGc                                >  2:2579828/1‑150 (MQ=255)
                                                                                                                       aCGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTCTTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGa                              >  1:2917474/1‑150 (MQ=255)
                                                                                                                       aCGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGa                              >  2:2133892/1‑150 (MQ=255)
                                                                                                                       aCGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGa                              >  2:1895089/1‑150 (MQ=255)
                                                                                                                       aCGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGa                              >  2:3510322/1‑150 (MQ=255)
                                                                                                                        cGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGAc                             <  1:3137022/150‑1 (MQ=255)
                                                                                                                              agccagcTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGa                       <  1:3388579/150‑1 (MQ=255)
                                                                                                                               gccagcTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGAc                      >  2:1460094/1‑150 (MQ=255)
                                                                                                                                   gctgctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGACTTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACg                  >  2:3992515/1‑150 (MQ=255)
                                                                                                                                   gctgctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACg                  >  2:942127/1‑150 (MQ=255)
                                                                                                                                   gctgctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACg                  <  1:3203137/150‑1 (MQ=255)
                                                                                                                                     tgctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGc                >  1:5001090/1‑150 (MQ=255)
                                                                                                                                     tgctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGc                >  2:1080638/1‑150 (MQ=255)
                                                                                                                                      gctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa               <  1:578737/150‑1 (MQ=255)
                                                                                                                                      gctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa               <  1:4288574/150‑1 (MQ=255)
                                                                                                                                      gctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa               <  1:4288564/150‑1 (MQ=255)
                                                                                                                                       ctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCGCCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAg              <  2:3690521/150‑1 (MQ=255)
                                                                                                                                       ctCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAg              <  1:2524946/150‑1 (MQ=255)
                                                                                                                                         cGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAt            >  2:2838681/1‑150 (MQ=255)
                                                                                                                                           tgcagcagGCTGCACGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAtct          >  1:1411830/3‑150 (MQ=255)
                                                                                                                                            tCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAtctc         >  1:3059661/1‑150 (MQ=255)
                                                                                                                                                cagGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCa     >  1:4475382/1‑150 (MQ=255)
                                                                                                                                                cagGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCa     >  2:4065919/1‑150 (MQ=255)
                                                                                                                                                   gCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGcc  >  1:976856/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCC  >  NC_002947/2568898‑2569193

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