Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,284 C→G 6.5% E121Q (GAG→CAG)  PP_3414 ← methyl‑accepting chemotaxis transducer/sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,2840CG 6.5% 226.7 / 3.6 72E121Q (GAG→CAG) PP_3414methyl‑accepting chemotaxis transducer/sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base G (2/3);  total (32/40)
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

CTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCT  >  NC_002947/3865138‑3865427
                                                                                                                                                  |                                                                                                                                               
cTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGt                                                                                                                                               >  5:272082/1‑149 (MQ=255)
cTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGt                                                                                                                                               <  8:157850/149‑1 (MQ=255)
    aTCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAg                                                                                                                                           >  5:341457/1‑149 (MQ=255)
        cccGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACg                                                                                                                                       >  4:146741/1‑149 (MQ=255)
        cccGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACg                                                                                                                                       <  5:70043/149‑1 (MQ=255)
        cccGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACg                                                                                                                                       >  3:144995/1‑149 (MQ=255)
           ggCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCGCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACCGGc                                                                                                                                    <  2:157020/149‑1 (MQ=255)
                 ggggggTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGt                                                                                                                                <  5:10459/147‑1 (MQ=255)
                  gggggTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAg                                                                                                                              <  7:186903/148‑1 (MQ=255)
                   ggggTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCt                                                                                                                            <  6:253343/149‑1 (MQ=255)
                    gggTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTg                                                                                                                           <  4:123216/149‑1 (MQ=255)
                          aCTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGACACCTTGATTACGCTGGTGACTTTCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTc                                                                                                                      >  6:322548/1‑148 (MQ=255)
                            tCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTccag                                                                                                                   >  2:175836/1‑149 (MQ=255)
                            tCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTccag                                                                                                                   >  7:17178/1‑149 (MQ=255)
                                              cGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCtt                                                                                                                       <  6:233702/127‑1 (MQ=255)
                                              cGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCtt                                                                                                                       >  5:233702/1‑127 (MQ=255)
                                               ggCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGcc                                                                                                <  8:89903/149‑1 (MQ=255)
                                                  tGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTg                                                                                             >  4:69347/1‑149 (MQ=255)
                                                    aTGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGt                                                                                            <  8:17178/148‑1 (MQ=255)
                                                    aTGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGt                                                                                            <  6:111001/148‑1 (MQ=255)
                                                    aTGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGt                                                                                            <  8:243904/148‑1 (MQ=255)
                                                    aTGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTc                                                                                           >  6:239966/1‑149 (MQ=255)
                                                       gCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGa                                                                                         <  5:360137/148‑1 (MQ=255)
                                                       gCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGAt                                                                                        <  3:280923/149‑1 (MQ=255)
                                                        ccTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATa                                                                                       <  4:144995/149‑1 (MQ=255)
                                                              gCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGGGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTc                                                                                  >  3:321000/1‑148 (MQ=255)
                                                               cAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGt                                                                                            >  1:192179/1‑137 (MQ=255)
                                                               cAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGt                                                                                            <  2:192179/137‑1 (MQ=255)
                                                                 ggCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTAGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAAc                                                                              <  6:341457/149‑1 (MQ=255)
                                                                  gCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTc                                                                                  >  6:250535/1‑144 (MQ=255)
                                                                  gCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTc                                                                                  <  5:250535/144‑1 (MQ=255)
                                                                   cTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGt                                                                            <  3:156684/149‑1 (MQ=255)
                                                                       cGCTTTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCg                                                                        >  6:374865/1‑149 (MQ=255)
                                                                          tctcGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCt                                                                      >  2:333554/1‑148 (MQ=255)
                                                                              gTGCTCCTGCACGACCCGCTGGGTGATGTCCGACGCCACCTTGATTCCGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGa                                                                    <  4:243141/146‑1 (MQ=255)
                                                                              gTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAg                                                                 >  5:355457/1‑149 (MQ=255)
                                                                                  tCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTcca                                                                                                                    <  7:174166/94‑1 (MQ=255)
                                                                                  tCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTcca                                                                                                                    >  8:174166/1‑94 (MQ=255)
                                                                                  tCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCtt                                                             >  8:296044/1‑149 (MQ=255)
                                                                                       caccacCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccc                                                         <  1:175836/148‑1 (MQ=255)
                                                                                       caccacCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgc                                                           <  8:286491/146‑1 (MQ=255)
                                                                                       caccacCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgc                                                           >  7:286491/1‑146 (MQ=255)
                                                                                       caccacCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCAc                                                                                                                                        <  5:330689/69‑1 (MQ=38)
                                                                                       caccacCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCAc                                                                                                                                        >  6:330689/1‑69 (MQ=255)
                                                                                           accCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccctgccc                                                    >  5:284881/1‑149 (MQ=255)
                                                                                           accCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccctgccc                                                    <  3:146741/149‑1 (MQ=255)
                                                                                               gCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAAc                                                >  2:217009/1‑149 (MQ=255)
                                                                                                cTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACt                                               >  7:117242/1‑149 (MQ=255)
                                                                                                  gggTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTcg                                             >  2:2257/1‑149 (MQ=255)
                                                                                                     tGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCgcgc                                          <  1:217009/149‑1 (MQ=255)
                                                                                                          tCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAg                                     >  1:5582/1‑149 (MQ=255)
                                                                                                           ccGACGCCACCTTGAGTAGGCTGGTGACTTGCCGGTGGTGGTGCAGCCCGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCGCCGCCGTGGCGATACGCTCGAACGTACCGCTGATAGGCTGGCGCGGCCCCAACTCGCGccggcgg                                    >  4:370662/1‑144 (MQ=255)
                                                                                                            cGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGtt                                   <  8:117242/149‑1 (MQ=255)
                                                                                                            cGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGt                                    <  4:192784/148‑1 (MQ=255)
                                                                                                                gCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCATCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTg                               <  3:153123/149‑1 (MQ=255)
                                                                                                                gCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTg                               <  6:355457/149‑1 (MQ=255)
                                                                                                                gCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTg                               <  5:374865/149‑1 (MQ=255)
                                                                                                                 ccACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTgg                              <  6:272082/149‑1 (MQ=255)
                                                                                                                    ccTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTAc                           >  5:146734/1‑149 (MQ=255)
                                                                                                                        gATTACGCGGGTGACTTGCCGGTGCTGGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTCTTAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCg                        <  1:267804/148‑1 (MQ=255)
                                                                                                                         aTTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCgg                       <  2:229133/148‑1 (MQ=255)
                                                                                                                          ttACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGgcg                     >  1:138172/1‑149 (MQ=255)
                                                                                                                            aCGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGgcgc                    <  3:370662/148‑1 (MQ=255)
                                                                                                                                tGGTGACTTGCCTGTGCGCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCttttt               <  8:120583/149‑1 (MQ=255)
                                                                                                                                tGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCttttt               >  3:292851/1‑149 (MQ=255)
                                                                                                                                  gTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTg              <  7:296044/148‑1 (MQ=255)
                                                                                                                                  gTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGc             <  8:27325/149‑1 (MQ=255)
                                                                                                                                       ttGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGt        >  8:17759/1‑149 (MQ=255)
                                                                                                                                         gCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTcc      >  7:95965/1‑149 (MQ=255)
                                                                                                                                            ggTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGc   <  5:239966/149‑1 (MQ=255)
                                                                                                                                            ggTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGAGCTTTTTGCATAGTCCGGc   >  2:223346/1‑149 (MQ=255)
                                                                                                                                             atGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCt  <  8:18317/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCT  >  NC_002947/3865138‑3865427

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: