Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 421,955 A→C 10.8% L666V (TTG→GTG)  PP_0348 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947421,9550AC10.8% 386.2 / 7.4 165L666V (TTG→GTG) PP_0348hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (77/66);  new base C (20/2);  total (97/68)
Fisher's exact test for biased strand distribution p-value = 8.23e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

CTGGAACCGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAA  >  NC_002947/421807‑422102
                                                                                                                                                    |                                                                                                                                                   
cTGGAACCGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAAc                                                                                                                                                    <  2:1091862/150‑1 (MQ=255)
 tGGAACCGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAAcg                                                                                                                                                   <  2:4546369/150‑1 (MQ=255)
  ggAACCGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCCAcgc                                                                                                                                                  >  2:3991746/1‑150 (MQ=255)
      ccGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCgg                                                                                                                                              >  1:4935822/1‑150 (MQ=255)
      ccGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCgg                                                                                                                                              <  2:1230174/150‑1 (MQ=255)
      ccGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCAGCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCgg                                                                                                                                              >  2:461266/1‑150 (MQ=255)
        gagaAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCTATGCCGTCCAGGATATCGATTTCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGtt                                                                                                                                            <  2:239275/150‑1 (MQ=255)
        gagaAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGtt                                                                                                                                            <  2:1394514/150‑1 (MQ=255)
        gagaAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGtt                                                                                                                                            >  2:4947515/1‑150 (MQ=255)
        gagaAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGtt                                                                                                                                            <  2:2112419/150‑1 (MQ=255)
         agaAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTg                                                                                                                                           <  2:779019/150‑1 (MQ=255)
         agaAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTg                                                                                                                                           <  1:4938215/150‑1 (MQ=255)
         agaAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTg                                                                                                                                           <  2:1647958/150‑1 (MQ=255)
           aaaaCTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTggg                                                                                                                                         <  1:3437182/150‑1 (MQ=255)
            aaaCTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGgc                                                                                                                                        >  1:4316999/1‑150 (MQ=255)
            aaaCTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTggtc                                                                                                                                        <  2:4707548/150‑3 (MQ=255)
               cTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCCACGCGCGGTTGGgcagc                                                                                                                                     >  1:1199505/1‑150 (MQ=255)
                       tGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCt                                                                                                                             >  1:4796802/1‑150 (MQ=255)
                         cATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCtga                                                                                                                           <  2:4328159/150‑1 (MQ=255)
                         cATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCtga                                                                                                                           <  2:3381980/150‑1 (MQ=255)
                         cATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCtga                                                                                                                           <  1:1558492/150‑1 (MQ=255)
                          aTCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCtgat                                                                                                                          <  1:1890765/150‑1 (MQ=255)
                              ttGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGaaa                                                                                                                      >  1:4565796/1‑150 (MQ=255)
                              ttGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGaaa                                                                                                                      <  2:1503215/150‑1 (MQ=255)
                               tGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAAt                                                                                                                     >  2:296529/1‑150 (MQ=255)
                               tGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAAt                                                                                                                     <  1:3535689/150‑1 (MQ=255)
                                 aGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATcc                                                                                                                   <  2:2735700/150‑1 (MQ=255)
                                      gCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCaa                                                                                                              >  2:1956654/1‑150 (MQ=255)
                                      gCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCaa                                                                                                              >  2:1956693/1‑150 (MQ=255)
                                        tAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAAc                                                                                                            >  2:4891250/1‑150 (MQ=255)
                                        tAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCa                                                                                                                  <  2:4505041/144‑1 (MQ=255)
                                        tAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAAc                                                                                                            >  1:3080719/1‑150 (MQ=255)
                                        tAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAAc                                                                                                            >  1:5428267/1‑150 (MQ=255)
                                        tAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAAc                                                                                                            >  1:5428151/1‑150 (MQ=255)
                                           gcgcCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCCACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGcc                                                                                                         >  2:559110/1‑150 (MQ=255)
                                            cgcCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCa                                                                                                        <  1:2174170/150‑1 (MQ=255)
                                              ccTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGc                                                                                                      >  1:1986793/1‑150 (MQ=255)
                                                tGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTc                                                                                                    >  1:3319603/1‑150 (MQ=255)
                                                  aTAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAg                                                                                                  >  2:1604793/1‑150 (MQ=255)
                                                     ggAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAgcgc                                                                                               >  1:2528563/1‑150 (MQ=255)
                                                      gAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGGAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGcc                                                                                              >  1:596911/1‑150 (MQ=255)
                                                      gAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGcc                                                                                              >  1:2428001/1‑150 (MQ=255)
                                                        aTGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCAt                                                                                            <  2:3278832/150‑1 (MQ=255)
                                                           ccGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGAc                                                                                         >  2:3713961/1‑150 (MQ=255)
                                                           ccGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGAc                                                                                         >  1:3905328/1‑150 (MQ=255)
                                                            cGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGAcc                                                                                        >  1:806748/1‑150 (MQ=255)
                                                            cGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGAcc                                                                                        >  2:4930352/1‑150 (MQ=255)
                                                                   ccccTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCtgat                                                                                                                          <  1:1115852/109‑1 (MQ=255)
                                                                    cccTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGc                                                                                <  2:4335648/150‑1 (MQ=255)
                                                                     ccTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCa                                                                               >  1:3774823/1‑150 (MQ=255)
                                                                      cTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCAt                                                                              >  2:1901539/1‑150 (MQ=255)
                                                                       ttCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCAtt                                                                             >  2:3573836/1‑150 (MQ=255)
                                                                       ttCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGAATGCCATGCGGCCAATCGACAGAACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGACAGCTCAGCGCCATCACCGACTCAGCAtt                                                                             >  2:3760267/1‑150 (MQ=255)
                                                                          gACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGTCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAACAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCAc                                                                          >  1:2031189/1‑150 (MQ=255)
                                                                          gACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCAc                                                                          >  1:2217446/1‑150 (MQ=255)
                                                                             aTCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCCCCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTc                                                                       >  2:1052113/1‑150 (MQ=255)
                                                                              tCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTcc                                                                      <  1:4983475/150‑1 (MQ=255)
                                                                                ctctTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCt                                                                    <  1:14299/150‑1 (MQ=255)
                                                                                 tctTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTg                                                                   >  1:1741813/1‑150 (MQ=255)
                                                                                 tctTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTg                                                                   >  1:2693821/1‑150 (MQ=255)
                                                                                  ctTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTgg                                                                  >  1:4331073/1‑150 (MQ=255)
                                                                                   ttGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGgt                                                                 >  1:5205953/1‑150 (MQ=255)
                                                                                   ttGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGgt                                                                 >  2:4015920/1‑150 (MQ=255)
                                                                                       ttGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGgtcgtc                                                             <  1:3083792/150‑1 (MQ=255)
                                                                                         gAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCAt                                                           >  1:326883/1‑150 (MQ=255)
                                                                                             gCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCCACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTg                                                       >  1:2331118/1‑150 (MQ=255)
                                                                                               cTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTgtc                                                     <  2:2505552/150‑1 (MQ=255)
                                                                                               cTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTgtc                                                     >  1:1348362/1‑150 (MQ=255)
                                                                                                tCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTgtcg                                                    >  1:5192237/1‑150 (MQ=255)
                                                                                                  gATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTgtcggc                                                  >  1:82287/1‑150 (MQ=255)
                                                                                                       cGTCCAGGATATCGATTGCCATGCGGCCAATCGACACCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTgtcggcgtcgg                                             <  2:1724117/150‑1 (MQ=255)
                                                                                                         tCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGcc                                           >  2:3982704/1‑150 (MQ=255)
                                                                                                          ccAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCt                                          <  1:972813/150‑1 (MQ=255)
                                                                                                           cAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTg                                         <  2:4383893/150‑1 (MQ=255)
                                                                                                           cAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTg                                         >  2:1681227/1‑150 (MQ=255)
                                                                                                             ggATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTggg                                       <  1:1598038/150‑1 (MQ=255)
                                                                                                                  tCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATc                                  >  1:2945838/1‑150 (MQ=255)
                                                                                                                   cGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCa                                 <  2:3319603/150‑1 (MQ=255)
                                                                                                                     aTTGCCATGCGGCAAATCGACAGCACAAGCAACGCGCGGTTGGGCAGCTGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAgc                               <  1:1025630/150‑1 (MQ=255)
                                                                                                                          cATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCg                          >  1:3248001/1‑150 (MQ=255)
                                                                                                                          cATGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCg                          >  1:3377194/1‑150 (MQ=255)
                                                                                                                           aTGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTggg                                       <  2:585046/136‑1 (MQ=255)
                                                                                                                            tGCGGCCAATCGACAGCACCAGCACCGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCggg                        >  2:5033532/1‑150 (MQ=255)
                                                                                                                               ggCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCt                     >  2:2702065/1‑150 (MQ=255)
                                                                                                                                   aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGc                 >  2:2232324/1‑150 (MQ=255)
                                                                                                                                    aTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGcc                >  2:4622464/1‑150 (MQ=255)
                                                                                                                                     tCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCa               <  2:2531455/150‑1 (MQ=255)
                                                                                                                                     tCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCa               >  2:1257413/1‑150 (MQ=255)
                                                                                                                                     tCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCa               <  2:2051101/150‑1 (MQ=255)
                                                                                                                                     tCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCa               >  2:1257730/1‑150 (MQ=255)
                                                                                                                                     tCGACAGCACCAGCAACGCGCGGGTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACGTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCa               <  2:2883083/150‑1 (MQ=255)
                                                                                                                                         cAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgta           <  1:4502092/150‑1 (MQ=255)
                                                                                                                                            cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgtagta        <  1:4511649/150‑1 (MQ=255)
                                                                                                                                             aCCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGtagtag       >  1:4097018/1‑150 (MQ=255)
                                                                                                                                               cAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTc     <  1:1829706/150‑1 (MQ=255)
                                                                                                                                                aGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCg    >  1:2878150/1‑150 (MQ=255)
                                                                                                                                                 gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCga   >  2:448580/1‑150 (MQ=255)
                                                                                                                                                 gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCga   >  2:4470957/1‑150 (MQ=255)
                                                                                                                                                  cAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaa  <  2:91563/150‑1 (MQ=255)
                                                                                                                                                  cAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaa  <  2:3016840/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGGAACCGAGAAAACTGGTACCTGCATCGTTGAGGTGGCTAAGCGCCTGATAGGAATGCCGCATCACCCCTTCGACATCCTCTTGGTTGAAGGCCTCGATGCCGTCCAGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAA  >  NC_002947/421807‑422102

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