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.C99.6% 762.5 / ‑4.1 238coding (395/987 nt)yrbHputative isomerase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (128/109);  total (128/110)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 263 total aligned reads displayed.

CTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑C‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGA  >  NC_000913/3339211‑3339387
                                                                                           |                                                                                       
 tATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c‑cc                                                                                     >  1:151414/1‑93 (MQ=255)
 tATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c‑cc                                                                                     >  2:4854291/1‑93 (MQ=255)
  aTCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c‑ccc                                                                                    <  1:860653/93‑2 (MQ=255)
   tctcTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c‑cccg                                                                                   >  2:2202494/1‑91 (MQ=255)
    ctctAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c‑cccgg                                                                                  <  1:2657315/93‑4 (MQ=255)
       tAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGgaga                                                                                <  2:5955504/93‑1 (MQ=255)
       tAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGgaga                                                                                <  1:4263334/93‑1 (MQ=255)
         aCTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGagc                                                                              >  2:123555/1‑93 (MQ=255)
          ctctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGagca                                                                             >  1:1011094/1‑93 (MQ=255)
          ctctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGagca                                                                             <  1:335962/93‑1 (MQ=255)
           tctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGagcag                                                                            >  2:4959795/1‑93 (MQ=255)
             tGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAgcagca                                                                          <  1:2558285/93‑1 (MQ=255)
             tGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAgcagca                                                                          <  1:197186/93‑1 (MQ=255)
              ggTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCAt                                                                         <  1:1630442/93‑1 (MQ=255)
              ggTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCAt                                                                         <  2:3932621/93‑1 (MQ=255)
              ggTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCAt                                                                         >  2:5293450/1‑93 (MQ=255)
               gTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATg                                                                        >  1:2894496/1‑93 (MQ=255)
               gTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATg                                                                        >  1:4032059/1‑93 (MQ=255)
                tGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATgg                                                                       >  2:1781118/1‑93 (MQ=255)
                tGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATgg                                                                       <  2:2131040/93‑1 (MQ=255)
                   aTCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGgcgc                                                                    >  2:5364001/1‑93 (MQ=255)
                     ccAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGgcgcgc                                                                  <  2:1499511/93‑1 (MQ=255)
                     ccAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGgcgcgc                                                                  >  1:5617631/1‑93 (MQ=255)
                       aGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGgcgcgcgc                                                                <  1:643077/93‑1 (MQ=255)
                        gCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgcc                                                               >  1:5292461/1‑93 (MQ=255)
                        gCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgcc                                                               >  2:4282808/1‑93 (MQ=255)
                        gCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgcc                                                               >  1:4371376/1‑93 (MQ=255)
                        gCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgcc                                                               >  1:3326508/1‑93 (MQ=255)
                        gCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGGACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgcc                                                               >  1:4433636/1‑93 (MQ=255)
                        gCGAAACCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgcc                                                               <  2:6028486/93‑1 (MQ=255)
                          gAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgccgc                                                             <  2:4937392/93‑1 (MQ=255)
                          gAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGcgccgc                                                             >  1:1813803/1‑93 (MQ=255)
                           aaaTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCa                                                            >  2:309039/1‑93 (MQ=255)
                           aaaTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCa                                                            >  2:2051570/1‑93 (MQ=255)
                            aaTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAg                                                           >  2:3478668/1‑93 (MQ=255)
                             aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGa                                                          <  2:5545795/93‑1 (MQ=255)
                             aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGa                                                          <  2:5721314/93‑1 (MQ=255)
                             aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGa                                                          <  1:5105569/93‑1 (MQ=255)
                              tCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGAt                                                         <  1:4217004/93‑1 (MQ=255)
                              tCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGAt                                                         >  1:4732299/1‑93 (MQ=255)
                                aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                       <  2:1332587/93‑1 (MQ=255)
                                 cGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGAtgtg                                                      <  2:1640061/93‑1 (MQ=255)
                                  ggCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGGCGCAGATGTGc                                                     >  2:3227386/1‑93 (MQ=255)
                                  ggCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGc                                                     >  1:6075044/1‑93 (MQ=255)
                                   gCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCa                                                    >  2:5453972/1‑93 (MQ=255)
                                   gCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCa                                                    >  1:5715949/1‑93 (MQ=255)
                                   gCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCa                                                    >  1:5252816/1‑93 (MQ=255)
                                   gCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCa                                                    <  1:38686/93‑1 (MQ=255)
                                    ccTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCAt                                                   >  1:547371/1‑93 (MQ=255)
                                    ccTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCAt                                                   >  2:737475/1‑93 (MQ=255)
                                        aaTTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgt                                               <  1:4190176/93‑1 (MQ=255)
                                          ttCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgtgt                                             <  1:1310190/93‑1 (MQ=255)
                                           tCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgtgtg                                            >  1:1762745/1‑93 (MQ=255)
                                            ccAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTgtgtgt                                           >  2:5145751/1‑93 (MQ=255)
                                             cAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGtt                                          >  1:3282247/1‑93 (MQ=255)
                                             cAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGtt                                          >  1:5416641/1‑93 (MQ=255)
                                              aGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTa                                         >  1:619165/1‑93 (MQ=255)
                                              aGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTa                                         >  2:4717990/1‑93 (MQ=255)
                                               gTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaa                                        >  1:507930/1‑93 (MQ=255)
                                               gTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaa                                        >  1:2273611/1‑93 (MQ=255)
                                                 gCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAg                                      <  1:80854/93‑1 (MQ=255)
                                                   ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                    <  2:5036225/93‑1 (MQ=255)
                                                   ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                    >  2:3831897/1‑93 (MQ=255)
                                                   ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                    <  1:5630853/93‑1 (MQ=255)
                                                   gtacgCGTCTTCACGTACAGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                    <  1:74293/89‑1 (MQ=255)
                                                    tAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAg                                   <  2:1120505/93‑1 (MQ=255)
                                                      aGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCg                                 <  1:2254076/93‑1 (MQ=255)
                                                      aGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCg                                 <  1:5923254/93‑1 (MQ=255)
                                                        cGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaa                               <  2:3988265/93‑1 (MQ=255)
                                                        cGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaa                               <  2:5196975/93‑1 (MQ=255)
                                                         gTCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaaa                              <  1:6031931/93‑1 (MQ=255)
                                                          tCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaag                             >  2:2845598/1‑93 (MQ=255)
                                                          tCTTCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaag                             <  1:4173167/93‑1 (MQ=255)
                                                            ttCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaa                           >  2:770893/1‑93 (MQ=255)
                                                            ttCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaa                           >  1:5845271/1‑93 (MQ=255)
                                                             tCACGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaag                          <  1:2608669/93‑1 (MQ=255)
                                                               ccGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGcc                        <  2:3034338/92‑1 (MQ=255)
                                                               aCGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCGTGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGcc                        >  2:2792787/1‑93 (MQ=255)
                                                               aCGTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaaggc                        >  1:2321151/1‑91 (MQ=255)
                                                                 gTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTg                      >  2:4356295/1‑93 (MQ=255)
                                                                 gTACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTg                      <  2:1613274/93‑1 (MQ=255)
                                                                  tACCGTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGt                     >  1:5186197/1‑93 (MQ=255)
                                                                      gTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGt                 <  2:1975136/93‑1 (MQ=255)
                                                                      gTTAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGt                 >  2:696090/1‑93 (MQ=255)
                                                                        tAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTa               <  2:4915283/93‑1 (MQ=255)
                                                                        tAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTa               >  1:3586997/1‑93 (MQ=255)
                                                                        tAATCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTa               <  1:3666004/93‑1 (MQ=255)
                                                                           tCTGCATCTCCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAggg            >  2:6147197/1‑93 (MQ=255)
                                                                           tCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAggg            <  1:4808985/93‑1 (MQ=255)
                                                                           tCTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAggg            <  2:4388501/93‑1 (MQ=255)
                                                                            cTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGc           <  2:239909/93‑1 (MQ=255)
                                                                            cTGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGc           <  1:913786/93‑1 (MQ=255)
                                                                            cTGCATCACCGGTCG‑CACCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGc           <  1:3082129/93‑1 (MQ=255)
                                                                             tGCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCt          <  1:3706080/93‑1 (MQ=255)
                                                                              gCATCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTg         >  1:2373059/1‑93 (MQ=255)
                                                                                aTCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGc       >  2:4565497/1‑93 (MQ=255)
                                                                                aTCACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGc       >  1:99073/1‑93 (MQ=255)
                                                                                  cACCGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCAc     >  2:3422295/1‑93 (MQ=255)
                                                                                    ccGGTCGCC‑CCGGAGAGCAGCATGGCTCTCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCg   <  1:1041214/93‑1 (MQ=255)
                                                                                     cGGTCGCC‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGa  >  1:6225655/1‑93 (MQ=255)
                                                                                           |                                                                                       
CTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑C‑CCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGA  >  NC_000913/3339211‑3339387

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