Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,626 G→C 8.6% intergenic (‑37/‑381) flhA ← / → PP_4345 flagellar export pore protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,936,6260GC 8.6% 203.2 / 7.8 67intergenic (‑37/‑381)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base C (4/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCC  >  NC_002947/4936478‑4936763
                                                                                                                                                      |                                                                                                                                         
gCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAggg                                                                                                                                           >  7:84982/1‑149 (MQ=34)
  aaCATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAggggg                                                                                                                                         >  8:82965/1‑149 (MQ=34)
            ttGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCg                                                                                                                               >  4:167099/1‑149 (MQ=255)
              gCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTg                                                                                                                             >  3:177298/1‑149 (MQ=39)
              gCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTGGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGCCGCCCAGGGGGCGCCTGGGCGt                                                                                                                              <  5:94553/148‑1 (MQ=255)
               ccAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTgg                                                                                                                            >  7:159883/1‑149 (MQ=40)
                  aCATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGtcaaag                                                                                                                         <  3:167099/149‑6 (MQ=255)
                   cATCACCAGCAGCAACAGCTGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCATAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCttt                                                                                                                        <  2:194585/149‑1 (MQ=255)
                   cATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCttt                                                                                                                        <  3:214950/149‑1 (MQ=255)
                     tCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGc                                                                                                                      <  7:129193/149‑1 (MQ=255)
                          agcagcAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGAGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaag                                                                                                                   >  1:315802/1‑147 (MQ=255)
                               caACAGCGGCACCCCCAGGTTGCCCCTGCCGAGCCCGGCCAGCTTGTTGAGGGCGATGCTGATTAACTGAGTGCGAACCACCGCGAGACCTCTTGATGCAAAACTTTTACGCCCAGGCGGCGCCAGGGCCTGGCTTTGCaagatcgca                                                                                                             >  8:133806/1‑143 (MQ=14)
                               caACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGc                                                                                                                      <  7:133806/139‑1 (MQ=255)
                                aaCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTc                                                                                                           <  7:95235/149‑1 (MQ=255)
                                aaCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTc                                                                                                           <  4:177298/149‑1 (MQ=255)
                                  cAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTcc                                                                                                          >  2:103524/1‑148 (MQ=255)
                                     cGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACCGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAAc                                                                                                       <  1:225244/149‑1 (MQ=255)
                                         aCCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                  <  2:274766/149‑1 (MQ=255)
                                         aCCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                  >  2:55487/1‑149 (MQ=255)
                                          cccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGc                                                                                                 <  7:86390/149‑1 (MQ=255)
                                           ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                                                                                                <  8:84982/149‑1 (MQ=255)
                                           ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                                                                                                <  5:22073/149‑1 (MQ=255)
                                             ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                  <  2:57923/146‑1 (MQ=255)
                                             ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                  >  1:57923/1‑146 (MQ=255)
                                             ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGa                                                                                              <  6:153230/149‑1 (MQ=255)
                                                      cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGTTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCacac                                                                                     >  6:151580/1‑149 (MQ=255)
                                                      cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCacac                                                                                     >  1:311639/1‑149 (MQ=255)
                                                      cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGtc                                                                                                                             <  6:234792/109‑2 (MQ=255)
                                                      cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGtc                                                                                                                             >  5:234792/1‑108 (MQ=255)
                                                       ccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCacac                                                                                     >  3:91392/1‑148 (MQ=255)
                                                          gCCGGGCTCGGCCGGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaa                                                                                 <  6:167000/149‑1 (MQ=255)
                                                             gAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGa                                                                              <  6:178240/149‑1 (MQ=255)
                                                                   ggCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAAcc                                                                        >  1:47036/1‑149 (MQ=255)
                                                                            gttgCGGGCGTTGTTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGc                                                                <  5:151580/148‑1 (MQ=255)
                                                                                         cTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAAGTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTTCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGAc                                                  <  8:29792/149‑1 (MQ=255)
                                                                                         cTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGAc                                                  <  8:159883/149‑1 (MQ=255)
                                                                                              tAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                              <  7:82965/148‑1 (MQ=255)
                                                                                              tAACTGAGAGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                              <  2:315802/148‑1 (MQ=255)
                                                                                                aCTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                              >  7:86087/1‑146 (MQ=255)
                                                                                                aCTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                              <  8:86087/146‑1 (MQ=255)
                                                                                                aCTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTgg                                            <  2:47036/148‑1 (MQ=255)
                                                                                                 cTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCt                                          >  5:198696/1‑149 (MQ=255)
                                                                                                  tGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCt                                          >  6:279018/1‑148 (MQ=255)
                                                                                                    aGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCAt                                       >  1:190847/1‑149 (MQ=255)
                                                                                                        ctacccACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGagc                                   <  1:1114/145‑1 (MQ=255)
                                                                                                             cACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGagcaagc                               >  3:233925/1‑148 (MQ=255)
                                                                                                              aCCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                              >  6:318003/1‑132 (MQ=255)
                                                                                                              aCCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                              <  5:318003/132‑1 (MQ=255)
                                                                                                               ccGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAgcaagcaa                             >  8:52159/1‑148 (MQ=255)
                                                                                                                 gcgAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCa                          >  7:278935/1‑149 (MQ=255)
                                                                                                                      tCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCAcc                        >  6:64934/1‑146 (MQ=255)
                                                                                                                        ctctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGa                    >  3:118520/1‑148 (MQ=255)
                                                                                                                         tctTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGaa                   >  4:128211/1‑148 (MQ=255)
                                                                                                                           ttGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGaaa                  <  7:48810/147‑1 (MQ=255)
                                                                                                                             gATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGaaaaa                <  1:103524/147‑1 (MQ=255)
                                                                                                                               tGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga             <  5:64934/148‑1 (MQ=255)
                                                                                                                                gCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga             <  6:281704/147‑1 (MQ=255)
                                                                                                                                gCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga             >  5:281704/1‑147 (MQ=255)
                                                                                                                                 cAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGAc          >  7:191837/1‑149 (MQ=255)
                                                                                                                                 cAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCAACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGAc          >  2:25648/1‑149 (MQ=255)
                                                                                                                                  aaaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAATAACGAGAcc         >  5:248157/1‑149 (MQ=255)
                                                                                                                                  aaaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGAcc         <  2:190847/149‑1 (MQ=255)
                                                                                                                                   aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGa                    <  4:289249/137‑1 (MQ=255)
                                                                                                                                   aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGa                    >  3:289249/1‑137 (MQ=255)
                                                                                                                                   aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCt        >  6:191381/1‑149 (MQ=255)
                                                                                                                                      cTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGTGAGCACCACCGCGTTGa                                                   <  4:58411/103‑1 (MQ=255)
                                                                                                                                      cTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGTGAGCACCACCGCGTTGa                                                   >  3:58411/1‑103 (MQ=255)
                                                                                                                                       tttGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCg    >  7:161749/1‑149 (MQ=255)
                                                                                                                                         ccaCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGcc  >  7:68941/3‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                         
GCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCC  >  NC_002947/4936478‑4936763

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