Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,826:1 +C 6.5% coding (260/597 nt) PP_4825 → MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,8261.C 6.5% 249.0 / 11.1 62coding (260/597 nt)PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base . (29/29);  new base C (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

CCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGA  >  NC_002947/5489688‑5489970
                                                                                                                                           |                                                                                                                                                
ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                       <  1:1204679/150‑1 (MQ=255)
ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                       <  2:938043/150‑1 (MQ=255)
ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCGGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                       <  2:366349/150‑1 (MQ=255)
 cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                      >  1:702197/1‑150 (MQ=255)
                 gCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAg                                                                                                                      <  2:1516746/150‑1 (MQ=255)
                 gCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAg                                                                                                                      <  2:1400583/150‑1 (MQ=255)
                 gCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAg                                                                                                                      <  2:562811/150‑1 (MQ=255)
                  cAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCATGCTGTTCATTTCGGCCCTGGGCATGTCCCAGGGCAAGCCGGTGGTACAgg                                                                                                                     >  2:155766/1‑150 (MQ=255)
                  cAGCCTTGGCCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAgg                                                                                                                     <  2:1110184/150‑1 (MQ=255)
                     gggaggTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAg                                                                                                                  >  1:1291504/5‑150 (MQ=255)
                      tttGGTCGCCAGCGTGTTGCTTTATCTGTCCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCTCATCGGTGCTGGCAGCGTCCTGTTCATTTCGGCCCTGGGAACCGCCCAGGGCAAGCCGGCGGTACAGGCAGa                                                                                                                 <  1:44489/149‑1 (MQ=255)
                      cTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGa                                                                                                                 <  1:917767/150‑1 (MQ=255)
                             gCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGtt                                                                                                          >  2:869352/1‑150 (MQ=255)
                               cAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTca                                                                                                        >  2:991078/1‑150 (MQ=255)
                                   gtgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagcag                                                                                                    <  2:364989/150‑1 (MQ=255)
                                          tGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                             <  2:1647901/150‑1 (MQ=255)
                                           gTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                            >  1:574118/1‑118 (MQ=255)
                                           gTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                            <  2:574118/118‑1 (MQ=255)
                                            tACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGAcc                                                                                           >  1:710503/1‑150 (MQ=255)
                                               cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCaa                                                                                                                                   >  1:1300834/1‑107 (MQ=255)
                                               cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCaa                                                                                                                                   <  2:1300834/107‑1 (MQ=255)
                                                  ttCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAg                                                                                                                                        >  2:1590972/1‑99 (MQ=255)
                                                  ttCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAg                                                                                                                                        <  1:1590972/99‑1 (MQ=255)
                                                     gggcgggcGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                         <  1:454391/127‑1 (MQ=255)
                                                     gggcgggcGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                         >  2:454391/1‑127 (MQ=255)
                                                     gggcgggcGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTActg                                                                                  >  1:1222037/1‑148 (MQ=255)
                                                         gggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAgg                                                                                                                     >  1:903514/1‑111 (MQ=255)
                                                         gggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAgg                                                                                                                     <  2:903514/111‑1 (MQ=255)
                                                         gggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGa                                                                              <  2:1291504/150‑1 (MQ=255)
                                                           gcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCa                                                                                                                                         >  1:406208/1‑89 (MQ=255)
                                                           gcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCa                                                                                                                                         <  2:406208/89‑1 (MQ=255)
                                                             gATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCAt                                                                          >  1:287347/1‑150 (MQ=255)
                                                               tCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTc                                                                        <  1:1030947/150‑1 (MQ=255)
                                                                 ttCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCg                                                                      >  1:729121/1‑150 (MQ=255)
                                                                   cGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                    <  2:710503/150‑1 (MQ=255)
                                                                     gAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGa                                                                  >  1:540769/1‑150 (MQ=255)
                                                                      aaTTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGAc                                                                 >  1:142467/1‑150 (MQ=255)
                                                                        ttttCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagcag                                                                                                    >  1:71172/1‑113 (MQ=255)
                                                                        ttttCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagcag                                                                                                    <  2:71172/113‑1 (MQ=255)
                                                                               cATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGccc                                                        <  2:903612/150‑1 (MQ=255)
                                                                                  cACCGCCGAGGCCTTCCGCATCGGTGCTGGCAGCCTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCgg                                                     <  2:1079624/150‑1 (MQ=255)
                                                                                   aCCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGc                                                    >  1:50718/1‑150 (MQ=255)
                                                                                               ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGcc                                        >  2:668742/1‑150 (MQ=255)
                                                                                                     aTCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCGGGCACCATCGATGCCcagcgg                                  >  2:93405/1‑148 (MQ=255)
                                                                                                     aTCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgctg                                  >  1:1388100/1‑150 (MQ=255)
                                                                                                     aTCGGCGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACTATCGGTGCCctgctg                                  >  1:834020/1‑150 (MQ=255)
                                                                                                     aTCGGCGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgctg                                  <  1:971314/150‑1 (MQ=255)
                                                                                                       cGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGt                                >  1:1161850/1‑150 (MQ=255)
                                                                                                         gTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGa                              <  2:50718/150‑1 (MQ=255)
                                                                                                         gTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGa                              >  1:614428/1‑150 (MQ=255)
                                                                                                          tGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGAt                             <  2:1490817/150‑1 (MQ=255)
                                                                                                             tGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATggg                          <  1:93405/150‑1 (MQ=255)
                                                                                                                  gCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGa                     <  1:991078/150‑1 (MQ=255)
                                                                                                                         gTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagccagcc              >  2:1395848/1‑150 (MQ=255)
                                                                                                                          ttCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGACCCAGCCg             >  2:1308300/1‑150 (MQ=255)
                                                                                                                                 ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgc                                    <  2:193889/118‑1 (MQ=255)
                                                                                                                                 ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgc                                    >  1:193889/3‑120 (MQ=255)
                                                                                                                                 ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTg        >  1:740651/3‑148 (MQ=255)
                                                                                                                                 ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTg        <  2:740651/146‑1 (MQ=255)
                                                                                                                                 cGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTggg      <  2:1342215/150‑1 (MQ=255)
                                                                                                                                 cGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTggg      >  1:1342215/1‑150 (MQ=255)
                                                                                                                                     ccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCCGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGa  >  2:1058430/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                
CCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGA  >  NC_002947/5489688‑5489970

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