Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 138,700 G→A 6.0% R534Q (CGG→CAG)  PP_0131 → membrane‑anchored phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947138,7000GA 6.0% 245.5 / 3.1 82R534Q (CGG→CAG) PP_0131membrane‑anchored phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (45/32);  new base A (2/3);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGGGTAATTGCCGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTT  >  NC_002947/138552‑138847
                                                                                                                                                    |                                                                                                                                                   
cGGGTAATTGCCGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGcggcgg                                                                                                                                                    <  2:1997341/150‑1 (MQ=255)
cGGGTAATTGCCGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGcggcgg                                                                                                                                                    <  2:363864/150‑1 (MQ=255)
   gTAATTGCCGAAGGAATCGAGTCCACGCAGGAGGCGGTGCTGCTGAACAGGGAGGGGGTCAATTATGGCCCGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGTTGACTCGTGGGCACCTGCCGCGGCGGgtg                                                                                                                                                 >  1:1816327/1‑150 (MQ=255)
      aTTGCCGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAg                                                                                                                                              <  2:457889/150‑1 (MQ=255)
      aTTGCCGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAg                                                                                                                                              <  2:346347/150‑1 (MQ=255)
       ttGCCGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCCGCAGCGGCAGgtgccca                                                                                                                                             <  1:124578/150‑5 (MQ=255)
           cGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGcgcg                                                                                                                                         >  2:1174365/1‑150 (MQ=255)
           cGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGcgcg                                                                                                                                         >  1:1557386/1‑150 (MQ=255)
            gAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGt                                                                                                                                        >  2:2168131/1‑150 (MQ=255)
            gAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGt                                                                                                                                        <  1:1219228/150‑1 (MQ=255)
            gAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGt                                                                                                                                        <  1:2288451/150‑1 (MQ=255)
              aGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCg                                                                                                                                      >  2:224826/1‑150 (MQ=255)
                  aTCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGc                                                                                                                                  <  2:188477/150‑1 (MQ=255)
                   tCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCg                                                                                                                                 >  1:1249721/1‑150 (MQ=255)
                       gTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCa                                                                                                                             <  2:1043640/150‑1 (MQ=255)
                        ttCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAg                                                                                                                            >  2:1179366/1‑150 (MQ=255)
                           gAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGAcc                                                                                                                         >  2:1237295/1‑150 (MQ=255)
                           gAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGAcc                                                                                                                         <  1:2198235/150‑1 (MQ=255)
                              gACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAgg                                                                                                                                             >  2:605581/1‑127 (MQ=255)
                              gACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAgg                                                                                                                                             <  1:605581/127‑1 (MQ=255)
                                 cAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAAc                                                                                                                   <  2:1506834/150‑1 (MQ=255)
                                         gctgctGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTa                                                                                                           <  1:1174365/150‑1 (MQ=255)
                                         gctgctGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGCGCCTCTTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTa                                                                                                           >  2:663423/1‑150 (MQ=255)
                                           tgctgAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGa                                                                                                         >  2:2216726/1‑150 (MQ=255)
                                              tgAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACt                                                                                                      >  1:385445/1‑150 (MQ=255)
                                                aaCAGCGAGGGGGTCAATTATGGCCAGTGCTGGTTGTTTGCCCGGCCGTCGCTCGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTgg                                                                                                    <  2:1816327/150‑1 (MQ=255)
                                                aaCAGCGAGGGGGGCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGNGTGNAGGCNCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTgg                                                                                                    <  1:1818411/150‑1 (MQ=255)
                                                        gggggTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTg                                                                                            <  2:1249721/150‑1 (MQ=255)
                                                           ggTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTgg                                                                                                    <  1:2286048/139‑1 (MQ=255)
                                                           ggTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTgg                                                                                                    >  2:2286048/1‑139 (MQ=255)
                                                               aaTTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGc                                                                                     >  1:1598510/1‑150 (MQ=255)
                                                                 ttATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCcg                                                                                   >  1:730660/1‑150 (MQ=255)
                                                                 ttATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCCGCCGCGGCAGgtgcc                                                                                                                                               >  2:1613060/1‑88 (MQ=255)
                                                                 ttATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCCGCCGCGGCAGgtgcc                                                                                                                                               <  1:1613060/90‑3 (MQ=255)
                                                                       ccAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGaaaa                                                                             >  1:1054067/1‑150 (MQ=255)
                                                                         aGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAaca                                                                           <  1:1662557/150‑1 (MQ=255)
                                                                            gCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATc                                                                                                                              >  1:1183304/1‑96 (MQ=255)
                                                                            gCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATc                                                                                                                              <  2:1183304/96‑1 (MQ=255)
                                                                             cTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGCACTAGGTCCCCTGCCCCGGCCGCGCAAACCCGcc                                                                       >  2:419189/1‑150 (MQ=255)
                                                                             cTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGcc                                                                       >  1:824959/1‑150 (MQ=255)
                                                                              tGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGTCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGccc                                                                      >  2:346504/1‑150 (MQ=255)
                                                                              tGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGccc                                                                      >  2:1798479/1‑150 (MQ=255)
                                                                                    tttGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATgcgc                                                                <  1:2216726/150‑1 (MQ=255)
                                                                                        cccGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCg                                                                                                                                      >  1:1583704/1‑76 (MQ=255)
                                                                                        cccGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCg                                                                                                                                      <  2:1583704/76‑1 (MQ=255)
                                                                                        cccGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGcgcccgc                                                            >  1:618821/1‑150 (MQ=255)
                                                                                         tcagCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGcgcccgcc                                                           >  1:1616349/4‑150 (MQ=255)
                                                                                         ccGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGcgcccgcc                                                           >  1:1680419/1‑150 (MQ=255)
                                                                                            gCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTg                                                        >  1:82337/1‑150 (MQ=255)
                                                                                             ccGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGc                                                                                                                                  >  2:2313766/1‑75 (MQ=255)
                                                                                             ccGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGc                                                                                                                                  <  1:2313766/75‑1 (MQ=255)
                                                                                                cTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGc                                                    >  1:201734/1‑150 (MQ=255)
                                                                                                 tGAACGCCCGGCAATTTCCCGAGTCGACAATTATTGGACCCCTGCCTCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCAAGCGCCCGCCGTGCAGCt                                                   <  2:1191829/150‑1 (MQ=255)
                                                                                                        ccGGCAGCTTGCCGAGACGGCGCCCCGTGGGCACCTGCCGCGGCGGGTTTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAAt                                            >  2:616101/1‑150 (MQ=255)
                                                                                                         cGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAAtt                                           <  2:1598510/150‑1 (MQ=255)
                                                                                                           gCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGAGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTc                                         >  2:1811667/1‑150 (MQ=255)
                                                                                                            cAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTcc                                        <  1:2168131/150‑1 (MQ=255)
                                                                                                             aGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCt                                       <  2:1189198/150‑1 (MQ=255)
                                                                                                             aGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCt                                       <  1:1600409/150‑1 (MQ=255)
                                                                                                             aGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACCCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCt                                       >  1:1043564/1‑150 (MQ=255)
                                                                                                              gTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGCAGCGATCAGTCCGTCAACGCCATGTCGGTCTGGGCCCACTGCCCCTACCGCCAATCCATGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCtt                                      <  1:1798479/150‑1 (MQ=255)
                                                                                                               tttGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATACGCCGGCCGTGCATCTGCACAATTTCCTTg                                     >  2:462928/1‑150 (MQ=255)
                                                                                                                    cGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAg                                >  1:995177/1‑150 (MQ=255)
                                                                                                                     gAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAgc                               <  1:1237295/150‑1 (MQ=255)
                                                                                                                         tGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGa                           <  2:618821/150‑1 (MQ=255)
                                                                                                                         tGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGa                           <  1:1376981/150‑1 (MQ=255)
                                                                                                                            tGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGcc                        >  1:360367/1‑150 (MQ=255)
                                                                                                                                   tGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCcccagccc                 >  2:930432/1‑150 (MQ=255)
                                                                                                                                    gggCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCcccagccca                >  2:1622758/1‑150 (MQ=255)
                                                                                                                                      gCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGa                           >  2:2232498/1‑137 (MQ=255)
                                                                                                                                      gCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGa                           <  1:2232498/137‑1 (MQ=255)
                                                                                                                                      gCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCccagcccagc              <  2:1054067/150‑1 (MQ=255)
                                                                                                                                      gCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCccagcccagc              >  1:88168/1‑150 (MQ=255)
                                                                                                                                       cACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCa                                                                                                               <  1:2182993/52‑1 (MQ=255)
                                                                                                                                       cACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCa                                                                                                               >  2:2182993/1‑52 (MQ=255)
                                                                                                                                          ccgCCGCGGCAGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGc                                                                        <  2:356537/86‑1 (MQ=255)
                                                                                                                                          ccgCCGCGGCAGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGc                                                                        >  1:356537/3‑88 (MQ=255)
                                                                                                                                           tGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCGGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGcgccgc         >  2:1012779/1‑150 (MQ=255)
                                                                                                                                           tGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGa                           >  2:1698647/1‑132 (MQ=255)
                                                                                                                                           tGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGa                           <  1:1698647/132‑1 (MQ=255)
                                                                                                                                           tGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGcgccgc         >  1:798895/1‑150 (MQ=255)
                                                                                                                                           tGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGcgccgc         >  1:353298/1‑150 (MQ=255)
                                                                                                                                            gCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCgccgcc        >  2:317382/1‑150 (MQ=255)
                                                                                                                                             ccGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCt       >  1:827463/1‑150 (MQ=255)
                                                                                                                                              ctcggcggGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCtt      <  1:2303588/148‑1 (MQ=255)
                                                                                                                                                  gcggGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCtt  <  2:353298/150‑1 (MQ=255)
                                                                                                                                                  gcggGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCtt  <  1:1622758/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGGTAATTGCCGAAGGAATCGAGTTCGAGGACCAGGCGGTGCTGCTGAACAGCGAGGGGGTCAATTATGGCCAGGGCTGGTTGTTTGCCCGGCCGCTGAACGCCCGGCAGTTTGCCGAGTTGGTGACTCGTGGGCACCTGCCGCGGCGGGTGTAGGCGCGTCGAAGCGATCAGACCGGCAACGCCATGTAGAACTGGGTCCCCTGCCCCGGCCGCGAAAACACGCCCATGCGCCCGCCGTGCAGCTGCACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTT  >  NC_002947/138552‑138847

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