Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,693 C→G 16.3% P268R (CCG→CGG)  cca → fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,6930CG16.3% 301.5 / 44.5 115P268R (CCG→CGG) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base C (51/45);  new base G (8/11);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 117 total aligned reads displayed.

GCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACT  >  NC_000913/3202595‑3202789
                                                                                                  |                                                                                                
gCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcg                                                                                                 <  2:530167/100‑1 (MQ=255)
  tctcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgcg                                                                                               <  1:2177065/100‑1 (MQ=255)
    tcgaTGGCGGCGATGCTGAGTCCGCAGGGCGATGTCCGTTTCGCGACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  1:1178361/97‑1 (MQ=37)
    tcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAATGGGCTGACGCCGCCAGACCTCTGGCCGCGtc                                                                                             <  1:1727115/100‑1 (MQ=255)
    tcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  1:1255341/100‑1 (MQ=38)
     cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtca                                                                                            >  1:1017811/1‑100 (MQ=255)
      aaTGGCGGTCATGCTGAGTCCGCAGGACGATGTCCGTTGCGCACCTTTATGCCACGATCTCGGTGAACGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  2:1682537/98‑1 (MQ=25)
      aaTGGCGGCGATGCTGCGGCCGGAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  2:687565/98‑1 (MQ=37)
      aaTGGCGACGATGCTGAGCCCGCAGGTCTATGTCCGTTTCCCAACTTTATGCCACGACCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  2:619089/98‑1 (MQ=25)
       aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                                                          >  2:2534444/1‑100 (MQ=255)
        tGGCGGCGGTGCTTAGTCCGCATGTCGATGGCCGTTCCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  2:1047652/96‑1 (MQ=25)
        tGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAGCTTTATGCCACGAGCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  2:606481/96‑1 (MQ=37)
        tGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagcc                                                                                         <  1:1083581/100‑4 (MQ=37)
        tGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                             <  1:1178336/96‑1 (MQ=38)
         ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTcatcat                                                                                        >  1:938545/1‑100 (MQ=255)
          gcggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatg                                                                                       >  2:2068222/1‑100 (MQ=255)
           cggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatgg                                                                                      >  2:626483/1‑100 (MQ=255)
           cggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatgg                                                                                      >  2:323110/1‑100 (MQ=255)
            ggcgATGCTGAGTCCGCAGGTCGACGTCCGTTTCGGAACTTGAGGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTcatcat                                                                                        <  1:682874/97‑1 (MQ=255)
             gcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                    >  2:2510513/1‑100 (MQ=255)
               gATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                    <  2:1540645/98‑1 (MQ=255)
                   cTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGcc                                                                              >  1:1793422/1‑100 (MQ=255)
                         ccGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcgg                                                                        >  2:1474250/1‑100 (MQ=255)
                         ccGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcg                                                                         >  2:1188508/1‑99 (MQ=255)
                          cGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                       <  1:147270/100‑1 (MQ=255)
                           gCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                      >  2:1297206/1‑100 (MQ=255)
                           gCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                      <  1:805314/100‑1 (MQ=255)
                            cAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                     <  2:1407613/100‑1 (MQ=255)
                            cAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                     <  1:892060/100‑1 (MQ=255)
                            cAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                     >  1:1696107/1‑100 (MQ=255)
                             aGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtgt                                                                    <  1:1839103/100‑1 (MQ=255)
                              ggTCGATGTCCGTTTCGCAACTTTATGCCACGATCTAGGTAGAGGGCTGACGCCGCCAGAACTCTGGCCGCGTGATCATGCTCATGGCCCGGCGGgtgt                                                                    >  2:2456883/1‑99 (MQ=255)
                              ggTCGATGTCCGGTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtgt                                                                    <  2:1482146/99‑1 (MQ=255)
                               gTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                  <  1:1066100/100‑1 (MQ=255)
                               gTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                  <  1:2347280/100‑1 (MQ=255)
                               gTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                  <  1:2203649/100‑1 (MQ=255)
                               gTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                  <  1:1762546/100‑1 (MQ=255)
                                 cGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAg                                                                <  2:768693/100‑1 (MQ=255)
                                  gATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                               <  2:2168183/100‑1 (MQ=255)
                                   aTGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGtt                                                              <  1:1273495/100‑1 (MQ=255)
                                      tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                           >  1:2483582/1‑100 (MQ=255)
                                        cGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTgg                                                         <  2:2074889/100‑1 (MQ=255)
                                        cGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTgg                                                         <  1:237818/100‑1 (MQ=255)
                                        cGTTGCGCACCTTTATGCCCCGATCTCGGTAAAGGGCTGACGCCGCCAGAACGCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                           <  2:1211475/98‑1 (MQ=255)
                                          tttCGCAACTTTAGGCCACGATCTCGGTAAAGGGCTGTCGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaa                                                       <  2:1587705/100‑1 (MQ=255)
                                            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaaca                                                     >  1:1758900/1‑100 (MQ=255)
                                            tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaaca                                                     >  2:1285286/1‑100 (MQ=255)
                                             cGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                    >  1:2457308/1‑100 (MQ=255)
                                             cGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                    >  2:409149/1‑100 (MQ=255)
                                             cGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                    >  1:2361281/1‑100 (MQ=255)
                                             cGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                    >  1:414572/1‑100 (MQ=255)
                                             cGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaaca                                                     <  2:643714/99‑1 (MQ=255)
                                              gCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaac                                                   >  1:1546029/1‑100 (MQ=255)
                                              gCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                    >  1:169281/1‑99 (MQ=255)
                                                aaCTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTa                                                 <  2:2381121/100‑1 (MQ=255)
                                                  cTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATg                                               <  2:1045777/100‑1 (MQ=255)
                                                    ttATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTCGTGGAACAACTATGc                                              >  2:390613/1‑99 (MQ=255)
                                                    ttATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGcc                                             <  1:1161557/100‑1 (MQ=255)
                                                      aTGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAg                                           <  1:1090368/100‑1 (MQ=255)
                                                       tGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGCTAGTGGAACAACTATGCCAg                                           >  2:2286534/1‑99 (MQ=255)
                                                         ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGt                                        >  1:2145568/1‑100 (MQ=255)
                                                         ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGt                                        >  2:1966098/1‑100 (MQ=255)
                                                         ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGt                                        >  1:1063451/1‑100 (MQ=255)
                                                         ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGt                                        >  2:904836/1‑100 (MQ=255)
                                                         ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGt                                        >  2:1597498/1‑100 (MQ=255)
                                                           aCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCt                                      <  1:2477946/100‑1 (MQ=255)
                                                            cGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTg                                     <  1:1943460/100‑1 (MQ=255)
                                                             gATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTgc                                    >  2:82925/1‑100 (MQ=255)
                                                             gATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTgc                                    <  2:1310347/100‑1 (MQ=255)
                                                               nctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTgcgc                                  >  1:968239/2‑100 (MQ=255)
                                                                  cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTg                                     <  2:2437637/94‑1 (MQ=255)
                                                                   ggTAAAGGGCTGATGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGc                              <  1:844218/100‑1 (MQ=255)
                                                                   ggTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGc                              >  1:2194233/1‑100 (MQ=255)
                                                                   ggTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGc                              <  2:2035766/100‑1 (MQ=255)
                                                                    gTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGcc                             <  1:1276773/100‑1 (MQ=255)
                                                                      aaaGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTg                               <  2:523947/96‑1 (MQ=255)
                                                                         ggggTGACGCCGGCAGAACTCTGGCCGCTTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGTGTGCCAAATg                        <  2:2383938/100‑1 (MQ=255)
                                                                          tgatGACGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGa                       >  1:1056763/4‑100 (MQ=37)
                                                                           gCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGaa                      <  1:1605074/100‑1 (MQ=255)
                                                                            cTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGaaa                     <  2:1063443/100‑1 (MQ=255)
                                                                             tGACGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAAt                    >  1:1604251/1‑100 (MQ=38)
                                                                              gACGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAAtt                   >  2:1442541/1‑100 (MQ=38)
                                                                                cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg                 >  1:970883/1‑100 (MQ=255)
                                                                                cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg                 >  2:219805/1‑100 (MQ=255)
                                                                                cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg                 >  1:528752/1‑100 (MQ=255)
                                                                                cgcggcCAGAGTTCTGGCGGCGACATCATCGTCATGGCCCGGCGGGTGTTAAGTTAGCAGAACAACGATGCCAGCGTATGCGCGTGCCACACGAAAtcg                  >  2:1675804/1‑97 (MQ=255)
                                                                                cgccgcCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg                 >  2:33909/1‑100 (MQ=255)
                                                                                cgccgcCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg                 <  1:1898567/100‑1 (MQ=255)
                                                                                cgccgcCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTc                  >  2:1630900/1‑99 (MQ=255)
                                                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgc                <  1:681311/100‑1 (MQ=255)
                                                                                 gccgccAGAACTCTGGCCGCGTCATCATGGTCAGGGCCCGGCGGGTGTTGAGTTAGTGGAACAACTATGCTAGCGTCTGCGCGTGCCAAATGAAATTcgc                >  2:1543718/1‑100 (MQ=255)
                                                                                  ccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgcg               >  1:2508825/1‑100 (MQ=255)
                                                                                    gccAGAAGTCTGGCCGCGTCATCATGGTCATGGTCCGGCGGGTGTTAAGTTAGTGGAAAAGCTGTGCCAGCGTCTGCGCGTGCGACATGAAATTc                  <  2:1247226/95‑1 (MQ=255)
                                                                                    gccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAt             >  1:2140593/1‑100 (MQ=255)
                                                                                          aCTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGcc       <  2:1064910/100‑1 (MQ=255)
                                                                                           ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCa      <  2:2087315/100‑1 (MQ=255)
                                                                                              tGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGAc   >  2:1699567/1‑100 (MQ=255)
                                                                                               ggCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACt  >  2:1128196/1‑100 (MQ=255)
                                                                                               ggCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACt  >  1:895227/1‑100 (MQ=255)
                                                                                               ggCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACt  >  2:692943/1‑100 (MQ=255)
                                                                                                  |                                                                                                
GCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACT  >  NC_000913/3202595‑3202789

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