Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,973,308:1 +T 6.6% coding (609/1713 nt) PP_2600 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,973,3081.T 6.6% 235.9 / 10.8 61coding (609/1713 nt)PP_2600hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (24/33);  new base T (2/2);  total (26/35)
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.78e-01

GACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAA  >  NC_002947/2973161‑2973457
                                                                                                                                                    |                                                                                                                                                     
gACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCg                                                                                                                                                     <  1:313759/150‑1 (MQ=255)
 aCCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCgg                                                                                                                                                    >  1:987817/1‑150 (MQ=255)
    cACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGcc                                                                                                                                                 >  2:145344/1‑150 (MQ=255)
        ccccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCAACATCTTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAggt                                                                                                                                             >  1:918449/1‑150 (MQ=255)
            gcgcACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTa                                                                                                                                         >  1:787683/1‑150 (MQ=255)
             cgcACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGAGCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTAc                                                                                                                                        >  2:461097/1‑150 (MQ=255)
                       cACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTggcggc                                                                                                                              <  1:1139315/150‑1 (MQ=255)
                           cTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                          >  1:116355/1‑150 (MQ=255)
                            tCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGcc                                                                                                                             <  1:716938/146‑1 (MQ=255)
                            tCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGcc                                                                                                                             >  2:716938/1‑146 (MQ=255)
                                   gcagGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATCTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGc                                                                                                                  >  2:1099256/1‑150 (MQ=255)
                                         gTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCCCTCCGTTGACCAGGTCCAGGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGt                                                                                                            <  1:749989/150‑1 (MQ=255)
                                         gTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGt                                                                                                            <  2:987817/150‑1 (MQ=255)
                                          tttGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTg                                                                                                           <  2:417166/150‑1 (MQ=255)
                                              gCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTgcagc                                                                                                       <  1:145344/150‑1 (MQ=255)
                                               ccGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTgcagca                                                                                                      <  1:962972/150‑1 (MQ=255)
                                               ccGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTgcagca                                                                                                      >  1:724380/1‑150 (MQ=255)
                                                 ggCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAgg                                                                                                    <  2:724380/150‑1 (MQ=255)
                                                    gcCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCt                                                                                                 >  1:1069401/1‑150 (MQ=255)
                                                     ccGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCtt                                                                                                <  1:49952/150‑1 (MQ=255)
                                                       gACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTg                                                                                              <  1:1012479/150‑1 (MQ=255)
                                                        aCGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTgg                                                                                             <  2:856692/150‑1 (MQ=255)
                                                        aCGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTgg                                                                                             <  1:406865/150‑1 (MQ=255)
                                                              gCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGc                                                                                       <  2:539524/150‑1 (MQ=255)
                                                              gCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGc                                                                                       <  2:311189/150‑1 (MQ=255)
                                                              gCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGc                                                                                       <  2:116355/150‑1 (MQ=255)
                                                                aGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCgg                                                                                     <  2:950352/150‑1 (MQ=255)
                                                                  gTGCACACCCAGGTCCTGTGCCGCTTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGt                                                                                   <  2:151925/150‑1 (MQ=255)
                                                                  gTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGt                                                                                   <  2:787683/150‑1 (MQ=255)
                                                                          ccAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGa                                                                           <  1:847014/150‑1 (MQ=255)
                                                                              gTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGaa                                                                       <  1:1198287/150‑1 (MQ=255)
                                                                               tCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTg                                                                                                           <  1:799108/113‑1 (MQ=255)
                                                                               tCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTg                                                                                                           >  2:799108/1‑113 (MQ=255)
                                                                                ccTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAgg                                                                     <  1:964567/150‑1 (MQ=255)
                                                                                ccTGTGCCGATCTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAgg                                                                     <  1:1099256/150‑1 (MQ=255)
                                                                                    tgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcggc                                                                 >  1:406720/1‑150 (MQ=255)
                                                                                    tgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcggc                                                                 >  1:406701/1‑150 (MQ=255)
                                                                                      ccGATTTAGCCAGGCGCGCGACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTg                                                                                                           <  2:1110468/106‑1 (MQ=255)
                                                                                      ccGATTTAGCCAGGCGCGCGACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTg                                                                                                           >  1:1110468/1‑106 (MQ=255)
                                                                                              gCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCa                                                       >  1:254601/1‑150 (MQ=255)
                                                                                                       gcCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAATCGCTTGCacac                                                                                                              >  1:338238/1‑86 (MQ=255)
                                                                                                       gcCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAATCGCTTGCacac                                                                                                              <  2:338238/86‑1 (MQ=255)
                                                                                                               tGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTGCGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCAt                                      <  2:1069401/150‑1 (MQ=255)
                                                                                                                  cacaCCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGc                                    >  2:595136/1‑149 (MQ=255)
                                                                                                                  cacaCCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGc                                    <  1:595136/149‑1 (MQ=255)
                                                                                                                    cacCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCg                                          <  2:1186106/141‑1 (MQ=255)
                                                                                                                    cacCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCg                                          >  1:1186106/1‑141 (MQ=255)
                                                                                                                       cGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTggcg                                                                                                                                <  1:557752/52‑1 (MQ=255)
                                                                                                                       cGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTggcg                                                                                                                                >  2:557752/1‑52 (MQ=255)
                                                                                                                       cGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGa                              >  2:448039/1‑150 (MQ=255)
                                                                                                                       cGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGAACAGGCGGCCCCCAACGCCCTTGAACACGTCGAGCGGGGCTTTGGTGGAGTGGGTCTTGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCAGGa                              >  2:136848/1‑150 (MQ=255)
                                                                                                                                   gCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcgc                  <  2:584041/150‑1 (MQ=255)
                                                                                                                                   gCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcgc                  >  1:584041/1‑150 (MQ=255)
                                                                                                                                   acctGGCCCGCCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcgc                  >  1:669634/4‑150 (MQ=255)
                                                                                                                                       ggCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCAt              <  1:939656/150‑1 (MQ=255)
                                                                                                                                       ggCCCGCCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCAt              <  2:669634/150‑1 (MQ=255)
                                                                                                                                           cGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTc          >  1:237751/1‑150 (MQ=255)
                                                                                                                                             gccATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGcc                                   <  2:775492/123‑1 (MQ=255)
                                                                                                                                             gccATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGcc                                   >  1:775492/2‑124 (MQ=255)
                                                                                                                                               cATGGCGGGCCAGGTGGTAAAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGTGTCATGTTGAGGAAGGCGGCAGGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTgc      <  1:461097/150‑1 (MQ=255)
                                                                                                                                                 tGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTgcgc    >  1:425393/1‑150 (MQ=255)
                                                                                                                                                   gCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCaa  <  2:237751/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAA  >  NC_002947/2973161‑2973457

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