Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,756 G→T 5.6% A228E (GCG→GAG) ‡ PP_3503 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,7560GT 5.6% 1028.0 / 33.7 304A228E (GCG→GAG) ‡PP_3503sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (150/137);  new base T (10/7);  total (160/144)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 308 total aligned reads displayed.

GCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAG  >  NC_002947/3975609‑3975905
                                                                                                                                                   |                                                                                                                                                     
gccgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCt                                                                                                                                                     >  1:9015668/1‑150 (MQ=255)
 ccgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCtt                                                                                                                                                    <  1:7139278/150‑1 (MQ=255)
 ccgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCtt                                                                                                                                                    >  2:5323881/1‑150 (MQ=255)
  cgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTc                                                                                                                                                   >  2:6849247/1‑150 (MQ=255)
       gCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATAc                                                                                                                                              <  2:1320281/150‑1 (MQ=255)
         aCCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGc                                                                                                                                            >  2:436505/1‑150 (MQ=255)
           ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                                                                          >  2:762890/1‑150 (MQ=255)
           ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                                                                          <  1:6498923/150‑1 (MQ=255)
                 aTCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCtt                                                                                                                                    >  2:3808221/1‑150 (MQ=255)
                   ccACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTgc                                                                                                                                  >  1:5508238/1‑150 (MQ=255)
                        ggTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTgg                                                                                                                             <  1:4539331/150‑1 (MQ=255)
                         gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                                                                            <  1:2067148/150‑1 (MQ=255)
                           aTGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGa                                                                                                                          >  2:5605146/1‑150 (MQ=255)
                             gggCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGc                                                                                                                        <  1:1339056/150‑1 (MQ=255)
                              ggctggctGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCa                                                                                                                       <  1:1079321/150‑1 (MQ=255)
                               gctggctgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCAc                                                                                                                      >  1:755177/1‑150 (MQ=255)
                                ctggctgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCAcc                                                                                                                     <  1:8476848/150‑1 (MQ=255)
                                 tggctgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCg                                                                                                                    >  1:3899875/1‑150 (MQ=255)
                                    ctgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTa                                                                                                                 <  1:3110945/150‑1 (MQ=255)
                                     tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGa                                                                                                                                                 <  2:2419477/117‑1 (MQ=255)
                                     tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaa                                                                                                                >  1:7750255/1‑150 (MQ=255)
                                     tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaa                                                                                                                <  1:2609267/150‑1 (MQ=255)
                                     tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaa                                                                                                                >  2:5379586/1‑150 (MQ=255)
                                      gCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaaa                                                                                                               <  1:5991290/150‑1 (MQ=255)
                                      gCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaaa                                                                                                               <  1:6426503/150‑1 (MQ=255)
                                       ccacccccCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCACTTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCGCCGGTAAAg                                                                                                              >  2:5063613/1‑150 (MQ=255)
                                       ccacccacCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAg                                                                                                              >  2:2536175/1‑150 (MQ=255)
                                         acccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgc                                                                                                            >  1:2641778/1‑150 (MQ=255)
                                          cccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcg                                                                                                           <  2:3027704/150‑1 (MQ=255)
                                           ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                                                          >  1:4213606/1‑150 (MQ=255)
                                           ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                                                          >  2:2323094/1‑150 (MQ=255)
                                           ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                                                          <  2:6521236/150‑1 (MQ=255)
                                            cacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGcc                                                                                                         <  2:7312728/150‑1 (MQ=255)
                                              cccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCt                                                                                                       >  2:1455400/1‑150 (MQ=255)
                                                cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                                                                     <  1:7272719/150‑1 (MQ=255)
                                                cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                                                                     <  2:3574630/150‑1 (MQ=255)
                                                  cTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtctc                                                                                                   <  1:3624446/150‑1 (MQ=255)
                                                  cTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtctc                                                                                                   <  1:8806394/150‑1 (MQ=255)
                                                   tCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCa                                                                                                  <  2:8356026/150‑1 (MQ=255)
                                                    cGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCAt                                                                                                 >  2:1122008/1‑150 (MQ=255)
                                                     gATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATg                                                                                                <  2:461688/150‑1 (MQ=255)
                                                     cctgtgCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCACTTCTCATg                                                                                                >  2:3643012/3‑150 (MQ=255)
                                                           gccctgcTGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGaa                                                                                          >  2:4871751/7‑150 (MQ=255)
                                                           cTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGaa                                                                                          >  2:8427292/1‑150 (MQ=255)
                                                            ttCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAAc                                                                                         <  1:6628714/150‑1 (MQ=255)
                                                             tctcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACa                                                                                        >  1:5065898/1‑150 (MQ=255)
                                                             tctcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACa                                                                                        <  2:3715167/150‑1 (MQ=255)
                                                                 cTGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTc                                                                                    <  2:4217277/150‑1 (MQ=255)
                                                                  tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCACGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGACCAGGTCg                                                                                   >  2:6984955/1‑150 (MQ=255)
                                                                  tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGc                                                                                                                                            >  2:8225438/1‑93 (MQ=255)
                                                                  tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                                                                   >  1:93679/1‑150 (MQ=255)
                                                                  tGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCg                                                                                   >  2:2734974/1‑150 (MQ=255)
                                                                  tGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCg                                                                                   >  2:7927699/1‑150 (MQ=255)
                                                                   gCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCgg                                                                                  >  1:5154167/1‑150 (MQ=255)
                                                                   gCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCgg                                                                                  >  1:8127434/1‑150 (MQ=255)
                                                                    caacacaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGa                                                                                 <  1:9009803/144‑1 (MQ=255)
                                                                    cAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGa                                                                                 >  1:8720388/1‑150 (MQ=255)
                                                                    cAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGa                                                                                 <  1:7998741/150‑1 (MQ=255)
                                                                      acacaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGACt                                                                               >  1:1219504/5‑150 (MQ=255)
                                                                      aGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACt                                                                               <  2:1909238/150‑1 (MQ=255)
                                                                        aaaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTgg                                                                             >  2:606387/1‑150 (MQ=255)
                                                                        aaaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTgg                                                                             <  1:8188524/150‑1 (MQ=255)
                                                                         aaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGa                                                                            <  2:7488613/150‑1 (MQ=255)
                                                                         aaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGa                                                                            >  1:578925/1‑150 (MQ=255)
                                                                         aaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGa                                                                            <  2:6713128/150‑1 (MQ=255)
                                                                          aCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGACTCGAt                                                                           <  2:454412/150‑1 (MQ=255)
                                                                             cagcaAGTTTGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                        <  2:890841/150‑1 (MQ=255)
                                                                             cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                        >  1:2155657/1‑150 (MQ=255)
                                                                             cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                        <  1:4183688/150‑1 (MQ=255)
                                                                             cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                        <  2:8604180/150‑1 (MQ=255)
                                                                              agcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAg                                                                       >  2:8561027/1‑150 (MQ=255)
                                                                               gcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGt                                                                      >  2:5699667/1‑150 (MQ=255)
                                                                                         cTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGAt                                                            <  1:1569439/150‑1 (MQ=255)
                                                                                         cTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGAt                                                            >  2:187067/1‑150 (MQ=255)
                                                                                           gCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCg                                                          >  2:8274877/1‑150 (MQ=255)
                                                                                             aaTTCCAACGGCAGGTCGCCGATTTCATCAAGAAACACTGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcc                                                        <  2:7158074/150‑1 (MQ=255)
                                                                                              aTTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGccc                                                       <  1:4934512/150‑1 (MQ=255)
                                                                                               ttCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcccc                                                      <  1:1416901/150‑1 (MQ=255)
                                                                                                  cAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCa                                                   <  2:3394926/150‑1 (MQ=255)
                                                                                                           gTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTCGATCAACTGCGGGGCGATCCCCGAGCACTTGATcca                                          <  1:6824864/150‑3 (MQ=38)
                                                                                                             cGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGa                                        <  1:5264911/150‑1 (MQ=255)
                                                                                                             cGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGa                                        <  2:8804331/150‑1 (MQ=255)
                                                                                                              gCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGAt                                       >  1:1711616/1‑150 (MQ=255)
                                                                                                                   tttCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACg                                  >  1:4771845/1‑150 (MQ=255)
                                                                                                                    ttCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACgg                                 <  1:436505/150‑1 (MQ=255)
                                                                                                                      catcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgct                               <  2:1052405/150‑1 (MQ=255)
                                                                                                                        tcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgctgc                             >  1:2905357/1‑150 (MQ=255)
                                                                                                                          aaGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTc                           <  1:8347540/150‑1 (MQ=255)
                                                                                                                          aaGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTc                           >  2:3908462/1‑150 (MQ=255)
                                                                                                                            gAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGc                         <  2:1220053/150‑1 (MQ=255)
                                                                                                                               aCAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGc                      <  1:5854331/150‑1 (MQ=255)
                                                                                                                                      gccgccATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGa               >  1:8246876/1‑150 (MQ=255)
                                                                                                                                        cgccATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACt             >  2:290290/1‑150 (MQ=255)
                                                                                                                                        cgccATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACt             >  2:2093272/1‑150 (MQ=255)
                                                                                                                                          ccATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTgg           >  2:2481350/1‑150 (MQ=255)
                                                                                                                                            aTCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCg         <  1:7046991/150‑1 (MQ=255)
                                                                                                                                             tCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCgg        <  1:2381657/150‑1 (MQ=255)
                                                                                                                                             tCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCgg        <  2:3127442/150‑1 (MQ=255)
                                                                                                                                              cGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGt       <  1:7381522/150‑1 (MQ=255)
                                                                                                                                                gCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGc     <  2:5894616/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAG  >  NC_002947/3975609‑3975905

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: