Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,573,595 C→G 6.6% C117W (TGC→TGG) ‡ PP_2257 → aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,573,5950CG 6.6% 230.2 / 4.0 72C117W (TGC→TGG) ‡PP_2257aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base G (3/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTGGGGGC  >  NC_002947/2573454‑2573741
                                                                                                                                             |                                                                                                                                                  
gcgcTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAcg                                                                                                                                              <  1:3712/148‑1 (MQ=255)
  gcTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAcgcgc                                                                                                                                           >  7:226017/1‑149 (MQ=255)
        tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGTGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGGGGGTGACGAGTCGGTGCGCGTCTGCCCTACGCGCGAGc                                                                                                                                       >  4:276230/1‑147 (MQ=255)
         gAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAg                                                                                                                                    >  6:64129/1‑149 (MQ=255)
          aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGt                                                                                                                                   <  5:64129/149‑1 (MQ=255)
          aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGt                                                                                                                                   <  8:226017/149‑1 (MQ=255)
               cGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGccc                                                                                                                              >  8:365196/1‑149 (MQ=255)
                ggCAAGAGCTGGATGGGTATCGTCAAGAACCCCTCCCGTAACGGCGACAACTACTGGGTCAATGCCAACGTCACGCCCATCCTCGAAGGAGGCCTGGTAGTGGGTTACGAGTCGCTGAGCGTAAGCCCTACCCGCGAGCAAGTAGcact                                                                                                                             >  6:307434/1‑146 (MQ=255)
                 gCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGt                                                                                                                            <  5:71524/149‑1 (MQ=255)
                 gCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGt                                                                                                                            >  2:230588/1‑149 (MQ=255)
                   aaGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTAGTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGGGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                                                          >  6:114725/1‑149 (MQ=255)
                       gCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGaa                                                                                                                      <  5:243087/149‑1 (MQ=255)
                       gCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGaa                                                                                                                      <  6:367279/149‑1 (MQ=255)
                            aTGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGTGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCt                                                                                                                 <  4:149616/149‑1 (MQ=255)
                             tGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTg                                                                                                                >  2:144147/1‑149 (MQ=255)
                              gggTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTg                                                                                                                >  3:56900/1‑148 (MQ=255)
                                    cGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATgcgc                                                                                                         >  6:333822/1‑149 (MQ=255)
                                     gTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCTCGAGCAAGTAGCCCGTGcc                                                                                                                         <  2:309834/132‑1 (MQ=255)
                                     gTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCTCGAGCAAGTAGCCCGTGcc                                                                                                                         >  1:309834/1‑132 (MQ=255)
                                     gTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATgcgc                                                                                                         <  8:168028/148‑1 (MQ=255)
                                        aaGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCt                                                                                                     >  2:129919/1‑149 (MQ=255)
                                             ccGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTgcgc                                                                                                 >  7:62760/1‑148 (MQ=255)
                                               gCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAAc                                                                                              >  4:178900/1‑149 (MQ=255)
                                                   ccGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCGAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGc                                                                                           >  4:138046/1‑148 (MQ=255)
                                                   ccGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGcc                                                                                          <  8:90369/149‑1 (MQ=255)
                                                    cGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                                                          >  5:258552/1‑116 (MQ=255)
                                                    cGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                                                          <  6:258552/116‑1 (MQ=255)
                                                    cGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGc                                                                                           >  5:5226/1‑147 (MQ=255)
                                                          ggCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAggcggc                                                                                   <  8:62760/149‑1 (MQ=255)
                                                           gCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCa                                                                                  <  3:346897/149‑1 (MQ=255)
                                                            cGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCAt                                                                                 <  7:365196/149‑1 (MQ=255)
                                                                   tactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACgg                                                                                            >  4:265665/1‑131 (MQ=255)
                                                                   tactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACgg                                                                                            <  3:265665/131‑1 (MQ=255)
                                                                     ctGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcg                                                                         <  1:230588/148‑1 (MQ=255)
                                                                      tGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcgcc                                                                       >  2:233586/1‑149 (MQ=255)
                                                                        ggTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcgccgc                                                                     >  1:177225/1‑149 (MQ=255)
                                                                                 cTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCt                                                              >  1:371911/1‑147 (MQ=255)
                                                                                 cTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGa                                                            <  1:233586/149‑1 (MQ=255)
                                                                                    cGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGc                                                          <  2:351137/148‑1 (MQ=255)
                                                                                    cGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGc                                                          <  2:360388/148‑1 (MQ=255)
                                                                                        ccGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCt                                                     >  7:30186/2‑149 (MQ=255)
                                                                                        ccGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCt                                                     <  8:30186/148‑1 (MQ=255)
                                                                                            ccATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgcg                                                 >  3:328500/1‑149 (MQ=255)
                                                                                            ccATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgcg                                                 >  4:157482/1‑149 (MQ=255)
                                                                                             cATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgcgc                                                <  5:114725/149‑1 (MQ=255)
                                                                                               tCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCGAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAg                                              <  3:138046/149‑1 (MQ=255)
                                                                                                  tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc                                           >  6:354746/1‑149 (MQ=255)
                                                                                                  tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc                                           >  7:342457/1‑149 (MQ=255)
                                                                                                  tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc                                           >  3:360982/1‑149 (MQ=255)
                                                                                                  tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggc                                            >  4:342624/1‑148 (MQ=255)
                                                                                                     atggTGGGCGGGTTGTGTGTTACGAGTCGGTGCGCGTTTGCCCTTCGCTCGATCAAGTAGCCGGTGCCGAAGCCCTGTATGCGCGTCTGAGATACGGCGAGGCGGCATTGTCGGGGCCGCGCCTGTTGAGCTTGCTGGCGCAggccttg                                        <  5:307434/147‑1 (MQ=255)
                                                                                                           ggCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGgccttggccttg                                  >  7:114940/1‑149 (MQ=255)
                                                                                                            gCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCa                                               <  8:200085/135‑1 (MQ=255)
                                                                                                            gCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCa                                               >  7:200085/1‑135 (MQ=255)
                                                                                                              gggTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCa                               <  3:276230/149‑1 (MQ=255)
                                                                                                                   gTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTgg                           <  6:5226/148‑1 (MQ=255)
                                                                                                                      ggTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCg                        <  1:129919/148‑1 (MQ=255)
                                                                                                                        ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     <  2:371911/149‑1 (MQ=255)
                                                                                                                        ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     >  6:43590/1‑149 (MQ=255)
                                                                                                                        ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     >  2:137916/1‑149 (MQ=255)
                                                                                                                        ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     <  3:157482/149‑1 (MQ=255)
                                                                                                                        ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     >  5:57168/1‑149 (MQ=255)
                                                                                                                              gTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGgctgc               >  3:74016/1‑149 (MQ=255)
                                                                                                                                  tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     >  5:114678/3‑139 (MQ=255)
                                                                                                                                  tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     <  6:114678/137‑1 (MQ=255)
                                                                                                                                  tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     <  6:109326/137‑1 (MQ=255)
                                                                                                                                  tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                     >  5:109326/3‑139 (MQ=255)
                                                                                                                                  tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCATGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGc           >  2:101697/3‑149 (MQ=255)
                                                                                                                                       cgTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTgg      <  1:144147/149‑1 (MQ=255)
                                                                                                                                        gTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTggg     <  6:57168/149‑1 (MQ=255)
                                                                                                                                          aTGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTggggg   >  3:339639/1‑149 (MQ=255)
                                                                                                                                            gCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTGGGGGc  >  7:184928/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTGGGGGC  >  NC_002947/2573454‑2573741

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: