Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,031,922 T→A 5.7% L5M (TTG→ATG) ‡ PP_0889 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,031,9220TA 5.7% 324.9 / 7.3 105L5M (TTG→ATG) ‡PP_0889hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/52);  new base A (4/2);  total (51/54)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 108 total aligned reads displayed.

CAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGGCAG  >  NC_002947/1031774‑1032071
                                                                                                                                                    |                                                                                                                                                     
cAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCtt                                                                                                                                                      >  2:1661047/1‑150 (MQ=255)
cAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCtt                                                                                                                                                      <  1:772330/150‑1 (MQ=255)
     tctcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg                                                                                                                                                 <  1:260979/150‑1 (MQ=255)
       tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcgcg                                                                                                                                               >  2:122748/1‑150 (MQ=255)
         ggCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGaa                                                                                                                                             >  1:1688306/1‑150 (MQ=255)
         ggCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGaa                                                                                                                                             <  2:2223739/150‑1 (MQ=255)
           ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAgg                                                                                                                                           <  1:238805/150‑1 (MQ=255)
            cACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTAATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAggc                                                                                                                                          <  2:1280291/150‑2 (MQ=255)
             aCCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGc                                                                                                                                         >  2:1671256/1‑150 (MQ=255)
             aCCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgctcaca                                                                                                                                         >  1:2194395/1‑144 (MQ=255)
             aCCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgct                                                                                                                                             >  2:1476761/1‑144 (MQ=255)
             aCCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgct                                                                                                                                             <  1:1476761/146‑3 (MQ=255)
              ccGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                                        <  1:793632/150‑1 (MQ=255)
              ccGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                                        <  2:2216178/150‑1 (MQ=255)
                gCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGc                                                                                                                                      >  2:2192064/1‑150 (MQ=255)
                  aCGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                                    >  1:1443349/1‑150 (MQ=255)
                   cGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCaaa                                                                                                                                   <  1:1589770/150‑1 (MQ=255)
                    gggCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGCTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAg                                                                                                                                  <  2:1299824/150‑1 (MQ=255)
                         cTGGGCATGCACCGACGAAATCTTCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCNCCTNCAACTGTGAGCGCGNGCTTGCCCGCGAAGGGCTGCAAAGCAGTc                                                                                                                             <  1:122748/150‑1 (MQ=255)
                            ggCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCa                                                                                                                          >  2:2081691/1‑150 (MQ=255)
                            ggCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCa                                                                                                                          >  2:2062606/1‑150 (MQ=255)
                               aTGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                                       <  1:2171012/150‑1 (MQ=255)
                                  cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                                    >  1:443300/1‑150 (MQ=255)
                                  cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                                    <  2:828624/150‑1 (MQ=255)
                                  cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTCGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                                    >  1:176820/1‑150 (MQ=255)
                                   aCCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGACTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAgg                                                                                                                   >  2:1820701/1‑150 (MQ=255)
                                    ccGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGTTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGc                                                                                                                  <  2:981625/150‑1 (MQ=255)
                                    ccGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGc                                                                                                                  <  1:1715033/150‑1 (MQ=255)
                                     cGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCg                                                                                                                 >  2:1579947/1‑150 (MQ=255)
                                      gacggacgTCGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                                <  1:1284720/150‑1 (MQ=255)
                                           acgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACCGTAAGCGCGCACTTTTCCGCGTTCGGGGGCAAATCAGTCCCAATCAAGGCGAGcgcg                                                                                                           >  2:1204753/1‑150 (MQ=255)
                                             gTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGAt                                                                                                         >  1:195474/1‑150 (MQ=255)
                                                gCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGaa                                                                                                      >  2:1566078/1‑150 (MQ=255)
                                                 cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCGATGCGAGCTCGATGAAg                                                                                                     <  2:2435519/150‑1 (MQ=255)
                                                   gcgcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGt                                                                                                   >  2:109881/1‑150 (MQ=255)
                                                     gcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTAc                                                                                                 <  2:1524598/150‑1 (MQ=255)
                                                       aaGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTg                                                                                               >  2:632408/1‑150 (MQ=255)
                                                        aGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGt                                                                                              <  1:1147332/150‑1 (MQ=255)
                                                        aGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGt                                                                                              <  1:2316511/150‑1 (MQ=255)
                                                           tGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTgg                                                                                           <  1:2192064/150‑1 (MQ=255)
                                                           tGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCCCCGGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCGGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTgg                                                                                           <  1:191077/150‑1 (MQ=255)
                                                            gCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGt                                                                                          >  2:1278306/1‑150 (MQ=255)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                                 <  1:632408/150‑1 (MQ=255)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                                 >  2:320591/1‑150 (MQ=255)
                                                                        ggTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCt                                                                              >  2:1687704/1‑150 (MQ=255)
                                                                         gTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTg                                                                             >  2:1662858/1‑150 (MQ=255)
                                                                          ttCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTgg                                                                            <  2:443300/150‑1 (MQ=255)
                                                                             ggCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTg                                                                         <  1:1671256/150‑1 (MQ=255)
                                                                               cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTCCTGTTGGTCATCCTTGCGCTGGTTGCa                                                                       >  2:824456/1‑150 (MQ=255)
                                                                               cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                                       <  2:1477842/150‑1 (MQ=255)
                                                                                  tGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTgg                                                                                           >  2:230468/1‑127 (MQ=255)
                                                                                  tGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTgg                                                                                           <  1:230468/127‑1 (MQ=255)
                                                                                    aTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAg                                                                  <  1:366536/150‑1 (MQ=255)
                                                                                      tCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcg                                                                >  1:2049919/1‑150 (MQ=255)
                                                                                        cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGTAGTcc                                                                                                                            <  1:1376612/88‑1 (MQ=255)
                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGTATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCa                                                             >  1:1639911/1‑150 (MQ=255)
                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgctc                                                                                                                                            <  2:1644329/71‑4 (MQ=255)
                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgctc                                                                                                                                            >  1:1644329/1‑68 (MQ=255)
                                                                                           aTCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           <  1:1566078/150‑1 (MQ=255)
                                                                                           aTCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           <  1:1661047/150‑1 (MQ=255)
                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCt                                                        >  1:346844/1‑150 (MQ=255)
                                                                                                tgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGATAGCGCCAACGCGGa                                                      >  2:1320708/1‑150 (MQ=255)
                                                                                                   tgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           >  2:1168623/1‑142 (MQ=255)
                                                                                                   tgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           <  1:1168623/142‑1 (MQ=255)
                                                                                                      tGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGc                                                                                                                  >  2:1566225/1‑84 (MQ=255)
                                                                                                      tGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGc                                                                                                                  <  1:1566225/84‑1 (MQ=255)
                                                                                                      tGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTg                                                <  1:1062620/150‑1 (MQ=255)
                                                                                                        gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                              <  1:1250875/150‑1 (MQ=255)
                                                                                                        gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                              <  1:1687704/150‑1 (MQ=255)
                                                                                                         ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCtt                                                                                   <  1:2428257/112‑1 (MQ=255)
                                                                                                                  gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCa                                    <  2:195474/150‑1 (MQ=255)
                                                                                                                  gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCa                                    <  1:2062606/150‑1 (MQ=255)
                                                                                                                        aTCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCa                              <  2:1639911/150‑1 (MQ=255)
                                                                                                                         tCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCACTCCGGGCACCGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAg                             >  1:884505/1‑150 (MQ=255)
                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                                     <  2:233380/123‑1 (MQ=255)
                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGa                                          <  2:1294865/134‑1 (MQ=255)
                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGa                                          >  1:1294865/1‑134 (MQ=255)
                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                          >  1:113023/1‑150 (MQ=255)
                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                          >  1:2432951/1‑150 (MQ=255)
                                                                                                                            ccttcTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGATTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                          >  1:1664538/5‑150 (MQ=255)
                                                                                                                             cTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACg                         <  1:1278306/150‑1 (MQ=255)
                                                                                                                                cTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGGTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGAcc                      >  2:595442/1‑150 (MQ=255)
                                                                                                                                 tGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCa                                                                                                                                <  2:2080134/43‑1 (MQ=255)
                                                                                                                                 tGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCa                                                                                                                                >  1:2080134/1‑43 (MQ=255)
                                                                                                                                 tGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                               >  2:197066/1‑108 (MQ=255)
                                                                                                                                 tGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                               <  1:197066/108‑1 (MQ=255)
                                                                                                                                 tGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTg                                 >  2:123882/1‑138 (MQ=255)
                                                                                                                                 tGCAACTGTGAGCGCGGGCTTGCCAGCTAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTg                                 <  1:123882/138‑1 (MQ=255)
                                                                                                                                  gCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGTAGACTGGCCTGACCCGCAGTGGCAGAACGACCCt                    <  1:1662858/150‑1 (MQ=255)
                                                                                                                                  gCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCt                    >  2:281180/1‑150 (MQ=255)
                                                                                                                                  gCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCGCGATGAAGGTGCTGTTGGTCATCCTTGCCCTGGTTGCAGGTAGCGCCAAAGCTGAAACCTGGACATGTCGGCAGGGGCAGAACGACCCt                    >  2:180579/1‑150 (MQ=255)
                                                                                                                                         tgAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCAccc             >  1:922209/1‑150 (MQ=255)
                                                                                                                                         ccttcgcgGGCAAGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCAccc             >  2:351375/5‑150 (MQ=255)
                                                                                                                                           aGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTg                                                <  1:36999/113‑1 (MQ=255)
                                                                                                                                           aGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTg                                                >  2:36999/1‑113 (MQ=255)
                                                                                                                                           aGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTc           >  1:505040/1‑150 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGAc        <  2:2432951/150‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTg                   >  2:1722307/1‑137 (MQ=255)
                                                                                                                                                 ggCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCATTGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTgg     <  1:888594/150‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGGCAg  >  1:2272487/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGGCAG  >  NC_002947/1031774‑1032071

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: