Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,339,304 +C coding (398/987 nt) yrbH → putative isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,339,3011.C98.4% 747.3 / 3.5 242coding (395/987 nt)yrbHputative isomerase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base C (138/100);  total (139/103)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 264 total aligned reads displayed.

GCTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑CCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGAC  >  NC_000913/3339210‑3339388
                                                                                            |                                                                                        
gCTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c                                                                                         <  2:269854/93‑1 (MQ=255)
  tATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccc                                                                                       >  1:4127173/1‑93 (MQ=255)
  tATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccc                                                                                       >  1:4286164/1‑93 (MQ=255)
  tATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccc                                                                                       <  2:646502/93‑1 (MQ=255)
   aTCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑cccc                                                                                      >  2:1451707/1‑92 (MQ=255)
    tctcTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccccg                                                                                     >  1:2601786/1‑91 (MQ=255)
    tctcTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccccg                                                                                     >  1:1357160/1‑91 (MQ=255)
     ctctAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGTATCACCGGTCG‑ccccgg                                                                                    >  1:3184319/1‑90 (MQ=255)
      tctAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccccgga                                                                                   >  2:2663725/1‑89 (MQ=255)
        tAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑Ggaga                                                                                 >  2:1222692/1‑93 (MQ=255)
          aCTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTGTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGagc                                                                               >  1:2349363/1‑93 (MQ=255)
          aCTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGagc                                                                               >  1:4649374/1‑93 (MQ=255)
          aCTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGagc                                                                               <  1:3090540/93‑1 (MQ=255)
          aCTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGagc                                                                               >  1:2853029/1‑93 (MQ=255)
            tctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGagcag                                                                             >  1:3510229/1‑93 (MQ=255)
             ctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGagcagc                                                                            <  2:2248977/93‑1 (MQ=255)
              tGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAgcagca                                                                           <  1:3939541/93‑1 (MQ=255)
              tGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAgcagca                                                                           >  1:2172611/1‑93 (MQ=255)
               ggTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTCCCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCAt                                                                          <  2:1123299/93‑1 (MQ=255)
               ggTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCAt                                                                          <  1:4020333/93‑1 (MQ=255)
               ggTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCAt                                                                          >  2:4221328/1‑93 (MQ=255)
                   aaTCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑CCCTGGAGAGCAGCATGgcg                                                                      >  1:1753176/1‑93 (MQ=255)
                     tCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGgcgcg                                                                    >  2:1831525/1‑93 (MQ=255)
                     tCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGgcgcg                                                                    >  1:1938192/1‑93 (MQ=255)
                       cAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGgcgcgcg                                                                  >  2:780897/1‑93 (MQ=255)
                       cAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGgcgcgcg                                                                  <  1:225404/93‑1 (MQ=255)
                         gCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGcgcc                                                                >  1:2064116/1‑93 (MQ=255)
                          cGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGcgccg                                                               <  2:4395672/93‑1 (MQ=255)
                           gAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGcgccgt                                                              >  2:1507021/1‑92 (MQ=255)
                           gAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGcgccgc                                                              >  1:2298750/1‑93 (MQ=255)
                            aaaTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCa                                                             >  1:2050338/1‑93 (MQ=255)
                            aaaTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCa                                                             >  2:3151848/1‑93 (MQ=255)
                            aaaTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCa                                                             <  2:1940587/93‑1 (MQ=255)
                              aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGa                                                           >  2:73133/1‑93 (MQ=255)
                              aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGa                                                           >  2:3892546/1‑93 (MQ=255)
                               tCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGAt                                                          >  2:4180818/1‑93 (MQ=255)
                                 aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        >  2:2302334/1‑93 (MQ=255)
                                 aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        >  1:4063673/1‑93 (MQ=255)
                                 aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        <  1:1017736/93‑1 (MQ=255)
                                 aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        <  1:16906/93‑1 (MQ=255)
                                 aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        <  1:2189702/93‑1 (MQ=255)
                                 aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        <  2:3472975/93‑1 (MQ=255)
                                   ggCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGc                                                      <  2:288146/93‑1 (MQ=255)
                                      cTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATc                                                   <  1:23806/93‑1 (MQ=255)
                                      cTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATc                                                   >  2:334291/1‑93 (MQ=255)
                                      cTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATc                                                   <  2:2482618/93‑1 (MQ=255)
                                      cTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATc                                                   >  2:353053/1‑93 (MQ=255)
                                      cTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATc                                                   <  1:954368/93‑1 (MQ=255)
                                           ttCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgtgt                                              >  2:4280822/1‑93 (MQ=255)
                                            tCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgtgtg                                             <  1:3407138/93‑1 (MQ=255)
                                              cAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGtt                                           >  1:3843900/1‑93 (MQ=255)
                                              cAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGtt                                           >  1:2690024/1‑93 (MQ=255)
                                              cAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGtt                                           >  1:3369823/1‑93 (MQ=255)
                                                gggCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaa                                         >  2:2028450/3‑93 (MQ=255)
                                                gTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaa                                         >  1:4237898/1‑93 (MQ=255)
                                                gTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaa                                         <  2:1358737/93‑1 (MQ=255)
                                                 tGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaaa                                        <  1:3820296/93‑1 (MQ=255)
                                                 tGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaaa                                        <  2:1263729/93‑1 (MQ=255)
                                                 tGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaaa                                        >  2:2849581/1‑93 (MQ=255)
                                                 tGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGCCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaaa                                        >  1:1056205/1‑93 (MQ=255)
                                                   cTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGt                                      <  2:1771896/93‑1 (MQ=255)
                                                    ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                     <  1:4387939/93‑1 (MQ=255)
                                                    ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                     >  2:4543686/1‑93 (MQ=255)
                                                    ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                     >  2:3736532/1‑93 (MQ=255)
                                                     tAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAg                                    <  1:1436408/93‑1 (MQ=255)
                                                     tAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAg                                    <  2:4032332/93‑1 (MQ=255)
                                                      aaGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGc                                   >  2:2118126/1‑93 (MQ=255)
                                                        gCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGa                                 >  2:325751/1‑93 (MQ=255)
                                                        gCGTCTTCACGTACCGTTAATATGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGa                                 >  2:2984591/1‑93 (MQ=255)
                                                          gTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaaa                               >  2:3117207/1‑93 (MQ=255)
                                                          gTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaaa                               <  1:793897/93‑1 (MQ=255)
                                                          gTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaaa                               <  1:2763486/93‑1 (MQ=255)
                                                           tCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGCTGTGCATCTGTGTGTTAAAGTAGCGaaaa                              >  1:3535745/1‑92 (MQ=255)
                                                           tCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaag                              <  2:504478/93‑1 (MQ=255)
                                                           tCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaag                              >  1:1241680/1‑93 (MQ=255)
                                                             ttCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaa                            >  1:2199919/1‑93 (MQ=255)
                                                             ttCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaa                            <  1:5530/93‑1 (MQ=255)
                                                             ttCACGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaa                            >  2:1110246/1‑93 (MQ=255)
                                                                aCGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGcc                         >  1:3337777/1‑93 (MQ=255)
                                                                aCGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGcc                         <  2:1550517/93‑1 (MQ=255)
                                                                 cGTACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCt                        >  2:889672/1‑93 (MQ=255)
                                                                   tACCGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGt                      <  2:164979/93‑1 (MQ=255)
                                                                      cGTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCg                   >  1:4031264/1‑93 (MQ=255)
                                                                       gTTAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGt                  >  2:2535963/1‑93 (MQ=255)
                                                                         tAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTa                <  1:3201690/93‑1 (MQ=255)
                                                                         tAATCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTa                <  2:3639881/93‑1 (MQ=255)
                                                                          aaTCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAg               <  2:218390/93‑1 (MQ=255)
                                                                          aaTCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAg               <  2:240488/93‑1 (MQ=255)
                                                                          aaTCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAg               >  1:848536/1‑93 (MQ=255)
                                                                           aTCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAgg              >  1:759465/1‑93 (MQ=255)
                                                                            tCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTGCGTTAggg             >  1:1600179/1‑93 (MQ=255)
                                                                            tCTGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAggg             >  1:3059910/1‑93 (MQ=255)
                                                                              tGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCt           <  1:2464035/93‑1 (MQ=255)
                                                                              tGCATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCt           >  2:1180822/1‑93 (MQ=255)
                                                                               gCATCACCGGCCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTg          <  2:1720971/93‑1 (MQ=255)
                                                                                cATCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTgg         <  1:1461295/93‑1 (MQ=255)
                                                                                  tCACCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCa       >  1:128319/1‑93 (MQ=255)
                                                                                    aCCGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCAcc     >  2:3773123/1‑93 (MQ=255)
                                                                                     ccGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCg    >  1:98984/1‑93 (MQ=255)
                                                                                      cGGTCGCCCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGa   <  1:3588949/93‑1 (MQ=255)
                                                                                            |                                                                                        
GCTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑CCC‑GGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGAC  >  NC_000913/3339210‑3339388

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