Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,076,832 T→A 8.7% L69Q (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,8320TA 8.7% 171.3 / 7.8 58L69Q (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 = 6.65e-01

GTGCGCAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCT  >  NC_002947/1076684‑1076979
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                                 <  1:287146/126‑1 (MQ=255)
                           aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAg                                                                                                                                                 >  2:287146/1‑126 (MQ=255)
                           aTGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGcc                                                                                                                         >  1:1481899/1‑150 (MQ=255)
                                    tACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCt                                                                                                                                         >  2:1391103/1‑125 (MQ=255)
                                    tACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCt                                                                                                                                         <  1:1391103/125‑1 (MQ=255)
                                         gATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggcgg                                                                                                                           <  2:902228/134‑1 (MQ=255)
                                         gATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTggcgg                                                                                                                           >  1:902228/1‑134 (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                                                                                                                                                   <  1:1399587/86‑1 (MQ=255)
                                                                 cTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCtgctgctgc                                                                                                                                                   >  2:1399587/1‑86 (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                                                          <  2:625131/126‑1 (MQ=255)
                                                                                                                  gAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAg                                                          >  1:625131/1‑126 (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           >  2:863665/1‑150 (MQ=255)
                                                                                                                                         aGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCa           <  1:833084/150‑1 (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)
                                                                                                                                                  gctgcAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCt  <  2:1481899/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGCGCAGCCAGATGGGCCTGGTGCTGATGGAGTCGTACTGGATCACCGATTTGCCACAGCGTGCCTGGCAAGGTGTGGCCGGCCAGTTCGGTAGCCGCACCGAGCTGATTGCCGAAAACGAGAAGCTCAAGACCGAAGCCCTGCTGCTGCAGGGGCGCCTGCAGAAGCTGGCGGCCCTGACCGAGCAGAACGTGCGGCTGCGCGAGTTGCTCAACTCTTCGGCGCTGGTCAACGAAAAGGTTGAAGTGGCCGAGCTGATCGGTGTCGACCCCAACCCGTTTACCCACCGCATTCT  >  NC_002947/1076684‑1076979

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: