Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 729,698 A→G H298R (CAC→CGC)  rhsC → RhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913729,6980AG92.1% 192.9 / 11.4 89H298R (CAC→CGC) rhsCRhsC
Reads supporting (aligned to +/- strand):  ref base A (7/0);  new base G (44/38);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 1.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

GGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAC  >  NC_000913/729667‑729786
                               |                                                                                         
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             <  1:1255058/93‑1 (MQ=14)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             >  2:7443319/1‑93 (MQ=14)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             <  2:4383993/93‑1 (MQ=14)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             <  1:7748727/93‑1 (MQ=14)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             <  1:4487765/93‑1 (MQ=14)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             >  1:2892595/1‑93 (MQ=14)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             <  1:179653/93‑1 (MQ=14)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCCCTGGGGTGCTATTGCTCGGCGccc                             >  1:6353279/1‑93 (MQ=11)
ggCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGccc                             <  2:5307092/93‑1 (MQ=31)
 gCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCGcccc                            >  1:199625/1‑93 (MQ=12)
 gCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGcccc                            <  1:6212233/93‑1 (MQ=17)
 gCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGcccc                            <  2:6955996/93‑1 (MQ=17)
 gCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGcccc                            <  2:4267501/93‑1 (MQ=17)
 gCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCGcccc                            >  1:6702421/1‑93 (MQ=17)
  cATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcg                           <  1:3602536/93‑1 (MQ=18)
  cATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcg                           >  2:2997519/1‑93 (MQ=18)
  cATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCGCCCcg                           >  1:1811150/1‑93 (MQ=18)
   aTCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCGcccccc                          >  2:4857573/1‑91 (MQ=255)
   aTCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcgc                          >  1:3944344/1‑93 (MQ=18)
   aTCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcgc                          <  2:148382/93‑1 (MQ=18)
   aTCCGTCTGTCTGCCGTGTGGCTGACGCAAGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCGcccccc                          >  2:5272580/1‑91 (MQ=255)
    tCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcgcg                         >  2:2223763/1‑93 (MQ=18)
    tCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcgcg                         >  1:407013/1‑93 (MQ=18)
    tCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcgcg                         >  1:6349455/1‑93 (MQ=33)
    tCCGTCTGTCTGCCGTGTGGATGACGCACGAACCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGAACGCTATGGCTGGGCACCCcgcg                         >  2:6959016/1‑93 (MQ=11)
     ccGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCcgccg                        >  2:1448983/1‑91 (MQ=16)
     ccGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCgcgg                        >  2:768674/1‑93 (MQ=18)
     ccGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCgcgg                        <  2:6131742/93‑1 (MQ=21)
     ccGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGcccccggg                        >  2:4122310/1‑89 (MQ=17)
      cGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCgcggc                       >  1:5895160/1‑93 (MQ=18)
       gtctgtctGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCgcggcg                      >  1:5865212/1‑93 (MQ=18)
       gtctgtctGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCgcggcg                      >  2:8615544/1‑93 (MQ=18)
        tctgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCgcggcgg                     >  2:2026349/1‑92 (MQ=18)
        tctgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGa                     >  2:772491/1‑93 (MQ=18)
        tctgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGa                     >  1:8381420/1‑93 (MQ=18)
         ctgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGaa                    <  1:5589019/93‑1 (MQ=21)
          tgtgtgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAAc                   >  1:8619425/1‑93 (MQ=255)
          tgtctgCCGTGTGGCTGACTCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAAc                   <  1:4149649/93‑1 (MQ=16)
          tgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAAc                   <  2:239535/93‑1 (MQ=21)
          tgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAAc                   <  2:8444449/93‑1 (MQ=21)
          tgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGAAGCCCCGCGGCGAAc                   >  1:2197660/1‑93 (MQ=255)
          tgtctgCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAAc                   >  1:1864671/1‑93 (MQ=33)
           gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACt                  <  2:6844106/93‑1 (MQ=21)
           gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACt                  <  1:6379684/93‑1 (MQ=21)
           gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACt                  <  2:4943371/93‑1 (MQ=21)
           gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACt                  >  1:386164/1‑93 (MQ=18)
           gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACt                  <  2:3253575/93‑1 (MQ=18)
           gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACCGCCCGCGGCGAACt                  >  1:544647/1‑93 (MQ=11)
            tctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGGGGCGAACTg                 >  1:4054479/1‑93 (MQ=16)
            tctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTg                 >  2:7409184/1‑93 (MQ=18)
             ctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTgg                <  1:1698134/93‑1 (MQ=21)
             ctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTgg                <  2:4193934/93‑1 (MQ=21)
             ctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTgg                <  2:179357/93‑1 (MQ=21)
              tgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggc               <  2:5100225/93‑1 (MQ=18)
              tgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCCGCGAACTggc               >  1:2637333/1‑93 (MQ=255)
               gctgtgtGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcg              <  2:7232815/90‑1 (MQ=16)
               gCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcg              <  2:3145446/93‑1 (MQ=21)
               gCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcg              <  2:2140261/93‑1 (MQ=21)
               gCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCCGCGAACTggcg              >  1:7909317/1‑93 (MQ=255)
               gCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGCCGCCCAGCGGCTCACCggcg              >  1:3594923/1‑93 (MQ=11)
                 cGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcggg            >  2:5683625/1‑92 (MQ=16)
                 cGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcggc            <  1:3489346/93‑1 (MQ=25)
                 cGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcggc            <  2:1734079/93‑1 (MQ=21)
                 cGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcggc            <  2:7546251/93‑1 (MQ=21)
                  gtgtGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTggcgggg           >  2:3222062/1‑91 (MQ=16)
                  gtgtGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGgcggcg           >  1:7205433/1‑93 (MQ=25)
                  gtgtGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGgcggcg           >  1:4101435/1‑93 (MQ=25)
                  gtgtGGCTGACGCGCGACCCGGAATACCCGGAGAATTTAACTGCCGCGTCGCTGGTGCGCTATGGCTGGACGCCCCGCACCGAAATGgcggcg           >  2:824767/1‑93 (MQ=11)
                   tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCTGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGgggcgg          >  2:6906166/1‑93 (MQ=11)
                   tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCGCCCCGCGGCGAACTGGgggcgg          >  1:6390913/1‑93 (MQ=11)
                   tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGCCGCCCCGCGGCGAACTGGcggcgg          >  1:5985633/1‑93 (MQ=255)
                   tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGcggcgg          <  2:1014972/93‑1 (MQ=255)
                   tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGcggcgg          >  1:7233825/1‑93 (MQ=255)
                   tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGcggcgg          >  1:4336387/1‑93 (MQ=255)
                   tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGCGCGCTATGGCTGGAAGCCCCGCGGCGAACTGGcggcgg          >  2:6906164/1‑93 (MQ=255)
                    gtgtCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGgt         <  2:8407187/93‑1 (MQ=255)
                    gtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGGCCCCTGGGGCCCTTTGGCGGGACGCCCCCCCGCGAACCgggggggg          >  2:6315717/1‑93 (MQ=2)
                    gtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGgt         <  1:1133106/93‑1 (MQ=255)
                     tgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGgtg        >  1:6772561/1‑93 (MQ=37)
                      ggCTGACGCGCGACCCGGAATCCCCGGAGAATTTCCCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGgtgt       <  1:3965454/93‑1 (MQ=255)
                      ggCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGAGACCCCCGGCGAACCGGgggcggggg       >  1:4845389/1‑90 (MQ=11)
                      ggCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGgtgt       <  1:3285826/93‑1 (MQ=255)
                       gCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTa      >  1:971920/1‑93 (MQ=255)
                       gCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTa      <  2:3224987/93‑1 (MQ=255)
                        cTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTAt     <  2:4527609/93‑1 (MQ=255)
                        cTGACGCGCGACCAGGAATACCCGGAGAAGTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTAt     <  1:8642093/93‑1 (MQ=255)
                         tGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATg    <  1:6919786/93‑1 (MQ=255)
                          gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCGCCCCGCGGCGAACTGGCGGCGGTGTATGa   >  1:6202191/1‑93 (MQ=255)
                          gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGa   <  1:7972144/93‑1 (MQ=37)
                          gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGa   >  1:2261750/1‑93 (MQ=255)
                          gACGCGCGACCCGGAATACCCGGAGAATTTACCTACCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGa   >  2:5497646/1‑93 (MQ=255)
                           gcgcgcgACCCGGAATACCCGGAGAATTTCCCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAc  <  1:3885556/92‑1 (MQ=255)
                           aCGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAc  <  2:47633/93‑1 (MQ=255)
                           aCGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAc  <  2:877551/93‑1 (MQ=255)
                           aCGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCCGCGAACCGGGGGCGGTGGATGAc  >  2:5020115/1‑93 (MQ=11)
                           aCGCGCGACCCGGAATACCCGGAGAATGTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAc  <  1:6293941/93‑1 (MQ=255)
                               |                                                                                         
GGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGAC  >  NC_000913/729667‑729786

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: