Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 115,956 C→A 5.8% Q87K (CAA→AAA)  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,9560CA 5.8% 319.5 / 7.2 105Q87K (CAA→AAA) PP_0110CyoE‑like protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base C (53/44);  new base A (3/3);  total (57/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

Only 100 of 110 total aligned reads displayed.

CCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACACCGGCTT  >  NC_002947/115810‑116104
                                                                                                                                                  |                                                                                                                                                    
ccACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                                                                   <  1:778647/150‑1 (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)
                   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                                                                                                                                 >  1:1427142/1‑138 (MQ=255)
                              tgttgtTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGc                                                                                                                                 <  2:1427142/138‑1 (MQ=255)
                              tgttgtTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgttgt                                                                                                                     >  1:289163/1‑150 (MQ=255)
                              tgttgtTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgttgt                                                                                                                     <  2:1515534/150‑1 (MQ=255)
                                  gttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGc                                                                                                                 >  2:548674/1‑150 (MQ=255)
                                   ttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAg                                                                                                                                        <  2:663580/126‑1 (MQ=255)
                                   ttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAg                                                                                                                                        >  1:663580/1‑126 (MQ=255)
                                    tCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCt                                                                                                               <  2:1497863/150‑1 (MQ=255)
                                      ggCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGt                                                                                                                                            <  1:1574982/119‑1 (MQ=255)
                                      ggCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGt                                                                                                                                            >  2:1574982/1‑119 (MQ=255)
                                            cTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCCGCGAATCAGCCTGTGGCGATTGGTCGCGTTGATGCCCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTg                                                                                                       >  1:1023435/1‑150 (MQ=255)
                                             tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGa                                                                                                                                              <  2:2227195/110‑1 (MQ=255)
                                             tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGa                                                                                                                                              >  1:2227195/1‑110 (MQ=255)
                                             tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTgg                                                                                                      <  1:1009396/150‑1 (MQ=255)
                                                  aTCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCgctg                                                                                                 >  1:768987/1‑150 (MQ=255)
                                                       ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGc                                                                                                                                       <  1:1879321/107‑1 (MQ=255)
                                                       ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGc                                                                                                                                       >  2:1879321/1‑107 (MQ=255)
                                                       ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTgg                                                                                             <  2:438298/149‑1 (MQ=255)
                                                       ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTgg                                                                                             >  1:438298/1‑149 (MQ=255)
                                                           tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                        >  1:1286084/1‑150 (MQ=255)
                                                           tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                        <  2:768987/150‑1 (MQ=255)
                                                           tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTGGGCCAAAGGcc                                                                                                                                                <  2:1979195/94‑1 (MQ=255)
                                                           tgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTGGGCCAAAGGcc                                                                                                                                                >  1:1979195/1‑94 (MQ=255)
                                                            gtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCTAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATgg                                                                                       <  2:1400476/150‑1 (MQ=255)
                                                                gggTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCct                                                                                   >  2:1498249/1‑150 (MQ=255)
                                                                  gtggtgCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                                                                   >  1:2359064/1‑84 (MQ=255)
                                                                  gtggtgCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                                                                   <  2:2359064/84‑1 (MQ=255)
                                                                       gCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                                                                            <  2:1880412/150‑1 (MQ=255)
                                                                         ggCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGgt                                                                          <  1:2164392/150‑1 (MQ=255)
                                                                          gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                                                                            >  1:288024/1‑147 (MQ=255)
                                                                          gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                                                                            <  2:288024/147‑1 (MQ=255)
                                                                          gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGgtg                                                                         <  1:1033815/150‑1 (MQ=255)
                                                                            aGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGTGTAGTGTCCTGCTGCTGGTGtt                                                                       >  2:1159051/1‑150 (MQ=255)
                                                                            aGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGtt                                                                       <  2:1652652/150‑1 (MQ=255)
                                                                                gTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCAcc                                                                   >  1:1791086/1‑150 (MQ=255)
                                                                                 tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGgcgc                                                                                                   >  1:1464660/1‑117 (MQ=255)
                                                                                 tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGgcgc                                                                                                   <  2:1464660/117‑1 (MQ=255)
                                                                                   aaTCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAAc                                                                >  1:947792/1‑150 (MQ=255)
                                                                                     tCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACACACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGc                                                              >  2:1389191/1‑150 (MQ=255)
                                                                                      tctgtgGTGGATCGGCGCGTTGATGCGCTGATGGCAAGCACCCACAAACGGCCTTTGGCCCAAGGCAGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGcc                                                             <  2:1623438/148‑1 (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                                               >  1:2028209/1‑150 (MQ=255)
                                                                                                    gcgcATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCt                                               <  2:2252845/150‑1 (MQ=255)
                                                                                                       cATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGGCCTCACCGCCTGGCTGAc                                            >  1:166439/1‑150 (MQ=255)
                                                                                                         ttcaTGCGCTGATGGCACGCACTCACAAACGGCCTTTGGCACAAGGGCGAGTCGAGCCCCTGCCGGCGCTGTTGTGCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTAGTGTTCACCAACGCCCTCACCGCCTGGCTGAccc                                          <  2:1052878/147‑1 (MQ=255)
                                                                                                           gATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTCGTCCGCCGTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTg                                        <  1:1496006/150‑1 (MQ=255)
                                                                                                              gcgcTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGAGCTGTTGTTCGCCTTGGCCCTGGCgctg                                                                                                 <  2:2341957/90‑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                                                                                                                   <  1:2210103/69‑1 (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)
                                                                                                                             acCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTAt                      >  1:1661863/1‑150 (MQ=255)
                                                                                                                                cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGGTGCATTTTGCCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg                   <  2:1418910/150‑1 (MQ=255)
                                                                                                                                cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg                   >  2:238020/1‑150 (MQ=255)
                                                                                                                                cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg                   >  2:716426/1‑150 (MQ=255)
                                                                                                                                cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg                   >  2:844601/1‑150 (MQ=255)
                                                                                                                                cacaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCg                   >  1:510576/1‑150 (MQ=255)
                                                                                                                                 acaAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCATGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCgg                  >  2:1279948/1‑150 (MQ=255)
                                                                                                                                     tcGGCCTTGGGCCAAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                                                                   <  1:1936790/48‑1 (MQ=25)
                                                                                                                                     tcGGCCTTGGGCCAAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                                                                   >  2:1936790/2‑49 (MQ=25)
                                                                                                                                     tcGGCCTTGGGCCAAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGct              >  2:421986/2‑150 (MQ=255)
                                                                                                                                     aCGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                        <  1:1111356/76‑1 (MQ=255)
                                                                                                                                     aCGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                                                                        >  2:1111356/1‑76 (MQ=255)
                                                                                                                                     aCGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGct              >  1:406373/1‑150 (MQ=255)
                                                                                                                                      cGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGctc             <  1:844601/150‑1 (MQ=255)
                                                                                                                                       ggCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGctct            >  2:2086320/1‑150 (MQ=255)
                                                                                                                                       ggCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGctc             >  1:601580/1‑149 (MQ=255)
                                                                                                                                       ggCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGctc             <  2:601580/149‑1 (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)
                                                                                                                                               gCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACACCg      >  1:815303/1‑148 (MQ=255)
                                                                                                                                               gCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACACCg      <  2:815303/148‑1 (MQ=255)
                                                                                                                                                cccAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACACCGGCt   >  2:1847311/1‑150 (MQ=255)
                                                                                                                                                 ccAAGGCTGAGTCGAGCCTCTGCCGGAGCTGTTGTTCGACTTGGCCCTGGCGCTGCTCGGCATGGTTCTGCTGCTGGTGTTTACCAACGCCCTCACCGCCTGGCTGACCCTCGCCTCATTACTCGGCTATGCGGTGCTCTACACCGGCtt  >  2:159244/1‑150 (MQ=255)
                                                                                                                                                 ccAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGccttggccttgg                                                                                                      <  2:1851624/50‑1 (MQ=255)
                                                                                                                                                 ccAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGccttggccttgg                                                                                                      >  1:1851624/1‑50 (MQ=255)
                                                                                                                                                 ccAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTTCTCTACACCGGCtt  >  2:169631/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGCTGACCCTGGCCTCATTACTCGGTTATGCGGTGCTCTACACCGGCTT  >  NC_002947/115810‑116104

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: