Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,764,801 A→G 7.8% D261G (GAC→GGC)  ctpA → carboxy‑terminal‑processing protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,764,8010AG 7.8% 252.0 / 9.2 90D261G (GAC→GGC) ctpAcarboxy‑terminal‑processing protease
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base G (5/2);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

TCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAT  >  NC_002947/5764654‑5764948
                                                                                                                                                   |                                                                                                                                                   
tCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGAc                                                                                                                                                    <  7:185446/149‑1 (MQ=255)
  aCCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCa                                                                                                                                                  >  8:185352/1‑149 (MQ=255)
  aCCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCa                                                                                                                                                  >  6:160086/1‑149 (MQ=255)
   cccAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCTATCGGCAGTGGAAGTGCCCGACCAc                                                                                                                                                 >  2:218776/1‑149 (MQ=255)
      aGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTc                                                                                                                                              <  8:145428/149‑1 (MQ=255)
          ccAGGTCCAGACCGGTGACGAAGTGGTCAAGGCCCTGGCCAAGCTGAGCAAGGACAACGGAAAGAAGCTGCGCGGGTTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCTACCACTTCCt                                                                                                                                            <  6:196276/147‑1 (MQ=255)
          ccAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccac                                                                                                                                           <  5:291677/148‑3 (MQ=255)
              gTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccactgccg                                                                                                                                      >  6:235092/1‑142 (MQ=255)
                 aaGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAgg                                                                                                                                   >  2:277570/1‑149 (MQ=255)
                 aaGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAgg                                                                                                                                   >  2:203563/1‑149 (MQ=255)
                 aaGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAgg                                                                                                                                   >  4:48447/1‑149 (MQ=255)
                   gACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTc                                                                                                                                 <  7:85919/149‑1 (MQ=255)
                   gACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTc                                                                                                                                 >  3:199135/1‑149 (MQ=255)
                      cGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGa                                                                                                                              <  6:249021/149‑1 (MQ=255)
                         tGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGAGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCg                                                                                                                           >  8:218563/1‑149 (MQ=255)
                           aCGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTc                                                                                                                         <  1:203563/149‑1 (MQ=255)
                           aCGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTc                                                                                                                         >  3:21812/1‑149 (MQ=255)
                             gAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                       >  6:322252/1‑149 (MQ=255)
                              aaGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                       >  6:134775/1‑148 (MQ=255)
                              aaGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                       >  3:307279/1‑148 (MQ=255)
                              aaGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                       >  3:115494/1‑148 (MQ=255)
                               aGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCTGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                                     <  1:30980/149‑1 (MQ=255)
                               aGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                                     <  4:115494/149‑1 (MQ=255)
                                gTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTacac                                                                                                                    <  7:185352/149‑1 (MQ=255)
                                 tGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTc                                                                                                                                              >  7:318838/1‑122 (MQ=255)
                                 tGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTc                                                                                                                                              <  8:318838/122‑1 (MQ=255)
                                    gCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCaaa                                                                                                                <  5:322252/149‑1 (MQ=255)
                                    gCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCa                                                                                                                  >  2:59172/1‑147 (MQ=255)
                                     cAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAg                                                                                                               <  3:48447/149‑1 (MQ=255)
                                     cAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAg                                                                                                               <  7:258422/149‑1 (MQ=255)
                                     cAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAg                                                                                                               <  1:277570/149‑1 (MQ=255)
                                          cccTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCg                                                                                                           >  5:15407/1‑148 (MQ=255)
                                             tGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCATAGGTCTGATCGTCTACACCAAAGGcc                                                                                                            >  4:247845/1‑144 (MQ=255)
                                             tGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGcc                                                                                                            <  3:247845/144‑1 (MQ=255)
                                               gCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATc                                                                                                       <  4:109197/147‑1 (MQ=255)
                                               gCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATc                                                                                                       <  4:21812/147‑1 (MQ=255)
                                               gCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGc                                                                                                     >  8:221361/1‑149 (MQ=255)
                                                 cAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGcc                                                                                                    >  1:277292/1‑148 (MQ=255)
                                                  aaGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCaa                                                                                                  >  5:199191/1‑149 (MQ=255)
                                                        ctcAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCTACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCg                                                                                             <  7:218563/146‑1 (MQ=255)
                                                        cgcAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACATCAAAGGCCGCATCGCCAACTcc                                                                                              >  1:23431/1‑147 (MQ=255)
                                                           aaGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                       <  2:188334/119‑1 (MQ=255)
                                                           aaGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCt                                                                                         >  4:348283/1‑149 (MQ=255)
                                                           aaGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGAGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                       >  1:188334/1‑119 (MQ=255)
                                                                    ggCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTtctc                                                                                <  8:305/149‑1 (MQ=255)
                                                                     gCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCg                                                                               <  6:15407/149‑1 (MQ=255)
                                                                         gaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGcc                                                                            >  8:213822/1‑148 (MQ=255)
                                                                            gCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCAACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGAcc                                                                        >  1:308005/1‑149 (MQ=255)
                                                                                cgcgGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGc                                                                    <  2:277292/149‑1 (MQ=255)
                                                                                 gcgGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGTGCTGCGTTTCTCGGCCGACCCGGcc                                                                   <  3:53719/149‑1 (MQ=255)
                                                                                   ggggTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGa                                                                 <  5:134775/149‑1 (MQ=255)
                                                                                       tggtgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGc                                                              >  2:248650/1‑148 (MQ=255)
                                                                                         gtgCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTTGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCa                                                            <  4:39317/148‑1 (MQ=255)
                                                                                           gCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAgcg                                                         >  6:95518/1‑149 (MQ=255)
                                                                                                aCCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGc                                                    <  4:199135/149‑1 (MQ=255)
                                                                                                aCCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGc                                                    >  7:17615/1‑149 (MQ=255)
                                                                                                     cgcAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACATCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTg                                                 <  2:23431/147‑1 (MQ=255)
                                                                                                        aacaacCCAGGCGGCGTGCTGCAATCGGCATTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCt                                            >  1:34507/1‑149 (MQ=255)
                                                                                                        aacaacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCt                                            >  2:94283/1‑149 (MQ=255)
                                                                                                           aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggt                                         <  3:280347/149‑1 (MQ=255)
                                                                                                           aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggt                                         >  3:257874/1‑149 (MQ=255)
                                                                                                            acCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                        <  2:34507/149‑1 (MQ=255)
                                                                                                            acCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                        <  5:160086/149‑1 (MQ=255)
                                                                                                               cAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTa                                     >  2:206124/1‑149 (MQ=255)
                                                                                                               cAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTa                                     <  1:248650/149‑1 (MQ=255)
                                                                                                                aGGCGGCGTGCTTCAATCGGCAGTGTAAGTTGCAGACCAATTCCTGACAAAAGGTCTGATCTTCTACACCAAAGTCCGCATCGCCAACTCCGAGCTTCGTTTCTCGGCCGACCCGGCCGATGCCATCGCGGGAGTGCCGCTGGTGGTAc                                    <  6:202414/149‑1 (MQ=255)
                                                                                                                aGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTAc                                    >  5:246704/1‑149 (MQ=255)
                                                                                                                   cggcggGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGa                                 <  8:135959/149‑1 (MQ=255)
                                                                                                                   cggcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGa                                 >  5:309463/1‑149 (MQ=255)
                                                                                                                   cggcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGa                                 <  7:213822/149‑1 (MQ=255)
                                                                                                                     gcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATc                               <  1:337525/149‑1 (MQ=255)
                                                                                                                     gcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATc                               <  1:59172/149‑1 (MQ=255)
                                                                                                                       gTGCTGCAATCGGCGGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCaa                             >  5:365456/1‑149 (MQ=255)
                                                                                                                        tgctgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAAc                            <  3:83685/149‑1 (MQ=255)
                                                                                                                          ctgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACgg                          <  5:79266/149‑1 (MQ=255)
                                                                                                                          ctgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACg                           <  7:221361/148‑1 (MQ=255)
                                                                                                                          ctgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACg                           <  3:348283/148‑1 (MQ=255)
                                                                                                                           tgcAATCGGCAGTGGAAGTGGCCGACCACTTCCAGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTa                                     >  2:129858/1‑137 (MQ=255)
                                                                                                                               aTCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCa                     <  6:246704/149‑1 (MQ=255)
                                                                                                                                     tcaggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcat               >  5:79426/4‑149 (MQ=255)
                                                                                                                                     tcaggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcat               >  1:4475/4‑149 (MQ=255)
                                                                                                                                      caggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGCCCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGCTCAACGGTGGCAGTgcatc              >  4:244463/3‑149 (MQ=255)
                                                                                                                                      caggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatc              >  6:90623/3‑149 (MQ=255)
                                                                                                                                       tGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcg             <  1:206124/149‑1 (MQ=255)
                                                                                                                                          aaGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCTACAACGAGCTGCGTTTCTCGGCCGACCCCGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcggca          >  1:298347/1‑149 (MQ=255)
                                                                                                                                               gCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAACGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGcatcggcatcgg      >  8:161069/1‑148 (MQ=255)
                                                                                                                                               gCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGCGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGa     <  7:301479/149‑1 (MQ=255)
                                                                                                                                               gCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGt                         >  5:354272/1‑129 (MQ=255)
                                                                                                                                               gCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCAGCATCTCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGt                         <  6:354272/129‑1 (MQ=255)
                                                                                                                                                  gACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAt  >  4:21988/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAT  >  NC_002947/5764654‑5764948

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: