Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,250 G→C 20.2% G194G (GGG→GGC dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2500GC20.2% 419.9 / 104.2 183G194G (GGG→GGCdusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (70/76);  new base C (22/15);  total (92/91)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

Only 100 of 197 total aligned reads displayed.

GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGT  >  NC_000913/4262151‑4262349
                                                                                                   |                                                                                                   
ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                                                                     <  1:3359423/100‑1 (MQ=255)
 gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGt                                                                                                    <  1:2107022/100‑1 (MQ=255)
  aCAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTACGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGt                                                                                                    <  2:2254018/99‑1 (MQ=255)
   cAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTa                                                                                                  <  2:1722253/100‑1 (MQ=255)
     gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                <  1:903796/100‑1 (MQ=255)
     gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                >  1:3724959/1‑100 (MQ=255)
     gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGGGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                >  2:710010/1‑100 (MQ=255)
      cTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGc                                                                                               >  1:1240297/1‑100 (MQ=255)
       tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                              <  1:2503181/100‑1 (MQ=255)
       tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                              >  1:3192212/1‑100 (MQ=255)
       tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                              <  1:1589834/100‑1 (MQ=255)
            aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaa                                                                                          <  1:932157/99‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                        <  1:3793818/100‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                        <  1:1409549/100‑1 (MQ=255)
                ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcca                                                                                     <  2:3050657/100‑4 (MQ=255)
                ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACACACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaaa                                                                                     >  1:3324726/1‑100 (MQ=255)
                 tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                    <  2:1553644/100‑5 (MQ=255)
                   tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAcc                                                                                  >  1:1037810/1‑100 (MQ=255)
                   tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAcc                                                                                  >  2:1846100/1‑100 (MQ=255)
                   tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggc                                                                                  <  2:200312/100‑7 (MQ=38)
                    gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCg                                                                                 <  1:3533619/100‑1 (MQ=255)
                    gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggct                                                                                 <  1:1458221/100‑8 (MQ=38)
                       aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                              <  1:392806/100‑1 (MQ=255)
                        tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAACCCGTGa                                                                              <  2:1389983/99‑1 (MQ=255)
                        tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGag                                                                             <  2:58798/100‑2 (MQ=255)
                             tcaATACCGTTTCCGGCAAATGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                        >  2:2339051/1‑100 (MQ=255)
                              caATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccg                                                                       >  1:990752/1‑100 (MQ=255)
                                   ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                  >  1:1669318/1‑100 (MQ=255)
                                   ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                  >  1:299218/1‑100 (MQ=255)
                                    cGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATACACCCTCGTTCAGCCTGGCTTTCGGGGTACAGCCCGAAAGAAAACAGTGAAAACacgcag                                                                    >  2:3371379/1‑95 (MQ=255)
                                      tttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGa                                                               <  1:434742/100‑1 (MQ=255)
                                      ttgccGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGa                                                               <  2:2692021/97‑1 (MQ=255)
                                       ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                              <  1:110413/100‑1 (MQ=255)
                                       ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                              <  1:2749688/100‑1 (MQ=255)
                                        tCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGACATCCCGCCGCTCGAtt                                                             >  2:2231347/1‑100 (MQ=255)
                                        tCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAtt                                                             >  1:1282494/1‑100 (MQ=255)
                                        tCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAtt                                                             <  1:177118/100‑1 (MQ=255)
                                          cGGCAAAGGCGCGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                           >  2:762240/1‑100 (MQ=255)
                                          cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGACAGAAAACCGTGAAATCCCGCCGATCGCTTAt                                                           >  2:847764/1‑100 (MQ=255)
                                          cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGCTTAt                                                           >  1:3471/1‑100 (MQ=255)
                                           ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                          >  2:3738351/1‑100 (MQ=255)
                                           ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCAAGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                          >  2:2586193/1‑100 (MQ=255)
                                               aaGGCGGGTGTGAGAGGTTCATCATCCACGCATGTAAATCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgg                                                       <  2:3510314/99‑2 (MQ=255)
                                               aaGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgcg                                                      >  2:964229/1‑100 (MQ=255)
                                                 ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                    <  1:975385/100‑1 (MQ=255)
                                                 ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCGTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                    <  2:2446286/100‑1 (MQ=255)
                                                   cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGGCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATACGCgtgtg                                                  >  2:2513095/1‑100 (MQ=255)
                                                     aGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCtcgtgt                                                 >  2:2227005/1‑99 (MQ=255)
                                                       tgtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTa                                                <  2:2814805/98‑1 (MQ=255)
                                                       tgtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                              <  1:657833/100‑1 (MQ=255)
                                                       tgtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                              <  1:3367450/100‑1 (MQ=255)
                                                       tgtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                              <  1:463941/100‑1 (MQ=255)
                                                        gtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa                                             >  1:2962480/1‑100 (MQ=255)
                                                         tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCaa                                            >  1:1964020/1‑100 (MQ=255)
                                                          gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgggagtgaac                                           >  2:3005693/1‑91 (MQ=255)
                                                           agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                          <  2:2435571/100‑1 (MQ=255)
                                                           agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                          >  1:1941093/1‑100 (MQ=255)
                                                              tGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGaa                                       <  1:2913080/100‑1 (MQ=255)
                                                              tGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGa                                        <  1:3625272/99‑1 (MQ=255)
                                                              tGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGAGTGTATCAACTGaa                                       >  1:3658086/1‑100 (MQ=255)
                                                               gTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAg                                      <  1:2731261/100‑1 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                    >  1:1843516/1‑100 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                    >  2:3756105/1‑100 (MQ=255)
                                                                    tcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                 <  1:102672/100‑1 (MQ=255)
                                                                         cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                            >  2:1279346/1‑100 (MQ=255)
                                                                         cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                            >  1:2766388/1‑100 (MQ=255)
                                                                          acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTcc                           <  2:1922609/100‑1 (MQ=255)
                                                                            gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGc                         >  1:1911686/1‑100 (MQ=255)
                                                                             cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa                        >  2:329795/1‑100 (MQ=255)
                                                                               cgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                      <  1:1358143/100‑1 (MQ=255)
                                                                                gTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt                     >  2:463903/1‑100 (MQ=255)
                                                                                 tcacggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                    >  2:3106946/11‑100 (MQ=37)
                                                                                   aaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGAc                  <  1:1783878/100‑1 (MQ=255)
                                                                                    cggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                 >  1:689080/8‑100 (MQ=38)
                                                                                    cggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                 >  1:2803924/8‑100 (MQ=38)
                                                                                     gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGCTTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa                <  2:972462/100‑1 (MQ=255)
                                                                                     gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa                >  2:1969712/1‑100 (MQ=255)
                                                                                        tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt             >  1:2219743/1‑100 (MQ=255)
                                                                                         ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc            >  1:789323/3‑100 (MQ=255)
                                                                                         ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc            >  1:1502842/3‑100 (MQ=255)
                                                                                         ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc            >  1:55527/3‑100 (MQ=255)
                                                                                         ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc            >  2:2914339/1‑100 (MQ=255)
                                                                                         ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc            >  2:142760/1‑100 (MQ=255)
                                                                                          tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa                >  2:2564599/2‑95 (MQ=255)
                                                                                           cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          <  2:2160669/100‑1 (MQ=255)
                                                                                           cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          >  1:1059475/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAAGCCCGAAAGAAAACCGTGAAATCCCACCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg           >  2:2356160/1‑99 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          >  2:3639718/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          >  1:2650440/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          >  2:2745694/1‑100 (MQ=255)
                                                                                            tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt         <  2:162667/100‑1 (MQ=255)
                                                                                            tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGGGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt         >  2:517080/1‑100 (MQ=255)
                                                                                            tttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt         >  1:1736599/1‑100 (MQ=255)
                                                                                             ttCGGGGCTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCGACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAtn        <  2:2046250/100‑2 (MQ=255)
                                                                                              tCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTa       <  2:1037801/100‑1 (MQ=255)
                                                                                               cGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTaa      <  2:1240288/100‑1 (MQ=255)
                                                                                                ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAAc     >  2:1467678/1‑100 (MQ=255)
                                                                                                 gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg    >  2:949597/1‑100 (MQ=255)
                                                                                                  ggTTAAGCCCGAAAGAAAACCGTGAACTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg   >  1:2320612/1‑100 (MQ=255)
                                                                                                  ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg    <  2:3241266/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGT  >  NC_000913/4262151‑4262349

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