Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 115,951 T→G 6.7% L85W (TTG→TGG)  PP_0110 → CyoE‑like protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947115,9510TG 6.7% 316.2 / 6.7 107L85W (TTG→TGG) PP_0110CyoE‑like protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base T (52/47);  new base G (5/3);  total (57/50)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

TTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACACC  >  NC_002947/115803‑116099
                                                                                                                                                    |                                                                                                                                                    
ttCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGTTGTTTTTCGGCAACCTGGGCATCGGCCTGTGGGGTGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTg                                                                                                                                                     <  1:418606/150‑1 (MQ=255)
ttCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTg                                                                                                                                                     >  1:1722835/1‑150 (MQ=255)
 tCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTgg                                                                                                                                                    <  1:604287/150‑1 (MQ=255)
  ccTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGc                                                                                                                                                   <  1:2276989/150‑1 (MQ=255)
   cTGGCTATCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGcc                                                                                                                                                  >  1:1418910/1‑150 (MQ=255)
   cTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGcc                                                                                                                                                  >  1:1608894/1‑150 (MQ=255)
   cTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGcc                                                                                                                                                  >  2:1496006/1‑150 (MQ=255)
       ccACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                                                              <  1:778647/150‑1 (MQ=255)
        cACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCACCCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGc                                                                                                                                                   <  1:1351658/144‑1 (MQ=255)
        cACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGc                                                                                                                                                   >  2:1351658/1‑144 (MQ=255)
           ccGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGACTGTGTGCGGGTGGTGCGGCGGTGGGCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCg                                                                                                                                          >  1:1497863/1‑150 (MQ=255)
            cgggcgggCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTGGGCCAAAGGCCgt                                                                                                                                         <  1:1080481/150‑2 (MQ=255)
              ggcgggcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGt                                                                                                                                       <  1:1080731/150‑1 (MQ=255)
                   gcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAg                                                                                                                                   >  1:98635/1‑149 (MQ=255)
                     gTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGccc                                                                                                                                >  2:1319405/1‑150 (MQ=255)
                          cTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGGCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCGCCAACGGCCGTGGGACCACGGCGGAGTAGAGCCCCTGcc                                                                                                                           >  2:1033815/1‑150 (MQ=255)
                                 gTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtg                                                                                                                    <  2:1165983/150‑1 (MQ=255)
                                  tGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTAAATCATGTGGTGGCTCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGACTTTTTCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgt                                                                                                                   <  1:1397376/150‑1 (MQ=255)
                                   gCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgtt                                                                                                                  <  1:1055798/150‑1 (MQ=255)
                                     tgttgtTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGc                                                                                                                            <  2:1427142/138‑1 (MQ=255)
                                     tgttgtTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGc                                                                                                                            >  1:1427142/1‑138 (MQ=255)
                                     tgttgtTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgttgt                                                                                                                <  2:1515534/150‑1 (MQ=255)
                                     tgttgtTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgttgt                                                                                                                >  1:289163/1‑150 (MQ=255)
                                         gttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGc                                                                                                            >  2:548674/1‑150 (MQ=255)
                                          ttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAg                                                                                                                                   <  2:663580/126‑1 (MQ=255)
                                           tCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCt                                                                                                          <  2:1497863/150‑1 (MQ=255)
                                             ggCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGt                                                                                                                                       <  1:1574982/119‑1 (MQ=255)
                                             ggCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGt                                                                                                                                       >  2:1574982/1‑119 (MQ=255)
                                                 aCCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGc                                                                                                                                                   >  2:1108343/1‑103 (MQ=255)
                                                   cTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCCGCGAATCAGCCTGTGGCGATTGGTCGCGTTGATGCCCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTg                                                                                                  >  1:1023435/1‑150 (MQ=255)
                                                    tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGa                                                                                                                                         >  1:2227195/1‑110 (MQ=255)
                                                    tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGa                                                                                                                                         <  2:2227195/110‑1 (MQ=255)
                                                    tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTgg                                                                                                 <  1:1009396/150‑1 (MQ=255)
                                                         aTCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCgctg                                                                                            >  1:768987/1‑150 (MQ=255)
                                                              ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGc                                                                                                                                  >  2:1879321/1‑107 (MQ=255)
                                                              ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGc                                                                                                                                  <  1:1879321/107‑1 (MQ=255)
                                                              ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTgg                                                                                        <  2:438298/149‑1 (MQ=255)
                                                              ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTgg                                                                                        >  1:438298/1‑149 (MQ=255)
                                                                  tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                   <  2:768987/150‑1 (MQ=255)
                                                                  tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                   >  1:1286084/1‑150 (MQ=255)
                                                                  tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTGGGCCAAAGGcc                                                                                                                                           <  2:1979195/94‑1 (MQ=255)
                                                                  tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTGGGCCAAAGGcc                                                                                                                                           >  1:1979195/1‑94 (MQ=255)
                                                                   gtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCTAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATgg                                                                                  <  2:1400476/150‑1 (MQ=255)
                                                                    tgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGt                                                                                 >  1:2285541/1‑150 (MQ=255)
                                                                       gggTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCct                                                                              >  2:1498249/1‑150 (MQ=255)
                                                                         gtggtgCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                                                              <  2:2359064/84‑1 (MQ=255)
                                                                         gtggtgCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                                                              >  1:2359064/1‑84 (MQ=255)
                                                                              gCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                                                                       <  2:1880412/150‑1 (MQ=255)
                                                                              gCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                                                                       <  2:1608894/150‑1 (MQ=255)
                                                                                ggCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGgt                                                                     <  1:2164392/150‑1 (MQ=255)
                                                                                 gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                                                                       <  2:288024/147‑1 (MQ=255)
                                                                                 gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                                                                       >  1:288024/1‑147 (MQ=255)
                                                                                 gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGgtg                                                                    <  1:1033815/150‑1 (MQ=255)
                                                                                   aGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGTGTAGTGTCCTGCTGCTGGTGtt                                                                  >  2:1159051/1‑150 (MQ=255)
                                                                                       gTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCAcc                                                              >  1:1791086/1‑150 (MQ=255)
                                                                                        tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGgcgc                                                                                              >  1:1464660/1‑117 (MQ=255)
                                                                                          aaTCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAAc                                                           >  1:947792/1‑150 (MQ=255)
                                                                                            tCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACACACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGc                                                         >  2:1389191/1‑150 (MQ=255)
                                                                                             tctgtgGTGGATCGGCGCGTTGATGCGCTGATGGCAAGCACCCACAAACGGCCTTTGGCCCAAGGCAGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGcc                                                        <  2:1623438/148‑1 (MQ=255)
                                                                                               tgtgGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCt                                                      >  1:60324/1‑150 (MQ=255)
                                                                                                 tggtggATCGGCGCATTGATGCGCTGGTCGCACCCAAAAACTGACAACTTCGAGCTGAAGGCCGAGTCGAAGCCCTGCCCGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGGATGGACCTGCTGCTGGTGTTCACCAACGcccaga                                                    >  1:1623438/1‑147 (MQ=255)
                                                                                                  ggtggATCGGCGCATTGATGCGCTGATGGCAAGCACCCACAGACGGCCTTTGGCCCAAGGCCGAGTCCAGCACCTGCCGGCGCTGTTGTTCGCCTTGGCCATGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCAc                                                   <  1:520975/150‑1 (MQ=255)
                                                                                                    tggATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCg                                                 >  1:1840434/1‑150 (MQ=255)
                                                                                                     ggATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCActgc                                                >  2:1153957/1‑147 (MQ=255)
                                                                                                       aTCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCt                                              <  2:60324/150‑1 (MQ=255)
                                                                                                         cGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAAAGCCCTCACCGCCTgg                                            >  2:722552/1‑150 (MQ=255)
                                                                                                           gcgcATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCt                                          <  2:2252845/150‑1 (MQ=255)
                                                                                                           gcgcATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCt                                          >  1:2028209/1‑150 (MQ=255)
                                                                                                              cATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGGCCTCACCGCCTGGCTGAc                                       >  1:166439/1‑150 (MQ=255)
                                                                                                                ttcaTGCGCTGATGGCACGCACTCACAAACGGCCTTTGGCACAAGGGCGAGTCGAGCCCCTGCCGGCGCTGTTGTGCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTAGTGTTCACCAACGCCCTCACCGCCTGGCTGAccc                                     <  2:1052878/147‑1 (MQ=255)
                                                                                                                     gcgcTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGAGCTGTTGTTCGCCTTGGCCCTGGCgctg                                                                                            >  1:2341957/1‑90 (MQ=255)
                                                                                                                       gcTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGgt                                                                     <  1:1014397/111‑1 (MQ=255)
                                                                                                                       gcTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCCCCAACGCCCTCACCGCCTCGCTGACCCTGGCCTc                              >  2:409791/1‑150 (MQ=255)
                                                                                                                        cTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                                                              >  2:2210103/1‑69 (MQ=255)
                                                                                                                                acgcacCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCaa                                                            >  2:628629/1‑111 (MQ=255)
                                                                                                                                acgcacCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCaa                                                            <  1:628629/111‑1 (MQ=255)
                                                                                                                                acgcacCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCgg                     >  1:1768646/1‑150 (MQ=255)
                                                                                                                                acgcacCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCgg                     >  1:747550/1‑150 (MQ=255)
                                                                                                                                acgcaaCCTGAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCgg                     <  1:1308028/150‑1 (MQ=255)
                                                                                                                                    acCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTAt                 >  1:1661863/1‑150 (MQ=255)
                                                                                                                                       cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGGTGCATTTTGCCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg              <  2:1418910/150‑1 (MQ=255)
                                                                                                                                       cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg              >  2:844601/1‑150 (MQ=255)
                                                                                                                                       cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg              >  2:238020/1‑150 (MQ=255)
                                                                                                                                       cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg              >  2:716426/1‑150 (MQ=255)
                                                                                                                                       cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg              >  1:510576/1‑150 (MQ=255)
                                                                                                                                            tcGGCCTTGGGCCAAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                                                              >  2:1936790/2‑49 (MQ=25)
                                                                                                                                            tcGGCCTTGGGCCAAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                                                              <  1:1936790/48‑1 (MQ=25)
                                                                                                                                            tcGGCCTTGGGCCAAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGct         >  2:421986/2‑150 (MQ=255)
                                                                                                                                            aCGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                   >  2:1111356/1‑76 (MQ=255)
                                                                                                                                            aCGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                   <  1:1111356/76‑1 (MQ=255)
                                                                                                                                             cGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGctc        <  1:844601/150‑1 (MQ=255)
                                                                                                                                              ggCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGctct       >  2:2086320/1‑150 (MQ=255)
                                                                                                                                              ggCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGctc        >  1:601580/1‑149 (MQ=255)
                                                                                                                                                 cTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTaca    <  2:1286084/150‑1 (MQ=255)
                                                                                                                                                  tttGGCGCAATGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTacac   >  2:1299735/1‑150 (MQ=255)
                                                                                                                                                  tttGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                   >  2:1000975/1‑70 (MQ=255)
                                                                                                                                                  tttGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                   <  1:1000975/70‑1 (MQ=255)
                                                                                                                                                  tttGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTacac   >  2:12387/1‑150 (MQ=255)
                                                                                                                                                  tttGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTacac   >  1:2306947/1‑150 (MQ=255)
                                                                                                                                                   ttGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACAcc  <  2:289163/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACACC  >  NC_002947/115803‑116099

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