Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,990 A→T 14.0% I64F (ATC→TTC)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,9900AT14.0% 238.3 / 29.6 93I64F (ATC→TTC) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/45);  new base T (6/7);  total (41/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

TGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGA  >  NC_002947/3762841‑3763138
                                                                                                                                                      |                                                                                                                                                    
tGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGa                                                                                                                                                      <  1:1827954/150‑1 (MQ=255)
  gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTc                                                                                                                                                    >  1:640729/1‑150 (MQ=255)
      gCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATGTTCCCgg                                                                                                                                                <  1:1849061/150‑1 (MQ=255)
       cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGTCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATGTTCCCGGt                                                                                                                                               <  2:332527/150‑1 (MQ=255)
         tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCa                                                                                                                                             <  2:640729/150‑1 (MQ=38)
          cAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATGTTCCCGGTCAc                                                                                                                                            <  2:1227722/150‑1 (MQ=255)
           aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACg                                                                                                                                           >  1:613284/1‑150 (MQ=255)
            aCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACgg                                                                                                                                          >  1:2290689/1‑150 (MQ=255)
            aCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACgg                                                                                                                                          <  2:520670/150‑1 (MQ=255)
            aCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGAACATGATCGCGGTCACgg                                                                                                                                          <  1:2158188/150‑1 (MQ=255)
              gCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTc                                                                                                                                        <  1:978640/150‑1 (MQ=255)
                aGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGa                                                                                                                                      <  2:1336887/150‑1 (MQ=255)
                 ggCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAg                                                                                                                                     >  1:804817/1‑150 (MQ=255)
                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATcc                                                                                                                                                   >  1:2273488/1‑131 (MQ=255)
                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATcc                                                                                                                                                   <  2:2273488/131‑1 (MQ=255)
                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CCTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                                >  1:2075416/1‑150 (MQ=255)
                        aCGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGTCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTGGAGCGCCTgg                                                                                                                              >  1:1566213/1‑150 (MQ=255)
                         cGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGc                                                                                                                             <  1:2432483/150‑1 (MQ=255)
                           gCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGt                                                                                                                                               >  1:684560/1‑131 (MQ=39)
                           gCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGt                                                                                                                                               <  2:684560/131‑1 (MQ=38)
                            ccggtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGa                                                                                                                          <  2:1561467/147‑1 (MQ=255)
                              tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATc                                                                                                                        >  2:345669/1‑150 (MQ=255)
                                     cGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGt                                                                                                                 >  1:1402619/1‑150 (MQ=255)
                                      gaggagCGACGTGGCAATATCGAGGAAGAGTACGGCAAAGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                                                                                                                <  2:2230653/150‑1 (MQ=255)
                                       aggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgc                                                                                                               >  2:995817/1‑150 (MQ=255)
                                       aggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgc                                                                                                               <  1:2382675/150‑1 (MQ=255)
                                       aggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGAGCATGTCGttc                                                                                                               >  1:2434078/1‑148 (MQ=255)
                                            cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCaa                                                                                                          <  1:1798233/150‑1 (MQ=255)
                                                   gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATGTTCCCGGTCACgt                                                                                                                                          <  2:304105/111‑2 (MQ=255)
                                                   gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATGTTCCCGGTCACgt                                                                                                                                          >  1:304105/1‑110 (MQ=255)
                                                     aaTATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAcc                                                                                                 >  2:85567/1‑150 (MQ=255)
                                                        atCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCg                                                                                              >  1:1134180/1‑150 (MQ=255)
                                                         tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                                             >  2:1634477/1‑150 (MQ=255)
                                                           gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                                                    >  1:2124708/1‑93 (MQ=255)
                                                           gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                                                    <  2:2124708/93‑1 (MQ=255)
                                                               ccGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGATGCAACCACTTACCCCGCCTGCgg                                                                                       <  2:1623892/150‑1 (MQ=255)
                                                                cGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACt                                                                                                     >  2:929337/1‑135 (MQ=255)
                                                                cGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACt                                                                                                     <  1:929337/135‑1 (MQ=255)
                                                                cGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGc                                                                                      >  2:1535605/1‑150 (MQ=255)
                                                                    tACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACtt                                                                                  >  2:1265874/1‑150 (MQ=255)
                                                                      cGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                                <  2:494900/150‑1 (MQ=255)
                                                                         ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                             <  2:1986734/150‑1 (MQ=255)
                                                                               cgcCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTgctgct                                                                       <  1:85567/150‑1 (MQ=255)
                                                                                gcCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGctgctg                                                                      <  2:1865309/150‑1 (MQ=255)
                                                                                   gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGcc                                                                               <  2:1624606/138‑1 (MQ=255)
                                                                                   gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGcc                                                                               >  1:1624606/1‑138 (MQ=255)
                                                                                        tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                                              <  2:863343/150‑1 (MQ=255)
                                                                                        tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                                              <  2:1748854/150‑1 (MQ=255)
                                                                                        tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                                              <  2:2290689/150‑1 (MQ=255)
                                                                                        tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                                              <  2:1134180/150‑1 (MQ=255)
                                                                                          aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTAACGGTCGAGCGCCTGACGGAACCTGCCGTCCTCCCACACCGTACGCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGc                                                            >  2:744422/1‑150 (MQ=255)
                                                                                            gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCt                                                          >  2:1938356/1‑150 (MQ=255)
                                                                                             ccTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCtt                                                         <  1:1634477/150‑1 (MQ=255)
                                                                                              cTACAAGGCGGAAAAGAACAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                        <  1:1432821/150‑1 (MQ=255)
                                                                                               tACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTgg                                                       <  2:1948850/150‑1 (MQ=255)
                                                                                                   aGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcac                                                   >  1:1234993/1‑150 (MQ=255)
                                                                                                    ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGTAACATGATCCCGGTAAAGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                  >  1:53411/1‑150 (MQ=255)
                                                                                                    agcagAAAAGATCAAGCGGGCCTGGGCGGTGGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                  <  2:1751544/146‑1 (MQ=255)
                                                                                                      cGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGCAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgc                                                <  1:936131/150‑1 (MQ=255)
                                                                                                        gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcgc                                              >  1:889862/1‑150 (MQ=255)
                                                                                                        gAAAAGATAAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcgc                                              >  1:1627635/1‑150 (MQ=255)
                                                                                                         aaaaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                           >  1:1934260/1‑120 (MQ=255)
                                                                                                         aaaaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                           <  2:1934260/120‑1 (MQ=255)
                                                                                                         aaaaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                           >  1:116071/1‑120 (MQ=255)
                                                                                                         aaaaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                           <  2:116071/120‑1 (MQ=255)
                                                                                                             gATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                                        >  1:425772/1‑103 (MQ=255)
                                                                                                             gATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                                        <  2:425772/103‑1 (MQ=255)
                                                                                                             gATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGa                                         >  2:1003839/1‑150 (MQ=255)
                                                                                                                cAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                  >  2:792621/1‑138 (MQ=255)
                                                                                                                cAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                  <  1:792621/138‑1 (MQ=255)
                                                                                                                 aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                                     <  2:613284/150‑1 (MQ=255)
                                                                                                                 aaGCGCGCCTGCGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                                     <  1:1322502/150‑1 (MQ=255)
                                                                                                                  aGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGc                                    >  2:2204754/1‑150 (MQ=255)
                                                                                                                     gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                                 <  2:53411/150‑1 (MQ=255)
                                                                                                                       gcCTGGGCGTTCGA‑CGGGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                               <  1:399125/150‑1 (MQ=255)
                                                                                                                         cTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                             <  2:880399/150‑1 (MQ=255)
                                                                                                                           gggCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                        >  1:1141069/1‑121 (MQ=255)
                                                                                                                           gggCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                        <  2:1141069/121‑1 (MQ=255)
                                                                                                                            ggCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                                      >  1:859634/1‑138 (MQ=255)
                                                                                                                            ggCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                                      <  2:859634/138‑1 (MQ=255)
                                                                                                                               gTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGTCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCa                       <  1:781049/150‑1 (MQ=255)
                                                                                                                                   gaccGTGACCGGGATCATGTTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGa                    >  1:1977065/4‑150 (MQ=255)
                                                                                                                                   gA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt                   >  1:830214/1‑150 (MQ=255)
                                                                                                                                   gA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt                   <  1:995817/150‑1 (MQ=255)
                                                                                                                                        tGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCa               <  2:804817/150‑1 (MQ=255)
                                                                                                                                           ccGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAGAGAAGGATCGGGATCGCAccc            <  2:2434078/150‑1 (MQ=255)
                                                                                                                                            cGGGATCATGTTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCa           >  2:1001774/1‑150 (MQ=255)
                                                                                                                                            cGGGATCATGTTCCCGGTCA‑TGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAcc                                                                                                 <  1:911639/63‑1 (MQ=25)
                                                                                                                                            cGGGATCATGTTCCCGGTCA‑TGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAcc                                                                                                 >  2:911639/1‑63 (MQ=25)
                                                                                                                                            cGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCa           >  1:376104/1‑150 (MQ=255)
                                                                                                                                            cGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGACACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGGCAATCAGGATCGCACCCa           >  1:249167/1‑150 (MQ=255)
                                                                                                                                               gAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGt        <  2:1566213/150‑1 (MQ=255)
                                                                                                                                               gAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGt        <  2:889862/150‑1 (MQ=255)
                                                                                                                                                     gATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCCATCAGCATGGCACCCATGTCAAGGa  <  2:1402619/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGA  >  NC_002947/3762841‑3763138

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: ‑