Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 491,091 T→C 5.9% N444D (AAC→GAC) ‡ lptD ← LPS‑assembly protein LptD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947491,0910TC 5.9% 243.5 / 3.5 83N444D (AAC→GAC) ‡lptDLPS‑assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base T (37/41);  new base C (3/2);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGATC  >  NC_002947/490947‑491238
                                                                                                                                                |                                                                                                                                                   
gCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTcag                                                                                                                                                  >  2:186066/1‑148 (MQ=255)
  tGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTcagcag                                                                                                                                               >  1:109575/1‑149 (MQ=255)
     aaCGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGAACCAGCAgggt                                                                                                                                            <  5:50024/149‑3 (MQ=255)
      aCGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGTCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGAACCAGCAgggt                                                                                                                                            <  8:104562/148‑3 (MQ=255)
         gCAGCGTACCGTTCTGGGTATGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCgg                                                                                                                                        <  7:149174/149‑1 (MQ=255)
         gCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCgg                                                                                                                                        <  3:315044/149‑1 (MQ=255)
         gCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCgg                                                                                                                                        <  2:106574/149‑1 (MQ=255)
              gTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTc                                                                                                                                   <  7:186896/149‑1 (MQ=255)
              gTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTc                                                                                                                                   >  6:44844/1‑149 (MQ=255)
               tACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTcc                                                                                                                                  <  1:2292/149‑1 (MQ=255)
                aCCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCa                                                                                                                                 >  3:34186/1‑149 (MQ=255)
                   gTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCTTGCTGGTTCAGCAGGTCGCGGCTTTCCACGc                                                                                                                              >  2:294468/1‑149 (MQ=255)
                        gggTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATc                                                                                                                         <  2:185150/149‑1 (MQ=255)
                          gTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTg                                                                                                                       <  5:13368/149‑1 (MQ=255)
                          gTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCAACCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTg                                                                                                                       >  6:22207/1‑149 (MQ=255)
                             aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                    <  3:267015/149‑1 (MQ=255)
                              gCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCg                                                                                                                   >  4:260106/1‑149 (MQ=255)
                               cTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCg                                                                                                                   >  1:149704/1‑148 (MQ=255)
                               cTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCg                                                                                                                   >  1:145223/1‑148 (MQ=255)
                                tGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                    >  7:41905/1‑146 (MQ=255)
                                tGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                    <  8:41905/146‑1 (MQ=255)
                                 gCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACAATCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGa                                                                                                                 >  3:38425/1‑148 (MQ=255)
                                     cAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCGATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCa                                                                                                            <  7:81561/149‑1 (MQ=255)
                                      aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAg                                                                                                           >  5:139395/1‑149 (MQ=255)
                                      aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAg                                                                                                           <  2:335968/149‑1 (MQ=255)
                                          gATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTcc                                                                                                       >  7:39381/1‑149 (MQ=255)
                                          gATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTcc                                                                                                       >  6:368209/1‑149 (MQ=255)
                                           aTCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCt                                                                                                      <  6:272018/149‑1 (MQ=255)
                                               tACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGaa                                                                                                  <  4:34186/149‑1 (MQ=255)
                                                  ggCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATc                                                                                                                         <  3:309563/123‑1 (MQ=255)
                                                  ggCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATc                                                                                                                         >  4:309563/1‑123 (MQ=255)
                                                   gCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAg                                                                                              <  3:85124/149‑1 (MQ=255)
                                                      tGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAAc                                                                                           <  6:116175/149‑1 (MQ=255)
                                                      tGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAAc                                                                                           <  5:44844/149‑1 (MQ=255)
                                                        aTCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACg                                                                                          <  3:260106/148‑1 (MQ=255)
                                                          cTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACggg                                                                                        <  8:39381/148‑1 (MQ=255)
                                                            gCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCg                                                                                     >  8:183864/1‑149 (MQ=255)
                                                             cGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACggg                                                                                        <  6:139395/145‑1 (MQ=255)
                                                                ttctgtggCCATCTAGTAGTCATGCACATTCTACGGCGCTGTGTAGCTGTCGCCTCGGTTGGTCTACGCAGCCTTATGGTTCAGCAGGTCGCGGCTTTACACGCCAATCTGGTCGGATTGCTGGTCCTGGAAGTAGAACGGGTCGCTGa                                                                                 <  8:151341/142‑1 (MQ=255)
                                                                gaTGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGa                                                                                 >  1:216731/1‑149 (MQ=255)
                                                                gaTGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGa                                                                                 >  3:197913/1‑149 (MQ=255)
                                                                    gTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGt                                                                              <  7:183864/148‑1 (MQ=255)
                                                                     tCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGAtaga                                                                             >  7:151341/1‑145 (MQ=255)
                                                                        ccATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGgtg                                                                         <  2:149704/149‑1 (MQ=255)
                                                                         cATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGTAGGTCGCGGCTTTc                                                                                                                                   <  8:34178/90‑1 (MQ=255)
                                                                               gTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAg                                                                      <  5:368209/145‑1 (MQ=255)
                                                                                     aTGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGa                                                                                 <  3:285274/128‑1 (MQ=255)
                                                                                     aTGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGa                                                                                 >  4:285274/1‑128 (MQ=255)
                                                                                     aTGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCg                                                                            <  1:361218/133‑1 (MQ=255)
                                                                                     aTGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCg                                                                            >  2:361218/1‑133 (MQ=255)
                                                                                           aTTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGtcat                                                        >  5:240990/1‑147 (MQ=255)
                                                                                                  cgcgcgGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCAACAGGCGtt                                                >  2:95552/1‑148 (MQ=255)
                                                                                                   gcgcgGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAg                                                     <  4:256389/142‑1 (MQ=255)
                                                                                                   gcgcgGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAg                                                     >  3:256389/1‑142 (MQ=255)
                                                                                                       ggTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCa                                          <  6:240990/149‑1 (MQ=255)
                                                                                                        gtgtAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAg                                         <  5:22207/149‑1 (MQ=255)
                                                                                                             gCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCAc                                    <  4:38425/149‑1 (MQ=255)
                                                                                                                   gCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtgt                              >  8:11884/1‑149 (MQ=255)
                                                                                                                     cTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTTATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGtgc                            <  1:283570/149‑1 (MQ=255)
                                                                                                                     cTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGtgc                            <  4:197913/149‑1 (MQ=255)
                                                                                                                         gTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGTCACCCTTGtgctgc                         >  3:345261/1‑148 (MQ=255)
                                                                                                                             gTCAACGCACCTTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGtt                    <  1:294468/149‑1 (MQ=255)
                                                                                                                             gTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATGTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCAtgc                                           <  5:194911/126‑3 (MQ=255)
                                                                                                                             gTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATcc                                           >  6:194911/1‑126 (MQ=255)
                                                                                                                             gTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtgtg                             <  6:13937/140‑1 (MQ=255)
                                                                                                                             gTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtgtg                             >  5:13937/1‑140 (MQ=255)
                                                                                                                             gTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGtt                    >  6:109087/1‑149 (MQ=255)
                                                                                                                             gTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGtt                    >  8:127026/1‑149 (MQ=255)
                                                                                                                                aaCGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGTCACCCTTGTGCTGCCAGTTGa                  <  4:345261/148‑1 (MQ=255)
                                                                                                                                aaCGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGa                  <  7:127026/148‑1 (MQ=255)
                                                                                                                                 ccgcgaccTGCTGAACCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGAc                 >  7:110593/7‑148 (MQ=255)
                                                                                                                                 ccgcgaccTGCTGAACCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGAc                 >  7:105349/7‑148 (MQ=255)
                                                                                                                                 aCGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAg                                         <  7:133722/124‑1 (MQ=255)
                                                                                                                                 aCGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAg                                         >  8:133722/1‑124 (MQ=255)
                                                                                                                                     aCCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATcc            <  2:216731/149‑1 (MQ=255)
                                                                                                                                      accTGCTGAACCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCAACAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATc             >  7:46893/2‑147 (MQ=255)
                                                                                                                                       ccTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGGCCATCCAg          >  5:305784/1‑149 (MQ=255)
                                                                                                                                       ccTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAg          >  1:24379/1‑149 (MQ=255)
                                                                                                                                         tgctgGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAgcg        >  6:11364/1‑149 (MQ=255)
                                                                                                                                          gctgGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAgcgc       >  1:196425/1‑149 (MQ=255)
                                                                                                                                              gTTCAGCAGGTCGCGGCTTTCGACGCCAATCTGTTCTGATTCCAGGTCCTGTAAGTAGAACGTGTCGCTGATGTCGGTGTAGTCCACCt                                                               <  8:191853/89‑1 (MQ=255)
                                                                                                                                              gTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGAt   >  7:23562/1‑149 (MQ=255)
                                                                                                                                              gTTCAGCAGGTCGAGGCTTTCCACGACAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCt                                                               >  7:191853/1‑89 (MQ=255)
                                                                                                                                               ttCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGATc  <  1:95552/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGATC  >  NC_002947/490947‑491238

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