Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 715,009:1 +T 11.3% 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.T11.3% 214.0 / 25.5 62intergenic (+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 . (29/26);  new base T (4/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGC  >  NC_002947/714869‑715155
                                                                                                                                              |                                                                                                                                                   
gACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTagcag                                                                                                                                             >  1:46048/1‑149 (MQ=255)
gACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTagcg                                                                                                                                              >  6:230921/1‑148 (MQ=255)
  cGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGCGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCg                                                                                                                                          >  1:255626/1‑149 (MQ=255)
     gcgcGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGt                                                                                                                                        >  3:348908/1‑148 (MQ=255)
      cgcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCa                                                                                                                                      <  5:181553/149‑1 (MQ=255)
        cgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGCGCGGTCAtt                                                                                                                                    <  5:269161/149‑1 (MQ=255)
           agagAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCt                                                                                                                                 <  8:230880/149‑1 (MQ=255)
             agaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCa                                                                                                                               <  3:325700/149‑1 (MQ=255)
              gaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTAGCGCGGTCATTGcagga                                                                                                                              <  5:230921/149‑5 (MQ=255)
              gaACATCACCTTATCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAg                                                                                                                              <  7:63015/149‑1 (MQ=255)
                   tCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcct                                                                                                                         >  7:306998/1‑149 (MQ=255)
                   tCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcct                                                                                                                         <  2:46048/149‑1 (MQ=255)
                    cACCTTTTCAATGCCCAAAGAGCTGCTGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTAGCGCGGTCAtt                                                                                                                                    <  5:194068/137‑1 (MQ=255)
                    cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAATGACCGCTGCTAGCGCGGTCAtt                                                                                                                                    >  6:194068/1‑137 (MQ=255)
                          ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                   <  6:162495/148‑1 (MQ=255)
                          ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                   <  7:329760/148‑1 (MQ=255)
                                 ccgatagagGTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCTGTCATTGCTCAGGTCCTCCGTAAGc                                                                                                                 <  8:610/138‑1 (MQ=255)
                                 cccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGTGGCTGCTGTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGATCAGGTCCTCCGTAAGc                                                                                                                 >  7:610/1‑143 (MQ=255)
                                  ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGTCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                          <  8:183061/149‑1 (MQ=255)
                                   cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                         >  6:360816/1‑149 (MQ=255)
                                    aaaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                         >  5:180201/1‑148 (MQ=255)
                                     aaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGActc                                                                                                       <  7:303279/149‑1 (MQ=255)
                                            gCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCt                                                                                                <  4:3380/149‑1 (MQ=255)
                                             cgtgTTCTCTGGCTGATCTGACGGGCAGTACCGCGGACAATGCAGTGGCTGCTTTGGAGCCGTTGCGAGACACAAGACCTCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAATCATCCTGACTCTGCCGCtt                                                                                               <  4:348908/146‑1 (MQ=255)
                                             cGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGc                                                                                                 <  2:255626/147‑1 (MQ=255)
                                             cGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCtt                                                                                               <  3:199167/149‑1 (MQ=255)
                                                    gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACAAAAGACCGCTCCTGCAATGACCCCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGc                                                                                          >  8:151279/1‑147 (MQ=255)
                                                     ctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGCGGCTGCTTTGCAGCCCTTGCGCGACACAACACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAACCAACCTCACTCTGCCGCTTGATCCCac                                                                                       >  8:291972/1‑148 (MQ=255)
                                                        gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGca                                                                                     <  7:291972/148‑1 (MQ=255)
                                                                  cGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGatat                                                                          >  6:140626/1‑149 (MQ=255)
                                                                         tACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                   <  6:180201/149‑1 (MQ=255)
                                                                         tACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCCGc                                                                   >  4:176040/1‑149 (MQ=255)
                                                                           ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcaccac                                                                                 <  8:220725/133‑1 (MQ=255)
                                                                           ccgccgCCAATGCAGGGGCTGCTTGGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcaccac                                                                                 >  7:220725/1‑133 (MQ=255)
                                                                            cgccgcCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                          <  1:302283/107‑1 (MQ=255)
                                                                            cgccgcCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                          >  2:302283/1‑107 (MQ=255)
                                                                                 ccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAgcg                                                           <  5:360816/149‑1 (MQ=255)
                                                                                   aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTGCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCCGCACCTGATATCGGCTGCGAGCCgcggc                                                         >  2:45645/1‑149 (MQ=255)
                                                                                   aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCt                                                                     >  8:69807/1‑137 (MQ=255)
                                                                                   aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCt                                                                     <  7:69807/137‑1 (MQ=255)
                                                                                        aagggCTGCAAAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCTGGACATTGCTCAGGTCCTCCGTAAGCAACATGACTCTGCCGCTTGATGCCAGCACAACCTGATATCGGATGCGAGCAGCGGCGCGt                                                     >  1:286289/3‑149 (MQ=255)
                                                                                          gggCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCg                                                  >  7:71173/1‑149 (MQ=255)
                                                                                                   tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATAGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCACCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTc                                          >  1:239654/1‑149 (MQ=255)
                                                                                                    gCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                       >  6:225660/1‑118 (MQ=255)
                                                                                                    gCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                       <  5:225660/118‑1 (MQ=255)
                                                                                                    gCAGCCCTAGCGCGACACAAGCCCGCTAATGCAAAGACCGCGCGAGCAGCGGTCATTGCTCAGGGCCTCCGTAACCAACCTGCCTCTCCCCCTAGATCCcaccac                                                                                    >  8:153745/1‑105 (MQ=255)
                                                                                                          cTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGc                                                 >  1:282972/1‑134 (MQ=255)
                                                                                                          cTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGc                                                 <  2:282972/134‑1 (MQ=255)
                                                                                                          cTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTg                                  >  5:175718/1‑149 (MQ=255)
                                                                                                           ttGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCGTGATGCCAGCCCCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgc                                 >  7:168532/1‑149 (MQ=255)
                                                                                                             gcgcGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgcg                                <  8:306998/148‑1 (MQ=255)
                                                                                                             gcgcGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgc                                 <  2:239654/147‑1 (MQ=255)
                                                                                                                    acaAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAAc                          >  2:274110/1‑147 (MQ=255)
                                                                                                                           cGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTg                 >  2:164546/1‑149 (MQ=255)
                                                                                                                              tCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc              <  8:71173/149‑1 (MQ=255)
                                                                                                                              ctgagCAATGACCGCTGCTAGCGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc              >  5:345322/5‑149 (MQ=255)
                                                                                                                                    aaTGACCGCTGCTAGCGCGGTCATTGCTCAGGTCCTCCGAAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGt         >  4:142695/1‑148 (MQ=255)
                                                                                                                                    aaTGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAAc                                                                                                              >  4:136198/1‑47 (MQ=255)
                                                                                                                                    aaTGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAAc                                                                                                              <  3:136198/47‑1 (MQ=255)
                                                                                                                                          ccGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATggg   >  3:171810/1‑149 (MQ=255)
                                                                                                                                          ccGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATggg   >  8:67305/1‑149 (MQ=255)
                                                                                                                                          ccGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATggg   >  8:72187/1‑149 (MQ=255)
                                                                                                                                           cgcgCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGc  <  4:171810/149‑1 (MQ=255)
                                                                                                                                           cgcgCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGc  <  2:286289/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGC  >  NC_002947/714869‑715155

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