Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,398,163 G→C 7.6% I498M (ATC→ATG) ‡ PP_2100 ← two‑component system sensor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,398,1630GC 7.6% 279.2 / 8.1 89I498M (ATC→ATG) ‡PP_2100two‑component system sensor protein
Reads supporting (aligned to +/- strand):  ref base G (42/40);  new base C (4/3);  total (46/43)
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

ACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGGTGGTGATCACCGCATCGACAG  >  NC_002947/2398015‑2398308
                                                                                                                                                    |                                                                                                                                                 
acaaAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAACGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAgg                                                                                                                                                   >  8:192731/1‑149 (MQ=255)
   aaaGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATc                                                                                                                                                >  1:315246/1‑149 (MQ=255)
       cACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGac                                                                                                                                             >  2:86868/1‑148 (MQ=255)
          ggCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGacaaca                                                                                                                                         <  7:329520/149‑1 (MQ=255)
                ggCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCg                                                                                                                                   <  2:315246/149‑1 (MQ=255)
                 gCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCgg                                                                                                                                  <  1:140354/149‑1 (MQ=255)
                   cAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGATCCTCGacaaag                                                                                                                                         >  1:187539/1‑138 (MQ=255)
                   cAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGATCCTCGacaaag                                                                                                                                         <  2:187539/140‑3 (MQ=255)
                   cAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGa                                                                                                                                              >  8:138102/1‑135 (MQ=255)
                   cAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGa                                                                                                                                              <  7:138102/135‑1 (MQ=255)
                   cAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGt                                                                                                                                 >  6:138926/1‑148 (MQ=255)
                   cAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTATTCCGGTGAGATGCTCAGCATAGCTTTGTCGAGGATCGAGAACAGAGCCGGTg                                                                                                                                >  8:153844/1‑149 (MQ=255)
                      caccacCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGATCCTCGacaaagcttt                                                                                                                                     >  7:210361/1‑135 (MQ=255)
                      caccacCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGATCCTCGacaaagcttt                                                                                                                                     <  8:210361/141‑7 (MQ=255)
                       acctccTGCATGGCCAGGAGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAACTCCACGCCATCGGTTGTTTATTACCCCCAGCTCTTCCGGTGAGAGTCTCAGCAAAGCTTTGTCTAGGATCGACAACAGCGCCGGTGActg                                                                                                                            <  8:115723/149‑3 (MQ=255)
                       accaccTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACgg                                                                                                                            <  4:41770/149‑1 (MQ=255)
                        ccaccTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCGGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACggg                                                                                                                           <  3:108787/149‑1 (MQ=255)
                                   gCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATc                                                                                                                                                >  5:113571/1‑117 (MQ=255)
                                   gCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATc                                                                                                                                                <  6:113571/117‑1 (MQ=255)
                                     cAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCa                                                                                                              >  5:228386/1‑149 (MQ=255)
                                        gcggcTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATgg                                                                                                           >  4:82129/1‑149 (MQ=255)
                                             tGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAg                                                                                                      >  6:20346/1‑149 (MQ=255)
                                             tGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAg                                                                                                      >  4:128657/1‑149 (MQ=255)
                                               tAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGtt                                                                                                    >  5:79210/1‑149 (MQ=255)
                                                 gCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTgg                                                                                                  >  7:135143/1‑149 (MQ=255)
                                                  ccTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGc                                                                                                 <  3:192810/149‑1 (MQ=255)
                                                  ccTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGc                                                                                                 <  7:192731/149‑1 (MQ=255)
                                                   cTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCg                                                                                                >  8:144165/1‑149 (MQ=255)
                                                    tttGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGGTGCCATGGCGAAGTTGGCgg                                                                                               >  4:64272/1‑149 (MQ=255)
                                                       ggAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGtt                                                                                            >  7:14946/1‑149 (MQ=255)
                                                        gAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTc                                                                                           <  7:253643/149‑1 (MQ=255)
                                                           ttGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCt                                                                                        <  2:194372/149‑1 (MQ=255)
                                                              tCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACACCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGAGGCCATGGAGAAGTTGGCGGGTTCGCtcc                                                                                      >  3:46734/1‑146 (MQ=255)
                                                               cATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGcc                                                                                     >  4:114162/1‑148 (MQ=255)
                                                                    tCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATg                                                                                                            <  5:238041/120‑1 (MQ=255)
                                                                    tCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATg                                                                                                            >  6:238041/1‑120 (MQ=255)
                                                                    tCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCgg                                                                                               >  8:107829/1‑133 (MQ=255)
                                                                    tCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCgg                                                                                               <  7:107829/133‑1 (MQ=255)
                                                                    tCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGccgacc                                                                                 >  8:42286/1‑147 (MQ=255)
                                                                     cATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAg                                                                              <  7:144165/149‑1 (MQ=255)
                                                                      aTGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAg                                                                              <  6:182291/148‑1 (MQ=255)
                                                                          ctgccAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCg                                                                         >  4:298274/1‑149 (MQ=255)
                                                                           tgccAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCg                                                                         >  3:236957/1‑148 (MQ=255)
                                                                            gccAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGa                                                                       >  4:199301/1‑149 (MQ=255)
                                                                             ccAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGAt                                                                      <  8:135143/149‑1 (MQ=255)
                                                                             ccAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGAt                                                                      <  3:82129/149‑1 (MQ=255)
                                                                               aaaGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTg                                                                                                                                >  6:154230/1‑89 (MQ=255)
                                                                                aaGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGAc                                                                   <  1:86868/149‑1 (MQ=255)
                                                                                  gCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCa                                                                 <  3:114162/149‑1 (MQ=255)
                                                                                          aGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTggg                                                         >  8:29725/1‑149 (MQ=255)
                                                                                          aGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGGTGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTggg                                                         <  3:64272/149‑1 (MQ=255)
                                                                                           gCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGt                                                        <  4:36130/149‑1 (MQ=255)
                                                                                             ggTTGTTGATTACCCCCAGCTCTTCCGGGGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTg                                                      <  5:249993/149‑1 (MQ=255)
                                                                                               ttgttgATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCt                                                    <  6:228386/149‑1 (MQ=255)
                                                                                               ttgttgATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCt                                                    <  3:128657/149‑1 (MQ=255)
                                                                                                 gttgATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGTGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGAACGAGGTGCGGATGACTTGATATGGGtttattg                                                  >  1:304429/1‑144 (MQ=255)
                                                                                                        aCCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAggg                                            <  3:199301/148‑1 (MQ=255)
                                                                                                         cccccAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACAGGGCAATGGTGGCCATGGCGAAGTTGGCgg                                                                                               <  5:186478/96‑1 (MQ=255)
                                                                                                         cccccAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACAGGGCAATGGTGGCCATGGCGAAGTTGGCgg                                                                                               >  6:186478/1‑96 (MQ=255)
                                                                                                          ccccAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGcc                                                                                     <  3:193309/105‑1 (MQ=255)
                                                                                                          ccccAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGcc                                                                                     >  4:193309/1‑105 (MQ=255)
                                                                                                          ccccAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTATGGGTTGCTGGCAAGGTCGt                                         >  3:226011/1‑149 (MQ=255)
                                                                                                           cccAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAggg                                            >  2:266721/1‑145 (MQ=255)
                                                                                                           cccAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAggg                                            <  1:266721/145‑1 (MQ=255)
                                                                                                              aGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCa                                                                 <  3:32894/121‑1 (MQ=255)
                                                                                                              aGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCa                                                                 >  4:32894/1‑121 (MQ=255)
                                                                                                               gCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTg                                    <  8:14946/149‑1 (MQ=255)
                                                                                                                  ctTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCg                                 >  7:292236/1‑149 (MQ=255)
                                                                                                                     ccGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTc                               <  5:138926/148‑1 (MQ=255)
                                                                                                                        gTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGAt                            >  6:241354/1‑148 (MQ=255)
                                                                                                                        gTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTGGCTGGCCAGGGCGTGATTGGCGTCGAt                            >  8:109965/1‑148 (MQ=255)
                                                                                                                         tGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCa                          <  7:42286/149‑1 (MQ=255)
                                                                                                                           agaTGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAgg                        <  4:236957/149‑1 (MQ=255)
                                                                                                                               gCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTg                                        >  7:91799/1‑129 (MQ=255)
                                                                                                                               gCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTg                                        <  8:91799/129‑1 (MQ=255)
                                                                                                                                   aGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGGTGGTGAt                >  3:42620/1‑149 (MQ=255)
                                                                                                                                     cAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGa                                       <  2:82264/124‑1 (MQ=255)
                                                                                                                                      aaaGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGa                             >  4:86615/1‑133 (MQ=255)
                                                                                                                                      aaaGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGa                             <  3:86615/133‑1 (MQ=255)
                                                                                                                                      aaaGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGgtggtg                  >  2:134257/1‑144 (MQ=255)
                                                                                                                                      aaaGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGgtggtg                  <  1:134257/144‑1 (MQ=255)
                                                                                                                                      aaaGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTCATTGGCGTCGATCAGGGAGGTGAGCAc             >  4:236734/1‑149 (MQ=255)
                                                                                                                                        gctgttGTCGATCCTCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAgg                                             >  2:167447/5‑115 (MQ=255)
                                                                                                                                        gctgttGTCGATCCTCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAgg                                             <  1:167447/111‑1 (MQ=255)
                                                                                                                                         ctgttGTCGATCCTCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGGTGGTGATCACCg           >  8:154913/4‑148 (MQ=255)
                                                                                                                                             tGTCGAGGAGCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGGTGGTGATCACCGCAt        <  5:20346/147‑1 (MQ=255)
                                                                                                                                                 gAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGgtggtg                  <  3:87230/133‑1 (MQ=255)
                                                                                                                                                 gAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGgtggtg                  >  4:87230/1‑133 (MQ=255)
                                                                                                                                                 gAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGGTGGTGATCACCGCATCGACAg  >  5:35511/1‑149 (MQ=255)
                                                                                                                                                 gAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGGTGGTGATCACCGCATCGACAg  >  5:77707/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTGTCATCCTCATGCCTGCCAAAGCCACGCCAGCGGTTGTTGATTACCCCCAGCTCTTCCGGTGAGATGCTCAGCAAAGCTTTGTCGAGGATCGACAACAGCGCCGGTGACGGGGCAATGGTGGCCATGGCGAAGTTGGCGGGTTCGCTGCCGACCGAGGTGCGGATGACCAGTTCTGGGTTGCTGGCCAGGGCGTGATTGGCGTCGATCAGGGTGGTGATCACCGCATCGACAG  >  NC_002947/2398015‑2398308

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: