Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,609,479 C→T 5.5% T420M (ACG→ATG)  PP_4080 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,609,4790CT 5.5% 216.1 / 2.4 72T420M (ACG→ATG) PP_4080hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base T (2/2);  total (36/36)
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.01e-01

TCGAATCCTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGCAGTCGATGCTT  >  NC_002947/4609332‑4609628
                                                                                                                                                   |                                                                                                                                                     
tCGAATCCTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACgg                                                                                                                                                     >  1:1485182/1‑150 (MQ=255)
 cGAATCCTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGc                                                                                                                                                    <  2:1377163/150‑1 (MQ=255)
 cGAATCCTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGc                                                                                                                                                    <  1:1396789/150‑1 (MQ=255)
       cTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGATGGCCAAcgg                                                                                                                                              >  1:649683/1‑150 (MQ=255)
       cTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGATGGCCAAcgg                                                                                                                                              >  1:649294/1‑150 (MQ=255)
       cTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAAcgg                                                                                                                                              >  1:1242437/1‑150 (MQ=255)
       cTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAAcgg                                                                                                                                              >  2:436164/1‑150 (MQ=255)
        ttGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAAcggc                                                                                                                                             >  1:643341/1‑150 (MQ=255)
           aaTGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGATGGCCAAcggcttc                                                                                                                                          <  2:649683/150‑1 (MQ=255)
           aaTGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGATGGCCAAcggcttc                                                                                                                                          <  2:649294/150‑1 (MQ=255)
              gACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATTCGGGCTTCAGCTACTTCGAACTGGACCGCCCTGCCCCCGCCTGGGCCTGCTTGAAGTCGGCCACcgtcttcggc                                                                                                                                       <  2:1225463/150‑1 (MQ=255)
                             cAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGt                                                                                                                        <  1:436164/150‑1 (MQ=255)
                               aCTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGa                                                                                                                      <  2:705844/150‑1 (MQ=255)
                               aCTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGa                                                                                                                      >  2:1213380/1‑150 (MQ=255)
                                 cagcGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGt                                                                                                                    >  1:876056/4‑150 (MQ=255)
                                             cGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTg                                                                                                        >  2:1568141/1‑150 (MQ=255)
                                               ctgctgCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGaa                                                                                                      <  1:782986/150‑1 (MQ=255)
                                                    gcCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGCCCCCGCCTGGGCCTGCATGAAGACGGCCAACCGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCCGGGCTTGAACTGCa                                                                                                 >  2:99524/1‑150 (MQ=255)
                                                     ccccTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCtt                                                                                                         <  1:507353/141‑1 (MQ=255)
                                                     ccccTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCtt                                                                                                         >  2:507353/1‑141 (MQ=255)
                                                     ccccTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTAGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCaa                                                                                                >  2:529441/1‑149 (MQ=255)
                                                      cccTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGa                                                                                                       <  1:1194801/142‑1 (MQ=255)
                                                      cccTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGa                                                                                                       >  2:1194801/1‑142 (MQ=255)
                                                        cTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACggcttcggcttc                                                                                                                                    <  2:680810/111‑1 (MQ=255)
                                                        cTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACggcttcggcttc                                                                                                                                    >  1:680810/1‑111 (MQ=255)
                                                           ccGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTgg                                                                                          >  2:1346668/1‑150 (MQ=255)
                                                                   accacGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGc                                                                                                           <  2:6531/125‑1 (MQ=255)
                                                                   accacGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGc                                                                                                           >  1:6531/1‑125 (MQ=255)
                                                                    ccacGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCTAGGGCTTGAACTGCAGTTCTGGGCGATCAgg                                                                                 <  2:1009957/150‑1 (MQ=255)
                                                                    ccacGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAgg                                                                                 <  2:1242437/150‑1 (MQ=255)
                                                                     cacGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGCGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGCCTTAAACTGCAGTTCTGCGCGCTCAgcg                                                                                >  2:686620/1‑148 (MQ=21)
                                                                          gACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGaa                                                                           <  1:1346668/150‑1 (MQ=255)
                                                                              tGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGat                                                                       >  2:619041/1‑150 (MQ=255)
                                                                                cccTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACt                                                                                                    <  2:385620/119‑1 (MQ=255)
                                                                                cccTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACt                                                                                                    >  1:385620/1‑119 (MQ=255)
                                                                                cccTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAAt                                                                          >  1:729277/1‑145 (MQ=255)
                                                                                cccTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAAt                                                                          <  2:729277/145‑1 (MQ=255)
                                                                                 ccTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATg                                                                    <  1:619041/150‑1 (MQ=255)
                                                                                        tGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACa                                                             >  2:977206/1‑150 (MQ=255)
                                                                                        ggcgcgcCTCATCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACa                                                             <  2:1242786/149‑1 (MQ=255)
                                                                                             gCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGTCTTAAACTGCAGTGCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAac                                                        <  2:876056/150‑1 (MQ=255)
                                                                                              cTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATgg                                                                   <  1:1112941/138‑1 (MQ=255)
                                                                                              cTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATgg                                                                   >  2:1112941/1‑138 (MQ=255)
                                                                                                  aGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGtt                                                                                                                              >  2:249769/1‑75 (MQ=255)
                                                                                                  aGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGtt                                                                                                                              <  1:249769/75‑1 (MQ=255)
                                                                                                     tACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAaccctaccct                                                >  2:69380/1‑150 (MQ=255)
                                                                                                     tACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAaccctaccct                                                >  2:957114/1‑150 (MQ=255)
                                                                                                        ttCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCaa                                             <  1:957114/150‑1 (MQ=255)
                                                                                                           gAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACt                                                                                                    <  2:667487/92‑1 (MQ=255)
                                                                                                           gAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACt                                                                                                    >  1:667487/1‑92 (MQ=255)
                                                                                                             aCTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCa                                                                                                 >  2:1579425/1‑93 (MQ=255)
                                                                                                             aCTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCa                                                                                                 <  1:1579425/93‑1 (MQ=255)
                                                                                                                 gACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTg                                                                                           <  2:485585/95‑1 (MQ=255)
                                                                                                                 gACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTg                                                                                           >  1:485585/1‑95 (MQ=255)
                                                                                                                          gACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGc                           <  2:1485182/150‑1 (MQ=255)
                                                                                                                             cccGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAggg                                                                                >  2:1120006/1‑94 (MQ=255)
                                                                                                                             cccGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAggg                                                                                <  1:1120006/94‑1 (MQ=255)
                                                                                                                             cccGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGa                        >  1:332680/1‑150 (MQ=255)
                                                                                                                              ccGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTggg                                                                                         >  2:901658/1‑84 (MQ=255)
                                                                                                                              ccGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTggg                                                                                         <  1:901658/84‑1 (MQ=255)
                                                                                                                              ccGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGaa                       >  1:737461/1‑150 (MQ=255)
                                                                                                                                gCCTGGGCCGGCATGAAGACGGCCAACGGCTTCGGCTGCCACGTTGCAGGTGAGTTCCCAGGGCTTAAACTGCAGGTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCGTACCCTCAACGACGCGAAGCACGTTGCGGAAAc                     <  1:686620/150‑1 (MQ=255)
                                                                                                                                  cTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATg                                                                    <  2:1410801/101‑1 (MQ=255)
                                                                                                                                  cTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATg                                                                    >  1:1410801/1‑101 (MQ=255)
                                                                                                                                    gggCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGAt                 <  1:1568141/150‑1 (MQ=255)
                                                                                                                                     ggCCTGCCTGAAGACTGCCAACGGCTTCGGCTTCCACGTTTCAGGTCAATACGCACGGGATCGTCAACGGTGCTTGAACATAAGGGGTGAATGATATGGAAAACAGTAACCCTACCCGCAACGACGCGAAGCACGTTGCGGAAACGGAtt                <  1:529441/150‑1 (MQ=255)
                                                                                                                                         tGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGCa            >  2:681585/1‑150 (MQ=255)
                                                                                                                                          gCATGAAGACTGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGCAg           <  1:69380/150‑1 (MQ=255)
                                                                                                                                          gCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGCAg           >  1:9916/1‑150 (MQ=255)
                                                                                                                                           cATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGCAGt          <  1:977206/150‑1 (MQ=255)
                                                                                                                                                 gACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGc             >  1:1540679/1‑141 (MQ=255)
                                                                                                                                                 gACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGc             <  2:1540679/141‑1 (MQ=255)
                                                                                                                                                   cGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGCAGTCGATGCtt  <  1:1213380/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCGAATCCTTGAATGACCTGGTGCCTTTGCAACTGCCGGGCATCCCGCTGCTGCCCCTGCCGGTTGCACCACGCGACCTGCCCTTCCATGCGGGCTTCAGCTACTTCGAACTGGACCGCCGTGACCCCGCCTGGGCCTGCATGAAGACGGCCAACGGCTTCGGCTTCCACGTTGCAGGTGAGTTCCCAGGGCTTGAACTGCAGTTCTGGGCGATCAGGGGTGAATGATATGGAAAACAGTAACCCTACCCTCAACGACGCGAAGCACGTTGCGGAAACGGATTTGCAGTCGATGCTT  >  NC_002947/4609332‑4609628

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: