Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,510,121 +G intergenic (‑358/+151) yhfZ ← / ← trpS hypothetical protein/tryptophan tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,510,1201.G95.8% 798.6 / 25.1 263intergenic (‑357/+152)yhfZ/trpShypothetical protein/tryptophan tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (8/3);  new base G (144/108);  total (152/111)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 296 total aligned reads displayed.

AAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGG  >  NC_000913/3510028‑3510209
                                                                                              |                                                                                         
aaGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTcccc                                                                                            <  2:4363769/93‑1 (MQ=255)
    tAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtt                                                                                       >  1:6296433/1‑93 (MQ=255)
     aaGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgttt                                                                                      <  2:2679227/93‑2 (MQ=255)
      aGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtttc                                                                                     >  2:6037377/1‑91 (MQ=255)
       ggCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCa                                                                                     >  2:5205477/1‑93 (MQ=255)
       ggCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCa                                                                                     >  1:7228159/1‑93 (MQ=255)
       ggCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCa                                                                                     <  1:4957361/93‑1 (MQ=255)
        gcgcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCag                                                                                    >  2:3011101/1‑93 (MQ=255)
        gcgcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCag                                                                                    <  2:1580130/93‑1 (MQ=255)
         cgcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCaga                                                                                   <  1:1415666/93‑1 (MQ=255)
         cgcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCaga                                                                                   <  2:6843618/93‑1 (MQ=255)
          gcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCagag                                                                                  >  2:857449/1‑93 (MQ=255)
          gcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCagag                                                                                  >  2:1412889/1‑93 (MQ=255)
             aCAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtga                                                                               >  1:3267339/1‑93 (MQ=255)
             aCAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtga                                                                               >  1:3850130/1‑93 (MQ=255)
                gCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGTATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCCGCGcgatga                                                                            >  1:6601527/1‑93 (MQ=255)
                  gagaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATa                                                                          >  1:1512781/1‑93 (MQ=255)
                   agaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaa                                                                         >  2:3230698/1‑93 (MQ=255)
                   agaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaa                                                                         >  1:783159/1‑93 (MQ=255)
                    gaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaa                                                                        >  1:4150149/1‑93 (MQ=255)
                      tgtgGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAg                                                                      >  1:6701078/1‑93 (MQ=255)
                      tgtgGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACAGGGTTTCCCCGGTTTCAGAGTGATGATAAAAg                                                                      >  2:3085460/1‑93 (MQ=255)
                          gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaa                                                                  >  1:3186768/1‑93 (MQ=255)
                          gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaa                                                                  <  1:2299322/93‑1 (MQ=255)
                           aaaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaaa                                                                 <  1:3323754/93‑1 (MQ=255)
                             aCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAtt                                                               >  2:3950437/1‑93 (MQ=255)
                              cGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTg                                                              >  1:768229/1‑93 (MQ=255)
                               ggCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGc                                                             <  1:7072091/93‑1 (MQ=255)
                                gCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGcc                                                            >  1:5644396/1‑93 (MQ=255)
                                  gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTg                                                          >  2:7434105/1‑93 (MQ=255)
                                  gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTg                                                          >  2:5920742/1‑93 (MQ=255)
                                  gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                                         >  1:2723204/1‑93 (MQ=255)
                                   aGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                                         <  2:7300689/93‑1 (MQ=255)
                                   aGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                                         <  1:2375684/93‑1 (MQ=255)
                                    ggCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                                        <  1:2028785/93‑1 (MQ=255)
                                    ggCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                                        >  1:1762131/1‑93 (MQ=255)
                                    ggCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                                        <  1:3531100/93‑1 (MQ=255)
                                      cACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgc                                                      <  1:6604970/93‑1 (MQ=255)
                                       aCTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcg                                                     >  1:4686812/1‑93 (MQ=255)
                                       aCTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGCGATGATAAAAGCAAAATTGCCTGATgcg                                                     >  2:3399901/1‑93 (MQ=255)
                                        ttttacaTTTTTTCGGATTATTTCAAATGAGGATAAAAAAACCGGGTTTCTCCGGTTTCAGAGTGATGATAAAAGCAAAAGTTCCTGATgcgc                                                    >  2:7018331/5‑93 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                    >  1:7350769/1‑93 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                    >  1:2898674/1‑93 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                    <  1:404588/93‑1 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                    >  2:7569619/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                   >  1:6163099/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                   >  1:3795181/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                   >  1:4652340/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                   <  1:4387347/93‑1 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                   >  2:426814/1‑93 (MQ=255)
                                          tCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTa                                                  >  1:4760323/1‑93 (MQ=255)
                                             caTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGc                                               <  1:1272392/93‑1 (MQ=255)
                                             caTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGc                                               <  1:2234545/93‑1 (MQ=255)
                                             caTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGc                                               <  1:6463856/93‑1 (MQ=255)
                                              aTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCt                                              >  1:4288150/1‑93 (MQ=255)
                                               ttttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCtt                                             >  2:4093038/1‑93 (MQ=255)
                                                tttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTa                                            >  1:6006209/1‑93 (MQ=255)
                                                tttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTa                                            >  2:2360891/1‑93 (MQ=255)
                                                  tttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATc                                          <  1:2555961/93‑1 (MQ=255)
                                                  tttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGCGATGATAAAAGCAAAATTGCCTTATGCGCTACACTTATc                                          >  2:5272827/1‑93 (MQ=38)
                                                   ttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCa                                         <  1:7584857/93‑1 (MQ=255)
                                                    tCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAg                                        <  1:5866821/93‑1 (MQ=255)
                                                      ggATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                      >  1:6182603/1‑93 (MQ=255)
                                                      ggATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                      >  2:6200252/1‑93 (MQ=255)
                                                      ggATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                      >  1:3469612/1‑93 (MQ=255)
                                                        attattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCt                                    <  2:88759/93‑1 (MQ=255)
                                                        attattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCt                                    >  2:7455173/1‑93 (MQ=255)
                                                        attattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCt                                    >  1:6866026/1‑93 (MQ=255)
                                                          tattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                  >  1:1060632/1‑93 (MQ=255)
                                                           attGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                 >  1:5505407/1‑93 (MQ=255)
                                                           attGCAAATGAGAATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                 <  2:6499744/93‑1 (MQ=255)
                                                            ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAt                                >  2:7097479/1‑93 (MQ=255)
                                                            ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAt                                <  1:1095704/93‑1 (MQ=255)
                                                             tGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt                               <  2:2191935/93‑1 (MQ=255)
                                                              gCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAttt                              <  2:4503121/93‑1 (MQ=255)
                                                                aaaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTcc                            <  2:6761314/93‑1 (MQ=255)
                                                                 aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCt                           <  2:1094929/93‑1 (MQ=255)
                                                                  aTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCtt                          <  2:3685969/93‑1 (MQ=255)
                                                                  aTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCtt                          <  2:3526613/93‑1 (MQ=255)
                                                                   tGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTg                         >  1:1978574/1‑93 (MQ=255)
                                                                    gAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGc                        >  1:5554768/1‑93 (MQ=255)
                                                                         tAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAatat                   >  2:2895780/1‑93 (MQ=255)
                                                                              aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCa              >  2:4893092/1‑93 (MQ=255)
                                                                              aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCa              >  1:5032261/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt             >  2:5854990/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt             >  2:1776507/1‑93 (MQ=255)
                                                                                aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtt            >  2:2466858/1‑93 (MQ=255)
                                                                                 ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTa           >  2:4659110/1‑93 (MQ=255)
                                                                                 ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTa           >  1:5114369/1‑93 (MQ=255)
                                                                                  cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTAc          <  1:963598/93‑1 (MQ=255)
                                                                                  cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTAc          <  1:3625756/93‑1 (MQ=255)
                                                                                  cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTAc          <  1:5465482/93‑1 (MQ=255)
                                                                                    ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACtt        <  2:6012856/93‑1 (MQ=255)
                                                                                     gTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACttt       <  2:7609209/93‑1 (MQ=255)
                                                                                      tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTg      >  2:570472/1‑93 (MQ=255)
                                                                                      tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTg      <  1:902971/93‑1 (MQ=255)
                                                                                      ttccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGt     <  2:5112770/90‑1 (MQ=255)
                                                                                       tccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTa    <  1:3119859/91‑1 (MQ=255)
                                                                                        ccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAg   <  2:6671311/92‑1 (MQ=255)
                                                                                         cccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAgg  >  2:3338347/1‑93 (MQ=255)
                                                                                              |                                                                                         
AAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGG  >  NC_000913/3510028‑3510209

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: