Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 715,009:1 +T 7.1% intergenic (+75/+2) ispH → / ← PP_0607 4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947715,0091.T 7.1% 219.3 / 8.0 56intergenic (+75/+2)ispH/PP_06074‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base . (26/26);  new base T (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCAC  >  NC_002947/714863‑715157
                                                                                                                                                   |                                                                                                                                                     
gagCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCgcgc                                                                                                                                                     <  6:194652/149‑1 (MQ=255)
      gACGGGCGCGAAGAGAACATCACCTTTTCAATGACCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCTCGACACAAGACCGCTCCTGCAATGACCGCGCTagcag                                                                                                                                               >  8:176953/1‑149 (MQ=255)
          ggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCTGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGt                                                                                                                                           <  3:231670/149‑1 (MQ=255)
          ggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGt                                                                                                                                           <  5:242253/149‑1 (MQ=255)
          ggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGt                                                                                                                                           <  2:339593/149‑1 (MQ=255)
                    gaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAg                                                                                                                                 <  6:120441/149‑1 (MQ=255)
                     aaCATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGCGCGGTCATTGcaggag                                                                                                                                <  5:134309/149‑6 (MQ=255)
                      aCATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAgg                                                                                                                                <  8:26872/148‑1 (MQ=255)
                       cATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtc                                                                                                                              >  4:358827/1‑149 (MQ=255)
                        aTCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcc                                                                                                                             >  8:275520/1‑149 (MQ=255)
                        aTCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAACACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcc                                                                                                                             >  7:79507/1‑149 (MQ=255)
                                ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTa                                                                                                                       <  8:174082/147‑1 (MQ=255)
                                   aaTGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCaa                                                                                                                 <  6:117390/149‑1 (MQ=255)
                                    aTGCCCAAAGAGCTGCTGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAAc                                                                                                                <  7:176953/149‑1 (MQ=255)
                                         cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                      >  3:278727/1‑139 (MQ=255)
                                         cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                      <  4:278727/139‑1 (MQ=255)
                                         cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                           >  8:202034/1‑149 (MQ=255)
                                                 tGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGc                                                                                                   >  2:272010/1‑149 (MQ=255)
                                                      gTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGAt                                                                                              <  7:275520/149‑1 (MQ=255)
                                                          gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGc                                                                                            >  2:17821/1‑147 (MQ=255)
                                                              gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGca                                                                                       <  1:38529/148‑1 (MQ=255)
                                                                     tGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGa                                                                               >  5:42191/1‑149 (MQ=255)
                                                                       aCGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGata                                                                             >  6:262742/1‑149 (MQ=255)
                                                                           gCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGata                                                                             >  7:275864/1‑145 (MQ=255)
                                                                           gCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGata                                                                             <  8:275864/145‑1 (MQ=255)
                                                                           gCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                         >  7:306071/1‑149 (MQ=255)
                                                                               tACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCATAAGCATCCTGACTCTGCCGCTAGATGCCAGCACCACCTGATATCGGCTGc                                                                     >  7:225652/1‑149 (MQ=255)
                                                                                     cgccAAAGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagca                                                                <  1:272010/149‑1 (MQ=255)
                                                                                        cAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGaacgcgcc                                                            >  3:315322/1‑142 (MQ=255)
                                                                                             caagggCTGCAAAGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGt                                                       >  1:87855/4‑149 (MQ=255)
                                                                                               ggggCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTc                                                     <  7:202034/149‑1 (MQ=255)
                                                                                                  gctgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGca                                                  >  2:47218/1‑149 (MQ=255)
                                                                                                    tgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGcaca                                                >  4:283951/1‑149 (MQ=255)
                                                                                                     gctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCaccacc                                                                                  >  7:218598/1‑114 (MQ=255)
                                                                                                     gctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCaccacc                                                                                  <  8:218598/114‑1 (MQ=255)
                                                                                                          gCAGCCCTTGCGCGACACAAGAACGCTCCTGAAATGACCGCGCTAGCAGCGGTCATTTCTCAGGTCCTCCGAAAGCAACCTGACTCTGCCGCTTGAAGCCAGCACCACCAGATATCGGCTGCGAGCACCGGCGCGTTCGCACATTTcc                                           >  6:71673/1‑148 (MQ=255)
                                                                                                          gCAGCCCTTGCGCGACACAACACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTcc                                           <  8:79507/148‑1 (MQ=255)
                                                                                                             gCCCTTGCGCGACACAAGACCGCTCCTGCAATGACGGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTATGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCTGtt                                       <  7:326121/149‑1 (MQ=255)
                                                                                                               ccTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGtttt                                     <  8:306071/149‑1 (MQ=255)
                                                                                                               ccTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGtttt                                     <  4:260063/149‑1 (MQ=255)
                                                                                                                   gcgcGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgc                                   <  1:17821/147‑1 (MQ=255)
                                                                                                                    cgcgACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGtt                                                      <  6:343202/127‑1 (MQ=255)
                                                                                                                    cgcgACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGtt                                                      >  5:343202/1‑127 (MQ=255)
                                                                                                                             aGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCt                       <  2:87855/149‑1 (MQ=255)
                                                                                                                                 cGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTg                   >  7:221411/1‑149 (MQ=255)
                                                                                                                                         cAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGAc             <  3:283951/147‑1 (MQ=255)
                                                                                                                                         cAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGt           >  6:172662/1‑149 (MQ=255)
                                                                                                                                          aaTGACCGCTGCTAGCGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGt           >  2:129019/1‑148 (MQ=255)
                                                                                                                                          aaTGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAAc                            >  5:188774/1‑131 (MQ=255)
                                                                                                                                          aaTGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAAc                            <  6:188774/131‑1 (MQ=255)
                                                                                                                                          aaTGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGt           >  2:125261/1‑148 (MQ=255)
                                                                                                                                          aaTGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGt           >  8:70668/1‑148 (MQ=255)
                                                                                                                                           aTGACCGCTGCTAGCGCGGTCATTGCTCTGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCa         >  1:204588/1‑149 (MQ=255)
                                                                                                                                            tGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCAt        <  3:333941/149‑1 (MQ=255)
                                                                                                                                            tGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCAt        >  2:399/1‑149 (MQ=255)
                                                                                                                                             gACCGCTGCTAGCGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATg       <  1:129019/149‑1 (MQ=255)
                                                                                                                                              aCCGCTCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATgg      <  1:125261/149‑1 (MQ=255)
                                                                                                                                                cgcgCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATgg      >  7:225751/1‑147 (MQ=255)
                                                                                                                                                cgcgCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATgg      <  8:225751/147‑1 (MQ=255)
                                                                                                                                                  cgcTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCAc  <  6:42191/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCAC  >  NC_002947/714863‑715157

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: