Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,416,554 A→T 5.4% intergenic (+87/+28) 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,5540AT 5.4% 462.7 / 11.4 149intergenic (+87/+28)pdxB/PP_2118erythronate‑4‑phosphate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (80/61);  new base T (3/5);  total (83/66)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 150 total aligned reads displayed.

ACGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTGC  >  NC_002947/2416409‑2416702
                                                                                                                                                  |                                                                                                                                                     
aCGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAgc                                                                                                                                                    >  8:18106/1‑149 (MQ=255)
 cGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAgcg                                                                                                                                                   >  2:68125/1‑149 (MQ=255)
         cgccgtAGTTCGCGCAGTTGGTCAGTGCCCTGGGTTGGGTGTTGGTCTAATCACGCCTTGTCTAGACCTTGCAGGTAGGTGTGGATACCTGTTGGGGCCGCCAAGCCCACTTGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGttt                                                                                                                                           <  5:241967/149‑2 (MQ=255)
         cgccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGttt                                                                                                                                           <  5:239007/149‑2 (MQ=255)
          gccgcAGTTGGCTCAGTTTGTCAGTGGCCTGGGTGGGGTGTTGGTCTAATCCACGCCTTGTCTTTTCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTg                                                                                                                                           <  1:305436/149‑1 (MQ=255)
          gccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGc                                                                                                                                          <  1:68125/149‑1 (MQ=255)
           ccgcAGTTGGCGCAGTTGTTCAGTGCCCTGGGGTGGGTGTTGGTCTAAGCACGCATTGTCTTGACCTTGCAGGTATGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCa                                                                                                                                         <  3:60490/149‑1 (MQ=255)
           ccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCa                                                                                                                                         <  4:97361/149‑1 (MQ=255)
            cgcagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTg                                                                                                                                           >  4:119849/1‑146 (MQ=255)
             gcagttgtcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCaga                                                                                                                                       <  8:273819/149‑1 (MQ=255)
             gcagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCaga                                                                                                                                       <  2:45093/149‑1 (MQ=255)
             gcagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCaga                                                                                                                                       <  3:123725/149‑1 (MQ=255)
             aactgcggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCaga                                                                                                                                       <  6:351150/143‑1 (MQ=255)
               agttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGtttttattg                                                                                                                                     <  4:156008/149‑8 (MQ=255)
                gttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGaa                                                                                                                                    <  7:146983/149‑1 (MQ=255)
                gttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGaa                                                                                                                                    <  3:200324/149‑1 (MQ=255)
                      gcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGc                                                                                                                                          >  8:156857/1‑137 (MQ=255)
                        agttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGa                                                                                                                            >  2:84836/1‑149 (MQ=255)
                          ttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGAt                                                                                                                           >  8:266572/1‑148 (MQ=255)
                             gTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCaga                                                                                                                                       <  1:242953/133‑1 (MQ=255)
                             gTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAAc                                                                                                                       >  7:236135/1‑149 (MQ=255)
                              tCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAAc                                                                                                                       >  2:329137/1‑148 (MQ=255)
                                aGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCtt                                                                                                                    >  1:355715/1‑149 (MQ=255)
                                         ggtggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCgg                                                                                                                                               >  5:353135/1‑113 (MQ=255)
                                         ggtggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCgg                                                                                                                                               <  6:353135/113‑1 (MQ=255)
                                              ggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTc                                                                                                       >  4:148850/1‑148 (MQ=255)
                                                      tCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCt                                                                                              <  4:329209/149‑1 (MQ=255)
                                                        tAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCtt                                                                                             >  5:140235/1‑148 (MQ=255)
                                                        tAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCtt                                                                                             >  1:135911/1‑148 (MQ=255)
                                                            gCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttctt                                                                                         >  4:161575/1‑149 (MQ=255)
                                                            gCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttctt                                                                                         >  2:351141/1‑149 (MQ=255)
                                                             cACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttcttc                                                                                        <  1:84836/149‑1 (MQ=255)
                                                             cACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTttcttc                                                                                        <  3:267752/149‑1 (MQ=255)
                                                              aCGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGa                                                                                                                            <  7:46798/112‑1 (MQ=255)
                                                              aCGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGa                                                                                                                            >  8:46798/1‑112 (MQ=255)
                                                                  cTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAggg                                                                                   <  1:351141/149‑1 (MQ=255)
                                                                  cTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAggg                                                                                   <  8:236135/149‑1 (MQ=255)
                                                                     gTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGcc                                                                                 <  3:158775/148‑1 (MQ=255)
                                                                      tCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCttt                                                                                                                   <  4:133926/113‑1 (MQ=255)
                                                                      tCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCttt                                                                                                                   >  3:133926/1‑113 (MQ=255)
                                                                      tCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCtt                                                                              <  5:256932/149‑1 (MQ=255)
                                                                        ttGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAAc                                                                                                                       >  4:200940/1‑107 (MQ=255)
                                                                           aCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCttt                                                                                            >  1:109797/1‑131 (MQ=255)
                                                                           aCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCttt                                                                                            <  2:109797/131‑1 (MQ=255)
                                                                              ttGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGTTCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCct                                                                      >  1:303791/1‑149 (MQ=255)
                                                                              ttGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCct                                                                      >  3:24136/1‑149 (MQ=255)
                                                                                gCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCctgc                                                                    >  2:53368/1‑149 (MQ=255)
                                                                                  aggtaggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCctgctg                                                                  >  7:53172/1‑149 (MQ=255)
                                                                                     taggtGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGat                                                               <  5:351765/149‑1 (MQ=255)
                                                                                        gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcat                                                            >  1:168220/1‑149 (MQ=255)
                                                                                        gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatca                                                             >  3:5839/1‑148 (MQ=255)
                                                                                        gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatc                                                              >  4:175537/1‑147 (MQ=255)
                                                                                         tgtgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcatc                                                           <  4:24136/149‑1 (MQ=255)
                                                                                         tgtgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcatc                                                           <  3:161575/149‑1 (MQ=255)
                                                                                           tgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTg                                                         >  6:202189/1‑149 (MQ=255)
                                                                                           tgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTg                                                         >  2:55697/1‑149 (MQ=255)
                                                                                                 ccTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGt                                                   >  1:21225/1‑149 (MQ=255)
                                                                                                  cTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGt                                                   >  7:259031/1‑148 (MQ=255)
                                                                                                   tGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGa                                                 >  2:346434/1‑149 (MQ=255)
                                                                                                    gTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcat                                                            >  1:332685/1‑137 (MQ=255)
                                                                                                        gggCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCgg                                             >  5:29591/1‑148 (MQ=255)
                                                                                                        gggCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCgg                                             >  3:311478/1‑148 (MQ=255)
                                                                                                         ggcggcCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTATGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGc                                            <  1:55299/149‑1 (MQ=255)
                                                                                                          gcggcCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCAc                                          >  7:135011/1‑149 (MQ=255)
                                                                                                              ccAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTcc                                                                                       <  7:293578/101‑1 (MQ=255)
                                                                                                              ccAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTcg                                      >  3:84231/1‑149 (MQ=255)
                                                                                                               cAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTcgc                                     >  7:331345/1‑149 (MQ=255)
                                                                                                               cAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTcgc                                     >  6:59285/1‑149 (MQ=255)
                                                                                                               cAAGCCAACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTcgc                                     >  7:108822/1‑149 (MQ=255)
                                                                                                                  gCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcc                                  >  5:126955/1‑149 (MQ=255)
                                                                                                                    ccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcccc                                >  3:326156/1‑149 (MQ=255)
                                                                                                                    ccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcccc                                >  3:362655/1‑149 (MQ=255)
                                                                                                                     cACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGAt                                                >  7:154846/1‑132 (MQ=255)
                                                                                                                      aCTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCtt                              <  8:8657/149‑1 (MQ=255)
                                                                                                                       ctggctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTg                             <  4:362655/149‑1 (MQ=255)
                                                                                                                        tggctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTgg                            >  2:235927/1‑149 (MQ=255)
                                                                                                                          gctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGc                           <  2:135911/148‑1 (MQ=255)
                                                                                                                          gctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGc                           <  4:43214/148‑1 (MQ=255)
                                                                                                                          gctggGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCg                          <  3:176861/149‑1 (MQ=255)
                                                                                                                             gggCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGccc                                 >  3:299378/1‑139 (MQ=255)
                                                                                                                             gggCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGccc                                 <  4:299378/139‑1 (MQ=255)
                                                                                                                              ggCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCga                      <  2:21225/149‑1 (MQ=255)
                                                                                                                               gCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgat                     >  5:219687/1‑149 (MQ=255)
                                                                                                                                ccAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgat                     >  5:29265/1‑148 (MQ=255)
                                                                                                                                 cattAAAAACCCGGCGCAAGCGCCGGGTTGCTGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTGTCTTCCAGGGCCTTGCAGGCCTGCTGGATCTTCTGCTGGGTGCTCGGCACTTCGCGCGCGTTTGCGTAgttga                   >  6:205313/4‑149 (MQ=255)
                                                                                                                                 cAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcccc                                <  2:127586/136‑1 (MQ=255)
                                                                                                                                  aaTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgatgat                  <  4:311478/149‑1 (MQ=255)
                                                                                                                                  aaTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgatgat                  >  3:105320/1‑149 (MQ=255)
                                                                                                                                   aTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGTATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                 <  8:108822/149‑1 (MQ=255)
                                                                                                                                   aTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCgg                                                    <  1:329582/114‑1 (MQ=255)
                                                                                                                                   aTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCgg                                                    >  2:329582/1‑114 (MQ=255)
                                                                                                                                   aTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                 <  8:114204/149‑1 (MQ=255)
                                                                                                                                   aTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                 >  1:184599/1‑149 (MQ=255)
                                                                                                                                     aaaaaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCg                >  3:82908/1‑148 (MQ=255)
                                                                                                                                      aaaaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGTGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGaa              >  6:41789/1‑149 (MQ=255)
                                                                                                                                        aaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAAcc            <  2:303791/149‑1 (MQ=255)
                                                                                                                                           ccGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTg                             >  7:160064/1‑129 (MQ=255)
                                                                                                                                           ccGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTg                             <  8:160064/129‑1 (MQ=255)
                                                                                                                                           ccGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCa         >  7:10391/1‑149 (MQ=255)
                                                                                                                                            cGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                 <  2:302857/140‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACGGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTGC  >  NC_002947/2416409‑2416702

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: