Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,337,042 C→A 6.3% L137I (CTC→ATC) ‡ gstA → glutathione S‑transferase GstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,337,0420CA 6.3% 253.5 / 5.1 79L137I (CTC→ATC) ‡gstAglutathione S‑transferase GstA
Reads supporting (aligned to +/- strand):  ref base C (33/41);  new base A (2/3);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

CGAAGAACCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAACTG  >  NC_002947/1336894‑1337187
                                                                                                                                                    |                                                                                                                                                 
cGAAGAACCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCt                                                                                                                                                  <  2:97223/150‑1 (MQ=255)
   agaaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCaa                                                                                                                                               <  1:589072/150‑1 (MQ=255)
   agaaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCaa                                                                                                                                               <  1:2216740/150‑1 (MQ=255)
    gaaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAg                                                                                                                                              <  2:253762/150‑1 (MQ=255)
     aaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAgg                                                                                                                                             <  2:2396768/150‑1 (MQ=255)
      aCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGATTGCCGCAGGCCGAACAACCTTTGGCCCGGCAGCCTCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGGGATTTCCCGCGCCCACACCTTGAGCAAGGt                                                                                                                                            >  1:420783/1‑150 (MQ=255)
       ccTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTg                                                                                                                                           >  1:935870/1‑150 (MQ=255)
        cTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGa                                                                                                                                          <  1:1105975/150‑1 (MQ=255)
         tGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGAt                                                                                                                                         <  2:1926529/150‑1 (MQ=255)
         tGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGAGCAAGGTgtg                                                                                                                                         <  2:420783/150‑3 (MQ=255)
          gccgccgccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATc                                                                                                                                        <  1:589562/150‑1 (MQ=255)
          gccgccgccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATc                                                                                                                                        <  1:2143255/150‑1 (MQ=255)
          gccgccgccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGAGCAAGGTgtgg                                                                                                                                        <  2:1240686/150‑4 (MQ=255)
             gccgccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGAt                                                                                                                                     <  2:1855208/150‑1 (MQ=255)
                     ccGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGAGCAAGGTgtg                                                                                                                                         >  1:2281365/1‑136 (MQ=255)
                     ccGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGAGCAAGGTgtg                                                                                                                                         <  2:2281365/138‑3 (MQ=255)
                                gCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATc                                                                                                                                        <  2:725878/128‑1 (MQ=255)
                                gCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATc                                                                                                                                        >  1:725878/1‑128 (MQ=255)
                                 cTGGCAATCGGCCAACGCAGGGCCGCTGGGCAACAGCCCGGCGGCCGCTCCCCTTGTCACGGTGTTCGGGGCTTCGTTCAACACCCACAAATTGACTTCCCGCGCCCACCCCTTGCTCGAGGTCATCGATGCCCGGTTGACCAATACGcc                                                                                                                 >  2:280554/1‑150 (MQ=255)
                                    gCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCAtt                                                                                                              >  1:826650/1‑150 (MQ=255)
                                     cTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTc                                                                                                             <  2:826650/150‑1 (MQ=255)
                                     cTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTc                                                                                                             <  1:2133801/150‑1 (MQ=255)
                                      tATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAg                                                                                                                                              >  2:2380383/1‑116 (MQ=255)
                                      tATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAg                                                                                                                                              <  1:2380383/116‑1 (MQ=255)
                                      tATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTcc                                                                                                            >  2:2242206/1‑150 (MQ=255)
                                            gggccaaaGGCCAGCGGGCCTGCGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACgg                                                                                                      >  1:857488/8‑150 (MQ=255)
                                             cgccgcAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACggg                                                                                                     <  1:2186979/150‑1 (MQ=255)
                                              gccgcAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGgc                                                                                                    >  2:412492/1‑150 (MQ=255)
                                              gccgcAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGgc                                                                                                    >  1:1345505/1‑150 (MQ=255)
                                              gccaaaGGCCAGCGGGCCTGCGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGgc                                                                                                    >  2:187898/6‑150 (MQ=255)
                                               ccgcAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAgg                                                                                                                                             <  2:1947407/108‑1 (MQ=255)
                                               ccgcAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAgg                                                                                                                                             >  1:1947407/1‑108 (MQ=255)
                                                   aGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGcac                                                                                               <  1:412492/150‑1 (MQ=255)
                                                     gCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCacac                                                                                             <  2:2369308/150‑1 (MQ=255)
                                                       ccGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTCCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCCcaccc                                                                                           >  1:623723/1‑150 (MQ=255)
                                                               cTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGc                                                                                   >  1:2144132/1‑150 (MQ=255)
                                                                  tGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGCTGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGa                                                                                <  2:941973/150‑1 (MQ=255)
                                                                  tGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGa                                                                                <  2:911342/150‑1 (MQ=255)
                                                                  cggCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGa                                                                                <  2:857488/149‑1 (MQ=255)
                                                                    gCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGt                                                                                                                                            <  1:913653/88‑1 (MQ=255)
                                                                    gCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGt                                                                                                                                            >  2:913653/1‑88 (MQ=255)
                                                                               gCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAactac                                                                   >  2:2409009/1‑150 (MQ=255)
                                                                               gCTCCCTTGGTCACGGTGTCCGGGGCTTCGTCCAATACCGCCGAAGTAATTTCCCGCGCCAACCCCTTGCTAAAGGTGATCGATGCCAAGTTGGCCAATCCGCCATTCCTGCCGGGAAGCACA‑CCACCATAGCGGACATTGCCAactact                                                                  <  1:490761/150‑1 (MQ=255)
                                                                                 tCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTc                                                                 <  1:187898/150‑1 (MQ=255)
                                                                                     cTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTAc                                                             <  2:1474359/150‑1 (MQ=255)
                                                                                         tCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAACACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCg                                                         >  1:101068/1‑150 (MQ=255)
                                                                                             ggTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCa                                                     >  1:2199824/1‑150 (MQ=255)
                                                                                                 ttCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGcc                                                 >  2:496912/1‑150 (MQ=255)
                                                                                                  tCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGccc                                                <  2:2087991/150‑1 (MQ=255)
                                                                                                      ggCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGa                                            >  1:586970/1‑150 (MQ=255)
                                                                                                       gCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGaa                                           >  2:446427/1‑150 (MQ=255)
                                                                                                            gTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCAACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCaa                                      >  2:1405860/1‑150 (MQ=255)
                                                                                                             ttCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAAc                                     >  1:1683263/1‑150 (MQ=255)
                                                                                                                aaTACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGc                                                        <  2:758875/128‑1 (MQ=255)
                                                                                                                aaTACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGc                                                        >  1:758875/1‑128 (MQ=255)
                                                                                                                 aTACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAg                                                                                                                                              >  2:1610052/1‑41 (MQ=255)
                                                                                                                 aTACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAg                                                                                                                                              <  1:1610052/41‑1 (MQ=255)
                                                                                                                  tACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTc                                <  1:496912/150‑1 (MQ=255)
                                                                                                                           aGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCaa                                      <  2:357902/135‑1 (MQ=255)
                                                                                                                           aGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCaa                                      >  1:357902/1‑135 (MQ=255)
                                                                                                                           aGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGcc                       <  2:724944/150‑1 (MQ=255)
                                                                                                                                 ttCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGc                 <  2:1345505/150‑1 (MQ=255)
                                                                                                                                   cccGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCAcacccac                                                                                         >  1:631471/1‑76 (MQ=255)
                                                                                                                                   cccGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCAcacccac                                                                                         <  2:631471/76‑1 (MQ=255)
                                                                                                                                   cccGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCa               <  1:2409009/150‑1 (MQ=255)
                                                                                                                                    ccGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAACGGACATTGCCAACTACACGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGCGCCCTCTGCCaa              >  2:388001/1‑150 (MQ=255)
                                                                                                                                     cgcgCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAAt             <  1:446427/150‑1 (MQ=255)
                                                                                                                                     cgcgCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCAACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAAt             <  1:1405860/150‑1 (MQ=255)
                                                                                                                                      gcgcCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTc                                                      >  1:1427296/1‑108 (MQ=255)
                                                                                                                                      gcgcCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTc                                                      <  2:1427296/108‑1 (MQ=255)
                                                                                                                                         cccACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTgc         <  1:2242206/150‑1 (MQ=255)
                                                                                                                                            acacCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCg                               <  1:333660/125‑1 (MQ=255)
                                                                                                                                            acacCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCg                               >  2:333660/1‑125 (MQ=255)
                                                                                                                                             cacCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGAAACGAAGAAGTACTCGTACATCGCTCACGCCAACGAAGGCAACGTGAAGCTGGAGCCCGATGCCAATGTGCGCaa     >  1:2223241/1‑150 (MQ=255)
                                                                                                                                              acCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAAc    >  2:472976/1‑150 (MQ=255)
                                                                                                                                                cTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTAc                                                             >  1:1388063/1‑91 (MQ=255)
                                                                                                                                                cTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTAc                                                             <  2:1388063/91‑1 (MQ=255)
                                                                                                                                                cTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAActg  >  1:31298/1‑150 (MQ=255)
                                                                                                                                                cTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAActg  <  1:472976/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGAAGAACCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAACTG  >  NC_002947/1336894‑1337187

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