Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,076,829 T→A 8.6% L68Q (CTG→CAG)  mreC → cell shape‑determining protein MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,076,8290TA 8.6% 170.0 / 6.0 58L68Q (CTG→CAG) mreCcell shape‑determining protein MreC
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base A (3/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CCGTGCGCAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATT  >  NC_002947/1076682‑1076977
                                                                                                                                                   |                                                                                                                                                    
ccGTGCGCAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCctgctgc                                                                                                                                                    >  1:1215764/1‑150 (MQ=255)
ccGTGCGCAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCctgctgc                                                                                                                                                    >  1:512811/1‑150 (MQ=255)
  gTGCGCAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCctgctgctg                                                                                                                                                  >  2:1539943/1‑150 (MQ=255)
     cgcAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCTCCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAg                                                                                                                                               <  1:705028/150‑1 (MQ=255)
         gCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGgc                                                                                                                                           <  1:328837/150‑1 (MQ=255)
           cAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCAGCAGCAgggcttc                                                                                                                                         >  1:917487/1‑146 (MQ=255)
            aGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCAGCAGCAgggcttcg                                                                                                                                        <  2:917487/150‑6 (MQ=255)
              aTGGGCCTGGTGCTGAGGGAGTGGTACTGGATCACCCATTTGCCACAGCGTTCCTGTCAAGGTGTGTCCGGTCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTg                                                                                                                                      <  2:512811/150‑1 (MQ=255)
                gggCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCa                                                                                                                                    <  2:454414/150‑1 (MQ=255)
                   ccTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGaa                                                                                                                                 <  2:270100/150‑1 (MQ=255)
                            gATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggcggc                                                                                                                        >  1:599723/1‑150 (MQ=255)
                             aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAg                                                                                                                                               >  2:287146/1‑126 (MQ=255)
                             aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAg                                                                                                                                               <  1:287146/126‑1 (MQ=255)
                             aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGcc                                                                                                                       >  1:1481899/1‑150 (MQ=255)
                                      tACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCt                                                                                                                                       >  2:1391103/1‑125 (MQ=255)
                                      tACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCt                                                                                                                                       <  1:1391103/125‑1 (MQ=255)
                                           gATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggcgg                                                                                                                         >  1:902228/1‑134 (MQ=255)
                                           gATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggcgg                                                                                                                         <  2:902228/134‑1 (MQ=255)
                                           gATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGaa                                                                                                         >  2:784821/1‑150 (MQ=255)
                                               gacGATTTGCAACCACGTGCCTTGCAAGGTGTAGCCGACGGGGCCGGTGGCCGCACCTAGCTGATTGCCGAAAACGAGAAGCACAAGAACGAAACCCTGCGGCGGCAGGGGGGCCTGCGGATGCTGGCCGCCCTGTCCGAGCAGAACGTg                                                                                                     >  2:1334343/3‑150 (MQ=255)
                                                 cGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCg                                                                                                   >  1:396700/1‑150 (MQ=255)
                                                   aTTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGc                                                                                                 >  1:61988/1‑150 (MQ=255)
                                                    tttGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAAc                                                                                                        <  1:1632428/142‑1 (MQ=255)
                                                    tttGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAAc                                                                                                        >  2:1632428/1‑142 (MQ=255)
                                                      tGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCAACGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTgc                                                                                              <  1:267755/150‑1 (MQ=255)
                                                         cacaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTgg                                                                                                                            >  1:744890/1‑117 (MQ=255)
                                                         cacaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTgg                                                                                                                            <  2:744890/117‑1 (MQ=255)
                                                          acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTTCTGCAGGCGCCCCTGCagcagc                                                                                                                               >  1:1245708/1‑113 (MQ=255)
                                                          acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTTCTGCAGGCGCCCCTGCagcagc                                                                                                                               <  2:1245708/113‑1 (MQ=255)
                                                          acaGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGa                                                                                          >  2:1167759/1‑150 (MQ=255)
                                                                   cTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCtgctgctgc                                                                                                                                                 >  2:1399587/1‑86 (MQ=255)
                                                                   cTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCtgctgctgc                                                                                                                                                 <  1:1399587/86‑1 (MQ=255)
                                                                    tGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAAc                                                                                >  2:609073/1‑150 (MQ=255)
                                                                         aGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGTACTTGCGGCTGCGCGAGTTGCTCAACTCTTc                                                                           <  2:599723/150‑1 (MQ=255)
                                                                           gtgtGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCgg                                                                         >  2:1156895/1‑150 (MQ=255)
                                                                                gccggccAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGATCACCAGCGAAGCCCTGCTGCTGCATTGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTg                                                                    >  1:200418/1‑150 (MQ=255)
                                                                                                    cACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCATCGGCGCTGGTCAACGAAAAGGTTGAAGt                                                >  2:383203/1‑150 (MQ=255)
                                                                                                     aCCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTg                                               <  1:770033/150‑1 (MQ=255)
                                                                                                     aCCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTg                                               >  2:1008662/1‑150 (MQ=255)
                                                                                                            tGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGc                                        <  1:1539943/150‑1 (MQ=255)
                                                                                                                  ccGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTTGCGGCACTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCg                                  >  2:330731/1‑150 (MQ=255)
                                                                                                                    gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAAc                                                              >  1:1589369/1‑120 (MQ=255)
                                                                                                                    gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAAc                                                              <  2:1589369/120‑1 (MQ=255)
                                                                                                                    gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAg                                                        >  1:625131/1‑126 (MQ=255)
                                                                                                                    gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAg                                                        <  2:625131/126‑1 (MQ=255)
                                                                                                                    gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGgt                                >  1:722689/1‑150 (MQ=255)
                                                                                                                    gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGgt                                >  2:833084/1‑150 (MQ=255)
                                                                                                                          gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGaa                          >  2:115981/1‑149 (MQ=255)
                                                                                                                          gagaAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGAc                          >  1:672989/1‑150 (MQ=255)
                                                                                                                                cTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGc                                                                                     <  1:107158/85‑1 (MQ=255)
                                                                                                                                cTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGc                                                                                     >  2:107158/1‑85 (MQ=255)
                                                                                                                                      aCCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACCCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGttt              <  2:192375/150‑1 (MQ=255)
                                                                                                                                       ccGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTa             <  2:1215764/150‑1 (MQ=255)
                                                                                                                                        cGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTAc            >  1:1200169/1‑150 (MQ=255)
                                                                                                                                         gcgccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAGACCGCTGAAATGGTCGTACTGATCTGGGTCGACCCCAACCCGTTTaac           >  2:15490/5‑148 (MQ=255)
                                                                                                                                           gccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGc                                             <  2:201891/112‑1 (MQ=255)
                                                                                                                                           gccccTGCAGCAGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGc                                             >  1:201891/3‑114 (MQ=255)
                                                                                                                                           aGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCa         <  1:833084/150‑1 (MQ=255)
                                                                                                                                           aGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCa         >  2:863665/1‑150 (MQ=255)
                                                                                                                                             cccTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCAc        >  2:128604/1‑149 (MQ=255)
                                                                                                                                             cccTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCAc        <  1:128604/149‑1 (MQ=255)
                                                                                                                                              ccTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAAcc                   >  2:1155405/1‑137 (MQ=255)
                                                                                                                                              ccTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCg      >  2:74446/1‑150 (MQ=255)
                                                                                                                                              actgctgctgCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAAcc                   <  1:1155405/136‑1 (MQ=255)
                                                                                                                                               ctgctgctgCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGc     >  2:171222/1‑150 (MQ=255)
                                                                                                                                               ctgctgctgCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGc     >  1:288568/1‑150 (MQ=255)
                                                                                                                                                  ctgctgcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCAtt  >  2:98304/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGTGCGCAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATT  >  NC_002947/1076682‑1076977

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