Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,416,555 A→T 5.9% intergenic (+88/+27) pdxB → / ← PP_2118 erythronate‑4‑phosphate dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,416,5550AT 5.9% 257.3 / 4.9 85intergenic (+88/+27)pdxB/PP_2118erythronate‑4‑phosphate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/41);  new base T (3/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

ACGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTG  >  NC_002947/2416409‑2416701
                                                                                                                                                  |                                                                                                                                                  
aCGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAgc                                                                                                                                                  >  3:243878/1‑149 (MQ=255)
     gAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGCCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCgg                                                                                                                                             <  1:318796/149‑1 (MQ=255)
      aGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCggg                                                                                                                                            <  2:164120/149‑1 (MQ=255)
      aGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCggg                                                                                                                                            <  5:319891/149‑1 (MQ=255)
        gcgcCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGTTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGtt                                                                                                                                          <  8:227447/149‑1 (MQ=255)
        gcgcCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGtt                                                                                                                                          <  4:23725/149‑1 (MQ=255)
         cgccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGttt                                                                                                                                         <  4:227532/149‑2 (MQ=255)
          gccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTg                                                                                                                                         >  2:42243/1‑148 (MQ=255)
          gccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTg                                                                                                                                         >  2:47759/1‑148 (MQ=255)
            cgcagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCag                                                                                                                                      <  4:153242/149‑1 (MQ=255)
             gcagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCGGGGTTGCaga                                                                                                                                     <  6:68668/149‑1 (MQ=255)
             gcagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCaga                                                                                                                                     <  2:88088/149‑1 (MQ=255)
                 ttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTTGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAAACCGGAGCAAGCGCCGGGTTGAAGAGAAt                                                                                                                                 >  2:229977/1‑149 (MQ=255)
                     cgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGc                                                                                                                             <  4:68244/149‑1 (MQ=255)
                     cgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGc                                                                                                                             <  8:76658/149‑1 (MQ=255)
                        agttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGa                                                                                                                          <  1:229977/149‑1 (MQ=255)
                             gTCAGTGCCCTGTGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGACCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAAc                                                                                                                     <  3:172585/149‑1 (MQ=255)
                             gTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCaa                                                                                                                      <  8:282171/148‑1 (MQ=255)
                                aGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCtt                                                                                                                  >  1:189926/1‑149 (MQ=255)
                                  tGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTg                                                                                                                >  1:299429/1‑149 (MQ=255)
                                       tgggtgggGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTAc                                                                                                           >  4:197868/1‑149 (MQ=255)
                                                  ttGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGacac                                                                                                <  2:299429/149‑1 (MQ=255)
                                                     gTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGc                                                                                             >  8:69589/1‑149 (MQ=255)
                                                     gTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGc                                                                                             >  6:275586/1‑149 (MQ=255)
                                                     gTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGc                                                                                             >  4:288870/1‑149 (MQ=255)
                                                      tCTAAGCACGGCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCGCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCt                                                                                            <  1:155063/149‑1 (MQ=255)
                                                      tCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCt                                                                                            >  1:99491/1‑149 (MQ=255)
                                                      tCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCt                                                                                            <  3:197868/149‑1 (MQ=255)
                                                          aGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttct                                                                                        <  1:47759/149‑1 (MQ=255)
                                                            cACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttcttc                                                                                      <  3:150143/149‑1 (MQ=255)
                                                            cACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttcttc                                                                                      >  2:163993/1‑149 (MQ=255)
                                                             aCGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttcttc                                                                                      >  1:48835/1‑148 (MQ=255)
                                                             aCGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTcc                                                                                     <  5:171112/149‑1 (MQ=255)
                                                                 cTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAggg                                                                                 <  3:290345/149‑1 (MQ=255)
                                                                      cgatctcTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTg                                                                            <  3:176125/143‑1 (MQ=255)
                                                                      cTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGCTCAACCTTTGCGTACGCGCTCGACACGGCGTGCTTCCAGGGCCTTg                                                                            >  4:195636/1‑149 (MQ=255)
                                                                         gACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAg                                                                         >  1:49876/1‑149 (MQ=255)
                                                                          aCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAgg                                                                        >  8:3395/1‑149 (MQ=255)
                                                                          aCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAg                                                                         <  3:195636/148‑1 (MQ=255)
                                                                             ttGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTg                                                                            >  4:176125/1‑142 (MQ=255)
                                                                               gCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGa                                                                                                                          >  3:212993/1‑94 (MQ=255)
                                                                               gCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGa                                                                                                                          <  4:212993/94‑1 (MQ=255)
                                                                                 aggtaggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCctgctg                                                                 <  7:69589/149‑1 (MQ=255)
                                                                                   gtaggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGa                                                               >  8:217646/1‑149 (MQ=255)
                                                                                   gtaggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGa                                                               <  4:226868/149‑1 (MQ=255)
                                                                                       gtgtGGTTACCTGTTTGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCttt                                                                                          >  5:286118/1‑118 (MQ=255)
                                                                                       gtgtGGTTACCTGTTTGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCttt                                                                                          <  6:286118/118‑1 (MQ=255)
                                                                                       gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcat                                                           >  5:67175/1‑149 (MQ=255)
                                                                                       gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcat                                                           <  2:156779/149‑1 (MQ=255)
                                                                                        tgtgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcatc                                                          <  2:99491/149‑1 (MQ=255)
                                                                                               aCCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCg                                                    <  8:249578/148‑1 (MQ=255)
                                                                                                 cTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGt                                                  >  4:267178/1‑148 (MQ=255)
                                                                                                   gTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGAt                                               >  7:115232/1‑149 (MQ=255)
                                                                                                    ttGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTc                                                     <  6:76290/142‑1 (MQ=255)
                                                                                                    ttGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTc                                                     >  5:76290/1‑142 (MQ=255)
                                                                                                    ttGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATc                                              <  1:42243/149‑1 (MQ=255)
                                                                                                    ttGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATc                                              <  2:189926/149‑1 (MQ=255)
                                                                                                       gggCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCgg                                            <  2:48835/148‑1 (MQ=255)
                                                                                                       gggCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCgg                                            >  7:129127/1‑148 (MQ=255)
                                                                                                             ccAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTcg                                     >  1:161642/1‑149 (MQ=255)
                                                                                                                 gCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcc                                 >  2:308127/1‑149 (MQ=255)
                                                                                                                  cccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGccc                                >  8:116677/1‑149 (MQ=255)
                                                                                                                   ccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcccc                               >  6:67170/1‑149 (MQ=255)
                                                                                                                    cACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCt                              >  7:206830/1‑149 (MQ=255)
                                                                                                                     aCTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCtt                             <  3:267178/149‑1 (MQ=255)
                                                                                                                      ctggctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTg                            <  2:252550/149‑1 (MQ=255)
                                                                                                                      ctggctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTg                            <  7:3395/149‑1 (MQ=255)
                                                                                                                       tggctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTgg                           >  4:211869/1‑149 (MQ=255)
                                                                                                                         gctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTgg                           <  7:116677/147‑1 (MQ=255)
                                                                                                                             ggCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCCTCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCga                     >  4:219865/1‑149 (MQ=255)
                                                                                                                               ccAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgat                    >  6:53020/1‑148 (MQ=255)
                                                                                                                                cAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCTTCgatga                  >  2:269204/1‑149 (MQ=255)
                                                                                                                                cAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgatga                  >  1:116873/1‑149 (MQ=255)
                                                                                                                                 aaTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgatgat                 <  5:153362/149‑1 (MQ=255)
                                                                                                                                  aTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                <  2:161642/149‑1 (MQ=255)
                                                                                                                                     aaaaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGaa             >  6:111905/1‑149 (MQ=255)
                                                                                                                                       aaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAAcc           >  4:317026/1‑149 (MQ=255)
                                                                                                                                       aaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAAcc           >  6:243936/1‑149 (MQ=255)
                                                                                                                                         cccGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCg          >  5:199495/1‑148 (MQ=255)
                                                                                                                                         cccGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAg                                                                         >  4:95276/1‑85 (MQ=255)
                                                                                                                                         cccGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAg                                                                         <  3:95276/85‑1 (MQ=255)
                                                                                                                                         cccGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGc         <  7:217646/149‑1 (MQ=255)
                                                                                                                                          ccGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCa        >  3:19346/1‑149 (MQ=255)
                                                                                                                                           cGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCAt       <  4:243878/149‑1 (MQ=255)
                                                                                                                                                cAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTg  >  4:115967/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTG  >  NC_002947/2416409‑2416701

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