Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,628 G→C 7.5% intergenic (‑39/‑379) 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,6280GC 7.5% 212.2 / 6.4 67intergenic (‑39/‑379)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (3/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

AACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCC  >  NC_002947/4936480‑4936763
                                                                                                                                                      |                                                                                                                                       
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                                                                                                           <  4:177298/149‑1 (MQ=255)
                              aaCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTc                                                                                                           <  7:95235/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                                                                                                                             >  5:234792/1‑108 (MQ=255)
                                                    cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGtc                                                                                                                             <  6:234792/109‑2 (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             >  5:281704/1‑147 (MQ=255)
                                                                                                                              gCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACga             <  6:281704/147‑1 (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)
                                                                                                                                                      |                                                                                                                                       
AACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCC  >  NC_002947/4936480‑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: