Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,630 G→C 7.5% intergenic (‑41/‑377) 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,6300GC 7.5% 208.8 / 6.4 67intergenic (‑41/‑377)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

TTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTACGG  >  NC_002947/4936490‑4936778
                                                                                                                                              |                                                                                                                                                    
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:55487/1‑149 (MQ=255)
                             aCCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                                                                                                                 <  2:274766/149‑1 (MQ=255)
                              cccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGc                                                                                                                <  7:86390/149‑1 (MQ=255)
                               ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                                                                                                               <  5:22073/149‑1 (MQ=255)
                               ccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                                                                                                               <  8:84982/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                                                             <  5:318003/132‑1 (MQ=255)
                                                                                                  aCCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGt                                                             >  6:318003/1‑132 (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                                   >  3:289249/1‑137 (MQ=255)
                                                                                                                       aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGa                                   <  4:289249/137‑1 (MQ=255)
                                                                                                                       aaaCTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCt                       >  6:191381/1‑149 (MQ=255)
                                                                                                                          cTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGTGAGCACCACCGCGTTGa                                                                  >  3:58411/1‑103 (MQ=255)
                                                                                                                          cTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGTGAGCACCACCGCGTTGa                                                                  <  4:58411/103‑1 (MQ=255)
                                                                                                                           tttGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCg                   >  7:161749/1‑149 (MQ=255)
                                                                                                                             ccaCGCCCAGGCGCCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGcc                 >  7:68941/3‑149 (MQ=255)
                                                                                                                                              gcCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCAt                                                      >  7:290666/1‑97 (MQ=255)
                                                                                                                                              gcCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCAt                                                      <  8:290666/97‑1 (MQ=255)
                                                                                                                                              gcCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTACgg  <  4:118520/149‑1 (MQ=255)
                                                                                                                                              gcCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTACgg  <  4:233925/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATGAACCGGGCGAGCACCACCGCGTTGACGGGTGGCTCATGAGCAAGCAATCACCAAGAAAAACGAGACCTCCCGCCCAATCGCCTCTACGG  >  NC_002947/4936490‑4936778

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: