Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,117,220 C→T 5.8% P285S (CCA→TCA)  PP_1884 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,117,2200CT 5.8% 199.2 / 2.4 69P285S (CCA→TCA) PP_1884hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base T (2/2);  total (33/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 = 6.97e-01

GGCTGGTCGCCATGGGGCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGGCGCGGTAGTGTTGAT  >  NC_002947/2117071‑2117369
                                                                                                                                                     |                                                                                                                                                     
ggCTGGTCGCCATGGGGCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGc                                                                                                                                                       >  1:2068694/1‑150 (MQ=255)
     gTCGCCATGGGGCGCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTg                                                                                                                                                  >  2:2040898/1‑150 (MQ=255)
       ccccATGGCGACCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGa                                                                                                                                                <  1:809/148‑1 (MQ=255)
        gCCATGGGGCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAg                                                                                                                                               >  2:2006266/1‑150 (MQ=255)
          cATGGGGCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTAGATAAGGCCTGGGGGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGtt                                                                                                                                             >  2:1302764/1‑150 (MQ=255)
              gggCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTc                                                                                                                                         <  1:1716529/150‑1 (MQ=255)
                gCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTc                                                                                                                                         >  1:1813941/1‑148 (MQ=255)
                gCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTc                                                                                                                                         <  2:1813941/148‑1 (MQ=255)
                 cccAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACg                                                                                                                                      >  2:1085917/1‑150 (MQ=255)
                  ccAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGa                                                                                                                                     >  1:1683932/1‑150 (MQ=255)
                    aGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGGGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAAt                                                                                                                                   <  1:1302764/150‑1 (MQ=255)
                           ggCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGttgaat                                                                                                                            <  1:515198/150‑1 (MQ=255)
                             ccGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGttgaattt                                                                                                                          <  2:906930/150‑1 (MQ=255)
                                cAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGTCATTGGAGTTTTTCACGAATGttgaatttgaa                                                                                                                       >  2:2316492/1‑150 (MQ=255)
                                cAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGTCATTGGAGTTTTTCACGAATGttgaatttgaa                                                                                                                       >  2:2318417/1‑150 (MQ=255)
                                   aaTATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTgaatttgaattt                                                                                                                    <  2:1683932/150‑1 (MQ=255)
                                   aaTATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTgaatttgaattt                                                                                                                    <  2:2290158/150‑1 (MQ=255)
                                         ccTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGt                                                                                                              >  2:1389383/1‑150 (MQ=255)
                                           tGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGc                                                                                                            <  2:1559662/150‑1 (MQ=255)
                                           tGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGc                                                                                                            >  1:100945/1‑150 (MQ=255)
                                             gCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAg                                                                                                          <  1:1313003/150‑1 (MQ=255)
                                              cccGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGa                                                                                                         >  2:1589970/1‑150 (MQ=255)
                                                  gAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAAc                                                                                                     <  1:2322037/150‑1 (MQ=255)
                                                   aaTACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGTCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACg                                                                                                    <  1:2316492/150‑1 (MQ=255)
                                                   aaTACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGTCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACg                                                                                                    <  1:2318417/150‑1 (MQ=255)
                                                   aaTACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACg                                                                                                    <  1:1039262/150‑1 (MQ=255)
                                                    aTACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACgg                                                                                                   <  1:2040898/150‑1 (MQ=255)
                                                    aTACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACgg                                                                                                   <  1:1085917/150‑1 (MQ=255)
                                                     tACCATTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGt                                                                                                  >  2:1138148/1‑150 (MQ=255)
                                                     caCCGTTACTTGATTCCCGATGACCACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTGGAATTTGAATTTGAACGTGCAGAAAACGGt                                                                                                  <  1:908863/149‑1 (MQ=255)
                                                         gTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATc                                                                                              <  2:100945/150‑1 (MQ=255)
                                                                gATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTCGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTCGAATTTGAACGTGCAGAAAACGGTGATCTGACAGc                                                                                       >  2:1593950/1‑150 (MQ=255)
                                                                  ttCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGa                                                                                     >  2:1931135/1‑150 (MQ=255)
                                                                     ccGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGtt                                                                                                                                <  1:382614/104‑1 (MQ=255)
                                                                     ccGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGtt                                                                                                                                >  2:382614/1‑104 (MQ=255)
                                                                      cGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGc                                                                                 >  1:1955902/1‑150 (MQ=255)
                                                                              aCGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGTGTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAg                                                                         <  1:145001/150‑1 (MQ=255)
                                                                              aCGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAg                                                                         <  1:1931135/150‑1 (MQ=255)
                                                                                  tGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGt                                                                                                  <  1:1604046/121‑1 (MQ=255)
                                                                                  tGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGt                                                                                                  >  2:1604046/1‑121 (MQ=255)
                                                                                   gATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGtgtg                                                                    >  2:1384929/1‑150 (MQ=255)
                                                                                    aTTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACt                                                                                   >  1:107473/1‑134 (MQ=255)
                                                                                    aTTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACt                                                                                   <  2:107473/134‑1 (MQ=255)
                                                                                         tCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATcgcc                                                              >  1:615450/1‑150 (MQ=255)
                                                                                            gCGACGTTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCa                                                           <  1:1384929/150‑1 (MQ=255)
                                                                                                cGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAgaga                                                       <  2:1059642/150‑1 (MQ=255)
                                                                                                        cTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGt                                               >  1:1915712/1‑150 (MQ=255)
                                                                                                        cTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGt                                               >  2:9664/1‑150 (MQ=255)
                                                                                                         ttGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTAACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTACAGGTGTGATCGCCGCAGAGAAAGAAAGtt                                              >  1:370185/1‑150 (MQ=255)
                                                                                                             tAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGGGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGt                                          >  1:1334771/1‑150 (MQ=255)
                                                                                                              aaGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGtt                                         <  1:2186308/150‑1 (MQ=255)
                                                                                                                   cTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTcc                                    <  2:615450/150‑1 (MQ=255)
                                                                                                                       gTGGCGTTGGCTACCGTGTTAAGCAAACTGACACTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTAc                                >  1:311696/1‑150 (MQ=255)
                                                                                                                         ggCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTa                              <  2:2068694/150‑1 (MQ=255)
                                                                                                                         ggCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTAGGTTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGATAGAAAGTTTTGTTAGTCCTTACTa                              <  1:1593950/150‑1 (MQ=255)
                                                                                                                          gCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTAt                             <  1:1389383/150‑1 (MQ=255)
                                                                                                                              tGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTa                                 <  2:1469308/142‑1 (MQ=255)
                                                                                                                              tGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTa                                 >  1:1469308/1‑142 (MQ=255)
                                                                                                                                  tACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTCCTCCTTACTACGacgccta                     >  2:1957660/1‑144 (MQ=255)
                                                                                                                                     cGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTg                  >  1:1657000/1‑150 (MQ=255)
                                                                                                                                       tgtTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTg                                                                                                             >  2:1304305/1‑57 (MQ=255)
                                                                                                                                       tgtTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTg                                                                                                             <  1:1304305/57‑1 (MQ=255)
                                                                                                                                       tgtTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGgc                >  2:2059437/1‑150 (MQ=255)
                                                                                                                                             gCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGGCGCGGTa          >  1:2021249/1‑150 (MQ=255)
                                                                                                                                             gCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGGCGCGGTa          <  2:2021249/150‑1 (MQ=255)
                                                                                                                                                 aCTGCCATTGGAGTTTTTCACGAATGTTGAGTTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATcgcc                                                              >  2:1406192/1‑94 (MQ=255)
                                                                                                                                                 aCTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGCACGTTCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTcc                                    <  1:1343233/120‑1 (MQ=255)
                                                                                                                                                 aCTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATcgcc                                                              <  1:1406192/94‑1 (MQ=255)
                                                                                                                                                 aCTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTcc                                    >  2:1343233/1‑120 (MQ=255)
                                                                                                                                                 aCTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGGCGCGGTAgtgt      <  1:9664/150‑1 (MQ=255)
                                                                                                                                                     ccATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGGCGCGGTAGTGTTGAt  >  2:58539/1‑150 (MQ=255)
                                                                                                                                                     ccATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGGCGCGGTAGTGTTGAt  <  2:1915712/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCTGGTCGCCATGGGGCCCAGTAACGGGCCGCAGAATATGCCTGGCCCGGAATACCGTTACTTGATTCCCGATGACAACGGTGATTACTCGGCGACGGTACTGCTTGATAAGGCCTGGGTGGCGTTGGCTACCGTGTTAAGCAAACTGCCATTGGAGTTTTTCACGAATGTTGAATTTGAATTTGAACGTGCAGAAAACGGTGATCTGACAGCGACTGCCATTGCAGGTGTGATCGCCGCAGAGAAAGAAAGTTTTGTTAGTCCTTACTATGAGTCTTACTGGCGCGGTAGTGTTGAT  >  NC_002947/2117071‑2117369

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