Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,616,606 G→C 6.6% A171A (GCG→GCC)  bepA → OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,616,6060GC 6.6% 200.2 / 5.6 61A171A (GCG→GCC) bepAOM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGAT  >  NC_000913/2616462‑2616747
                                                                                                                                                |                                                                                                                                             
ttCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                                                            <  3:254409/148‑1 (MQ=255)
   cGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                          >  3:286939/1‑147 (MQ=255)
     aTAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGt                                                                                                                                      <  7:30921/149‑1 (MQ=255)
      tAACGAAAGTCAACTGGCTTCACTAATGGCGCACGAAATCTCCCACGTCAACCAACGTCACCTGGCGCGAGCGATGTAAGATCAGCAGCGCAGCGAGCCGCTGACATGGGTCGGCGCGTTAGGAACTATTTTACTGGCGATGCCCAGTc                                                                                                                                     >  5:360534/1‑149 (MQ=255)
       aaCGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCa                                                                                                                                        <  1:99301/145‑1 (MQ=255)
       aaCGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCa                                                                                                                                        >  2:99301/1‑145 (MQ=255)
        aCGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTcca                                                                                                                                   >  7:7911/1‑148 (MQ=255)
         cGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGc                                                                                                                                          <  7:72198/141‑1 (MQ=255)
         cGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGc                                                                                                                                          >  8:72198/1‑141 (MQ=255)
         cGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGAGGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCg                                                                                                                                   >  5:310678/1‑148 (MQ=255)
          gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCa                                                                                                                                 <  8:378228/149‑1 (MQ=255)
                  aCTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGa                                                                                                                         <  8:9529/149‑1 (MQ=255)
                   cTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGAt                                                                                                                        <  6:331774/149‑1 (MQ=255)
                    tGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATg                                                                                                                       >  5:130650/1‑149 (MQ=255)
                        ttCAGTTATGGGGCACGTAATCTCACACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCTGCGGAGCGCGCCGCTGACCTGGGTCTGAGCGTTAGGTTCTATTTTACTGGCGATGGGCAGTCCGCAGGCGGGTATggcgg                                                                                                                   <  6:360534/149‑1 (MQ=255)
                        ttCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGGGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcgg                                                                                                                   >  2:186010/1‑149 (MQ=255)
                        ttCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcgg                                                                                                                   >  6:311213/1‑149 (MQ=255)
                         tCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcggc                                                                                                                  >  7:95567/1‑149 (MQ=255)
                            gTTATGGCGCACGAAATCTCCCCCGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGCCCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCt                                                                                                               >  8:203455/1‑149 (MQ=255)
                            gTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCt                                                                                                               >  6:42445/1‑149 (MQ=255)
                             ttATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCt                                                                                                               >  2:273411/1‑148 (MQ=255)
                                tGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGAc                                                                                                            <  1:16937/148‑1 (MQ=255)
                                    gcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcg                                                                                                                    <  8:203067/136‑1 (MQ=255)
                                    gcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcg                                                                                                                    >  7:203067/1‑136 (MQ=255)
                                       cGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTacac                                                                                                    <  4:238543/149‑1 (MQ=255)
                                           aTCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGc                                                                                                >  1:10769/1‑149 (MQ=255)
                                            tctcCCACGTCACCCAACGTCAACTGGCGCGATCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTATATTTTACTGGCTATGGCCAGTCCGCAGGCTGGGATGGCGGCGCTGACAGGTACACTGTCg                                                                                               >  8:374114/1‑149 (MQ=255)
                                                ccACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGaa                                                                                           >  5:365131/1‑149 (MQ=255)
                                                ccACGTCAACCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGaa                                                                                           >  2:4491/1‑149 (MQ=255)
                                                   cGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAcgc                                                                                        <  6:310678/149‑1 (MQ=255)
                                                       aCCCAACGTCACCTGGCGCGAGCGATGGAAGATCTGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCa                                                                                    <  6:130650/149‑1 (MQ=255)
                                                        cccAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAg                                                                                   <  8:95567/149‑1 (MQ=255)
                                                             cGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgat                                                                              <  1:186010/149‑1 (MQ=255)
                                                                 aCCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATc                                                                          <  5:42445/149‑1 (MQ=255)
                                                                          gagcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCgg                                                                                              >  4:135419/1‑120 (MQ=255)
                                                                          gagcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCgg                                                                                              <  3:135419/120‑1 (MQ=255)
                                                                          gagcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCAcc                                                                 >  6:82137/1‑149 (MQ=255)
                                                                           agcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCAccc                                                                <  1:273411/149‑1 (MQ=255)
                                                                            gcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCAcc                                                                 >  3:166555/1‑147 (MQ=255)
                                                                                 ggAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTc                                                                                     >  5:349558/1‑122 (MQ=255)
                                                                                 ggAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTc                                                                                     <  6:349558/122‑1 (MQ=255)
                                                                                  gAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCaaaa                                                         >  1:169309/1‑149 (MQ=255)
                                                                                   aaGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAAt                                                        >  6:106285/1‑149 (MQ=255)
                                                                                       tCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGaa                                                     >  8:55803/1‑148 (MQ=255)
                                                                                         agcagcGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACa                                                   <  4:286939/148‑1 (MQ=255)
                                                                                         agcagcGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAg                                                  >  2:196892/1‑149 (MQ=255)
                                                                                              cgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGc                                             <  5:311213/149‑1 (MQ=255)
                                                                                                   cgcgcCGCTGACCTGGGTCGGCGAGTTAGGTTCTATTTTACTGGCGATGGACAGTCCGCAGGCGGGGATGGCGGCGCTGACAGGTACACTGGCGGGAACGCGTCAGGTGATTATCATTTTCACCCAGAAAAATGAACAGGAAGCGGAcc                                        <  7:374114/149‑1 (MQ=255)
                                                                                                       ccgcTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAt                                    <  4:166555/149‑1 (MQ=255)
                                                                                                               cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTg                                  <  8:129506/143‑1 (MQ=255)
                                                                                                               cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTg                                  >  7:129506/1‑143 (MQ=255)
                                                                                                                           gTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACGGGAAGCGGACCGCATTGGTATTCAGGt                        <  5:259747/141‑1 (MQ=255)
                                                                                                                           gTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACGGGAAGCGGACCGCATTGGTATTCAGGt                        >  6:259747/1‑141 (MQ=255)
                                                                                                                                gTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTc            <  7:350900/148‑1 (MQ=255)
                                                                                                                                gTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCg           <  5:106285/149‑1 (MQ=255)
                                                                                                                                gTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCg           >  1:73964/1‑149 (MQ=255)
                                                                                                                                 ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTa                               >  2:285491/1‑128 (MQ=255)
                                                                                                                                 ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTa                               <  1:285491/128‑1 (MQ=255)
                                                                                                                                     tgcggaCTGGCCATCGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTc            >  2:89632/6‑143 (MQ=255)
                                                                                                                                     tgcggaCTGGCCATCGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTc            <  1:89632/138‑1 (MQ=255)
                                                                                                                                          aCTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGAt  <  1:4491/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGAT  >  NC_000913/2616462‑2616747

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