Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,218,196 A→G H290R (CAC→CGC)  ycgG → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,218,1960AG98.0% 342.5 / ‑4.5 98H290R (CAC→CGC) ycgGputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (46/50);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAATGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCACGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCCAG  >  NC_000913/1218099‑1218292
                                                                                                 |                                                                                                
aataatGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTcgcg                                                                                                <  2:2295172/100‑1 (MQ=255)
aataatGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTcgcg                                                                                                <  2:2295028/100‑1 (MQ=255)
 ataatGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCgg                                                                                               <  1:647683/100‑1 (MQ=255)
 ataatGGGGAGATCGTCCCTTTTTACCAACCTGTGATAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCgg                                                                                               <  1:744380/100‑1 (MQ=255)
  taatGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGt                                                                                              >  1:1075977/1‑100 (MQ=255)
  taatGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCACGGt                                                                                              >  1:2732547/1‑100 (MQ=255)
   aatGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTg                                                                                             >  1:2494504/1‑100 (MQ=255)
     tGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGa                                                                                           <  1:2172260/100‑1 (MQ=255)
      ggggAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGat                                                                                          >  1:2017455/1‑100 (MQ=255)
      ggggAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGat                                                                                          >  1:2796728/1‑100 (MQ=255)
       gggAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGAAATTGCGGGGAGTTGAGGTGTTCGCCCGCTGGAAACACCCTCGCGGTGGat                                                                                          >  1:2411938/1‑99 (MQ=255)
        ggAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGatat                                                                                        >  1:1197799/1‑100 (MQ=255)
        ggAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGatat                                                                                        >  1:2786558/1‑100 (MQ=255)
          agaTCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGatatat                                                                                      <  1:1149673/100‑1 (MQ=255)
           gaTCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGCTGTTAGCCCGCTGGAAACAACCTCACGGTGGctatata                                                                                     >  1:507460/1‑100 (MQ=255)
           gaTCGGCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGatatata                                                                                     <  2:1763698/100‑1 (MQ=255)
            aTCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGatatatat                                                                                    >  1:2314582/1‑100 (MQ=255)
                 cccTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCAccc                                                                               >  1:472359/1‑100 (MQ=255)
                  ccTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCcg                                                                              >  2:2686512/1‑100 (MQ=255)
                    tttttACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGGTGATATATATccaccgcg                                                                            <  1:1917403/100‑8 (MQ=255)
                      tttACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGc                                                                          <  1:1160314/100‑1 (MQ=255)
                      tttACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGc                                                                          <  1:1461224/100‑1 (MQ=255)
                       ttACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGGTGATATATATccaccgcgagg                                                                         <  1:2287639/100‑11 (MQ=255)
                         aCCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCAtt                                                                       >  1:2241798/1‑100 (MQ=255)
                          ccAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTATCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCAttt                                                                      <  1:917008/100‑1 (MQ=255)
                          ccAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCAttt                                                                      <  1:405862/100‑1 (MQ=255)
                           cAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTa                                                                     <  2:2747174/100‑1 (MQ=255)
                            aaCCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTAt                                                                    >  1:1843617/1‑100 (MQ=255)
                             aCCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTAtt                                                                   >  1:1997663/1‑100 (MQ=255)
                               cTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTcc                                                                 <  2:523428/100‑1 (MQ=255)
                               cTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTcc                                                                 <  1:2638552/100‑1 (MQ=255)
                                 gtgGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCAc                                                               >  2:2622150/1‑100 (MQ=255)
                                   ggTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACtt                                                             >  2:2056628/1‑100 (MQ=255)
                                    gTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTg                                                            >  1:1567783/1‑100 (MQ=255)
                                    gTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTg                                                            >  1:1188913/1‑100 (MQ=255)
                                    gTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTg                                                            <  1:2416345/100‑1 (MQ=255)
                                     tAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGc                                                           >  2:2069985/1‑100 (MQ=255)
                                           gTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGaaaa                                                     >  2:2019157/1‑100 (MQ=255)
                                            tCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGaaaaa                                                    <  1:2568577/100‑1 (MQ=255)
                                            tCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGaaaaa                                                    <  2:29609/100‑1 (MQ=255)
                                              gggAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGaaaaa                                                    >  2:1383057/1‑98 (MQ=255)
                                               ggAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATcc                                                 <  1:2736954/100‑1 (MQ=255)
                                               ggAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATcc                                                 >  1:92920/1‑100 (MQ=255)
                                               ggAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATcc                                                 >  1:187271/1‑100 (MQ=255)
                                                gAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCg                                                <  2:2785991/100‑1 (MQ=255)
                                                 aaGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCgg                                               >  2:122342/1‑100 (MQ=255)
                                                  aGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGa                                              <  1:2013307/100‑1 (MQ=255)
                                                  aGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGa                                              >  2:2048331/1‑100 (MQ=255)
                                                  aGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGa                                              <  2:541474/100‑1 (MQ=255)
                                                     gACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTa                                           <  1:1919704/100‑1 (MQ=255)
                                                     gACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTa                                           >  1:782105/1‑100 (MQ=255)
                                                       cATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAAt                                         >  2:1773866/1‑100 (MQ=255)
                                                        aTTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATc                                        <  1:804491/100‑1 (MQ=255)
                                                         ttGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCg                                       <  1:856918/100‑1 (MQ=255)
                                                           gCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGtt                                     >  1:492938/1‑100 (MQ=255)
                                                             ggggggTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTcc                                   >  2:2650459/1‑100 (MQ=255)
                                                              gggAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCg                                  <  1:1383062/100‑1 (MQ=255)
                                                              gggAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCg                                  <  2:2590032/100‑1 (MQ=255)
                                                                   ttGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATGTCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTAc                             <  2:1843604/100‑1 (MQ=255)
                                                                   ttGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTAc                             >  1:1716160/1‑100 (MQ=255)
                                                                   ttGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTAc                             >  1:1965434/1‑100 (MQ=255)
                                                                   ttGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTAc                             <  1:2141450/100‑1 (MQ=255)
                                                                    tGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACg                            >  1:2820312/1‑100 (MQ=255)
                                                                    tGAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACg                            >  2:749459/1‑100 (MQ=255)
                                                                     gAGGTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGc                           >  1:2068386/1‑100 (MQ=255)
                                                                       ggTGTTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCaa                         >  1:2617644/1‑100 (MQ=255)
                                                                        gtgtTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCaaa                        <  1:1045692/100‑1 (MQ=255)
                                                                         tgtTAGCCCGCTGGACACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAg                       >  2:998979/1‑100 (MQ=255)
                                                                         tgtTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAg                       <  2:2796707/100‑1 (MQ=255)
                                                                         tgtTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAg                       >  1:2282560/1‑100 (MQ=255)
                                                                          gtTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGc                      <  1:2311559/100‑1 (MQ=255)
                                                                          gtTAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGc                      >  1:729196/1‑100 (MQ=255)
                                                                           ttAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGcc                     <  1:2057071/100‑1 (MQ=255)
                                                                           ttAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGACATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGcc                     <  2:1626953/100‑1 (MQ=255)
                                                                            tAGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCt                    <  2:2241785/100‑1 (MQ=255)
                                                                             aGCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCtg                   >  2:64108/1‑100 (MQ=255)
                                                                              gCCCGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCtga                  >  1:2851847/1‑100 (MQ=255)
                                                                               cccGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCtgat                 <  1:1372890/100‑1 (MQ=255)
                                                                               cccGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCtgat                 <  1:868139/100‑1 (MQ=255)
                                                                                ccGCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCtgatg                <  2:187270/100‑1 (MQ=255)
                                                                                ccGCTGGAAACAACCTCGCGGTCGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCtgatg                >  2:497295/1‑100 (MQ=255)
                                                                                  gCTGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGaa              <  2:472357/100‑1 (MQ=255)
                                                                                    tGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATc            <  1:564152/100‑1 (MQ=255)
                                                                                    tGGAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATc            <  2:876890/100‑1 (MQ=255)
                                                                                     ggAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCa           <  1:2622169/100‑1 (MQ=255)
                                                                                     ggAAACAACCGCGCGGGGGATATATATCACGCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCa           <  1:2118438/100‑1 (MQ=255)
                                                                                      gAAACAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAg          >  1:2729969/1‑100 (MQ=255)
                                                                                       aaaCAACCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAgg         <  2:1997649/100‑1 (MQ=255)
                                                                                        aacaacCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGt        <  1:1113911/100‑1 (MQ=255)
                                                                                          caacCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTg      <  2:247875/100‑1 (MQ=255)
                                                                                          caacCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCCGATGAATCAGGTTg      <  1:2019171/100‑1 (MQ=255)
                                                                                           aacCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATTCATCAggcttt     <  1:2050633/100‑5 (MQ=255)
                                                                                            acCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGcc    <  2:2314567/100‑1 (MQ=255)
                                                                                            acCTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGcc    <  1:1394202/100‑1 (MQ=255)
                                                                                             ccTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCca   >  2:2628037/1‑100 (MQ=255)
                                                                                             ccTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCca   >  1:227542/1‑100 (MQ=255)
                                                                                             ccTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCca   <  1:2069998/100‑1 (MQ=255)
                                                                                              cTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCcag  >  1:2788272/1‑100 (MQ=255)
                                                                                              cTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCcag  >  2:229297/1‑100 (MQ=255)
                                                                                              cTCGCGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCcag  >  1:59322/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AATAATGGGGAGATCGTCCCTTTTTACCAACCTGTGGTAAATGGTCGGGAAGGGACATTGCGGGGAGTTGAGGTGTTAGCCCGCTGGAAACAACCTCACGGTGGATATATATCACCCGCGGCATTTATTCCACTTGCTGAAAAATCCGGATTAATCGTTCCGCTTACGCAAAGCCTGATGAATCAGGTTGCCAG  >  NC_000913/1218099‑1218292

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