Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,690,884 A→C 6.4% D75A (GAC→GCC)  PP_2357 → type I pili subunit CsuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,690,8840AC 6.4% 184.8 / 3.9 63D75A (GAC→GCC) PP_2357type I pili subunit CsuB
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base C (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

CTGGCCGCGCTGTTCGGCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCATGGTGGCATCGGCGA  >  NC_002947/2690742‑2691028
                                                                                                                                               |                                                                                                                                                
cTGGCCGCGCTGTTCGGCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAgg                                                                                                                                            <  4:318111/149‑1 (MQ=255)
   gCCGCGCTGTTCGGCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAggg                                                                                                                                           <  3:187734/147‑1 (MQ=255)
     cgcgCTGTTCGGCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGa                                                                                                                                       <  8:83646/149‑1 (MQ=255)
           gTTCCGCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGgcgc                                                                                                                                 >  3:171633/1‑149 (MQ=255)
            ttCGGCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGACTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGgcgcg                                                                                                                                <  4:356326/149‑1 (MQ=255)
                gCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTg                                                                                                                            >  1:59903/1‑149 (MQ=255)
                gCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTg                                                                                                                            >  5:75793/1‑149 (MQ=255)
                    gttgttCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATg                                                                                                                        >  6:350529/1‑149 (MQ=255)
                       gttCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGcc                                                                                                                     >  2:199355/1‑149 (MQ=255)
                       gttCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGcc                                                                                                                     >  2:397929/1‑149 (MQ=255)
                           aGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTgg                                                                                                                 >  2:121288/1‑149 (MQ=255)
                           aGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGAAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGAAGCGGGAGCTCAGGCCCTGGGACGTATCGACCGAGGGGCGACGGCGCGCCTGGATGGCCCTgg                                                                                                                 >  5:181209/1‑149 (MQ=255)
                              ggCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGgcgc                                                                                                              <  1:13949/149‑1 (MQ=255)
                                    aCCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCtt                                                                                                        >  5:131791/1‑149 (MQ=255)
                                      cGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGa                                                                                                      <  6:75793/149‑1 (MQ=255)
                                      cGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGa                                                                                                      <  7:185947/149‑1 (MQ=255)
                                      cGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGa                                                                                                      <  2:59903/149‑1 (MQ=255)
                                              cggATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAggg                                                                                                                                           >  1:104131/1‑104 (MQ=255)
                                              cggATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAggg                                                                                                                                           <  2:104131/104‑1 (MQ=255)
                                                  ttttCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGt                                                                                          >  2:283743/1‑149 (MQ=255)
                                                     tctcGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGtt                                                                                         >  4:367093/1‑147 (MQ=255)
                                                       tcGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGc                                                                                     >  4:266256/1‑149 (MQ=255)
                                                         gTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGcc                                                                                    <  6:308805/148‑1 (MQ=255)
                                                         gTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGc                                                                                     >  2:301820/1‑147 (MQ=255)
                                                         gTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCa                                                                                   <  3:266256/149‑1 (MQ=255)
                                                         gTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGAGGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAg                                                                                      <  1:283743/146‑1 (MQ=255)
                                                          tGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAg                                                                                  >  6:256245/1‑149 (MQ=255)
                                                           ggAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGc                                                                                 <  1:199355/149‑1 (MQ=255)
                                                               gTGCGGGTACTGGTGCAGCGTGGCTGCATGGCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTc                                                                              >  1:260564/1‑149 (MQ=255)
                                                                tGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTccg                                                                            <  5:31878/149‑1 (MQ=255)
                                                                 gCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTccgc                                                                           >  2:308000/1‑149 (MQ=255)
                                                                  cGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTccgcc                                                                          <  4:323170/149‑1 (MQ=255)
                                                                     gTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCcgcg                                                                       >  4:255184/1‑149 (MQ=255)
                                                                        cTGGTGCAGCGTCGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCt                                                                    >  6:90331/1‑149 (MQ=255)
                                                                         tGGTGCAGCGTGGCTGCAT‑GCTGGTCAAGCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTTTCGTCCTTGGGGCGACGGCGCGCCTGGATTGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCTGCGTCCGGCGCGCCtt                                                                   <  4:87462/149‑2 (MQ=255)
                                                                           gTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTgg                                                                  <  2:260564/148‑1 (MQ=255)
                                                                           gTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTg                                                                   >  7:164739/1‑147 (MQ=255)
                                                                           gTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGa                                                                 <  5:350529/149‑1 (MQ=255)
                                                                             gCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAgtg                                                                                                   <  6:98933/113‑1 (MQ=255)
                                                                              cAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTg                                                              <  1:397929/149‑1 (MQ=255)
                                                                                   tGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAAcc                                                         >  7:364231/1‑149 (MQ=255)
                                                                                             gCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGc                                                >  5:122387/1‑149 (MQ=255)
                                                                                              cTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGcc                                               <  6:122387/149‑1 (MQ=255)
                                                                                               tGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGccc                                              >  7:14747/1‑149 (MQ=255)
                                                                                               tGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGccc                                              >  8:402902/1‑149 (MQ=255)
                                                                                                 gTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGccc                                              <  5:111969/147‑1 (MQ=255)
                                                                                                 gTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCt                                             <  6:131791/148‑1 (MQ=255)
                                                                                                     aCCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAg                                        <  5:90331/149‑1 (MQ=255)
                                                                                                      ccAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGCCCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAgg                                       >  6:104214/1‑149 (MQ=255)
                                                                                                      ccAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGCCCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAgg                                       >  6:110388/1‑149 (MQ=255)
                                                                                                                     gggAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGAtggtg                        <  3:255184/149‑1 (MQ=255)
                                                                                                                        aGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGcc                     <  1:121288/149‑1 (MQ=255)
                                                                                                                         gCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCa                    <  8:164739/149‑1 (MQ=255)
                                                                                                                             aGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGAt                                                    <  3:259620/113‑1 (MQ=255)
                                                                                                                             aGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGAt                                                    >  4:259620/1‑113 (MQ=255)
                                                                                                                             aGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGc                  <  1:301820/147‑1 (MQ=255)
                                                                                                                             aGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCAt                >  2:243797/1‑149 (MQ=255)
                                                                                                                                 ccTGGGACGTATCGCCCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCAtggtg            <  5:104214/149‑1 (MQ=255)
                                                                                                                                 ccTGGGACGTATCGCCCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGGGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCAtggtg            <  5:110388/149‑1 (MQ=255)
                                                                                                                                      gACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCATGGTGGCATc       <  4:171633/149‑1 (MQ=255)
                                                                                                                                      gACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCATGGTGGCATc       <  8:364231/149‑1 (MQ=255)
                                                                                                                                            tCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCTTATCAGGTGCGCGTCGATGGTGGCCAGCATGGTGGCATCGGCGa  <  1:308000/148‑1 (MQ=255)
                                                                                                                                            tCGACCTAGGGGCGACGGCGCGCATGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCATGGTGGCATCGGCGa  <  5:256245/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CTGGCCGCGCTGTTCGGCCTGTTGTTCAGCGGCAGTACCGTGGCGGCGGATTTTCTCGTGGAAGTGCGGGTACTGGTGCAGCGTGGCTGCAT‑GCTGGTCAACCAGACCCGCGATGCGGGAGCTCAGGCCCTGGGACGTATCGACCTAGGGGCGACGGCGCGCCTGGATGGCCCTGGCGCGCCCTTGAGTGGTGTGCTGTTGAGCCAGCGTCCGCCGCGCCTGGAGTGCAACCCGGATACGCCCTATCAGGTGCGCGTCGATGGTGGCCAGCATGGTGGCATCGGCGA  >  NC_002947/2690742‑2691028

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: ‑