Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,319 Δ1 bp 16.9% coding (855/897 nt) ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,3190C.16.9% 569.6 / 117.5 195coding (855/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (77/85);  new base . (20/13);  total (97/98)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

Only 100 of 199 total aligned reads displayed.

TCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTC  >  NC_000913/3066220‑3066418
                                                                                                     |                                                                                                   
    aGCCACGTCTTCGACTTACACCGTATTTTAATTTCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGC‑ATATTGCTGAAG‑ATCGt                                                                                                <  1:1335625/100‑1 (MQ=38)
     gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGt                                                                                                >  2:1361463/1‑100 (MQ=255)
      ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCg                                                                                                 <  1:3663265/98‑1 (MQ=255)
      ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTg                                                                                               <  1:4226906/100‑1 (MQ=255)
        aCGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGc                                                                                             <  2:2087869/100‑1 (MQ=255)
         cGTCTTCGACTTACACCGTAATTTAAGTGCGACCGGCTCGAGGCGGGGTGAATGACTCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACg                                                                                            <  2:673647/100‑1 (MQ=12)
          gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                           <  1:1071880/100‑1 (MQ=255)
          gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgg                                                                                           <  2:1337173/100‑1 (MQ=21)
           tCTTCGACTTACACCGTAATTTAATTTCACCCGGCTCGATGCTGTGTTAATGACTCGGCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACg                                                                                            <  1:2951262/98‑1 (MQ=12)
           tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATTCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACg                                                                                            <  1:1657917/98‑1 (MQ=18)
           tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGt                                                                                          <  1:942718/100‑1 (MQ=255)
           tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACggg                                                                                          <  1:2167009/100‑2 (MQ=18)
            cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         <  2:3758605/100‑1 (MQ=255)
            cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         >  2:1132253/1‑100 (MQ=255)
             ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGAGGtat                                                                                        >  2:1572622/1‑100 (MQ=255)
              tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCCCGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtata                                                                                       >  2:4163442/1‑100 (MQ=255)
              tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgggacg                                                                                       <  1:3065858/100‑5 (MQ=18)
                gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtatacc                                                                                     <  2:2972646/100‑3 (MQ=255)
                gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                     <  1:19208/100‑1 (MQ=255)
                  cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACga                                                                                   >  2:3988400/1‑100 (MQ=255)
                  cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACg                                                                                    >  2:3636742/1‑99 (MQ=255)
                   ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCCTCCCGTCACGAT‑CTTCAGCAATA‑TGCTTAAGCATCGTGGCGGTATAACga                                                                                   <  1:337243/99‑1 (MQ=255)
                    tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatg                                                                                 <  2:328990/100‑1 (MQ=255)
                     acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatga                                                                                >  2:2302212/1‑100 (MQ=255)
                     acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatga                                                                                <  2:4272262/100‑1 (MQ=255)
                     acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatga                                                                                >  2:2730965/1‑100 (MQ=255)
                        ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                             >  2:3817547/1‑100 (MQ=255)
                        ccGTAATTTAATTGCAACCGGCTCGAGGCTGGGATACTGACTCTTCCAGTCACGAT‑CTTCCGAAATA‑TGCAGAAGCATCGTGGCGGGATAACCATGAtt                                                                              >  1:1083022/1‑99 (MQ=255)
                         cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                            <  2:490569/100‑1 (MQ=255)
                           taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                          <  2:2591694/100‑1 (MQ=255)
                           taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCAGTATAACGATGATTTCTg                                                                          <  1:4281909/100‑1 (MQ=255)
                           taatttaattGCAACCGGCTCGATGCTGGATTAATTACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACTATGATTTCTg                                                                          <  1:4093879/100‑1 (MQ=255)
                            aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                         <  2:2037611/100‑1 (MQ=255)
                              tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                       <  2:3586830/100‑1 (MQ=255)
                               ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                      >  2:2454157/1‑100 (MQ=255)
                                taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                      <  1:524716/99‑1 (MQ=255)
                                 aattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCGCGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATGACGATGATTTGTGATTGaa                                                                    >  1:1064302/1‑100 (MQ=255)
                                 aattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                    <  2:769865/100‑1 (MQ=255)
                                   ttGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATa                                                                  <  1:4058409/100‑1 (MQ=255)
                                     gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                >  2:3473605/1‑100 (MQ=255)
                                      cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                >  2:1215597/1‑99 (MQ=255)
                                      cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                <  2:2557048/99‑1 (MQ=255)
                                      cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                               >  2:1448451/1‑100 (MQ=255)
                                       aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                              <  1:392950/100‑1 (MQ=255)
                                       aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                              <  2:3383394/100‑1 (MQ=255)
                                        aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGat                                                             >  2:2343610/1‑100 (MQ=255)
                                        aCCGGCTCGATGCTGAGTTAGTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                 <  1:565380/96‑1 (MQ=255)
                                               cGATGCTGGGTTAATGACTCGTGCCGTCACGCT‑CTTCCGCAATA‑TGCCGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                      <  1:2975127/100‑1 (MQ=255)
                                               cGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                      >  1:1257663/1‑100 (MQ=255)
                                                 aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                    <  2:835516/100‑1 (MQ=255)
                                                  tGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCat                                                   <  2:2876793/100‑1 (MQ=255)
                                                    cTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCGGATTGAATAAAGATATTCCTGCatat                                                 >  2:3081739/1‑100 (MQ=255)
                                                      gggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGTGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                               <  1:2454146/100‑1 (MQ=255)
                                                          tAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                           >  2:1263499/1‑100 (MQ=255)
                                                          tAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                           >  1:1573102/1‑100 (MQ=255)
                                                           aaTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCgg                                          <  2:2936530/100‑1 (MQ=255)
                                                           aaTGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCgg                                          >  2:2823438/1‑100 (MQ=255)
                                                             tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGCTTTCTGATTGAATAAAGATATTCCTGCACATACATCGGcc                                        >  2:1901808/1‑100 (MQ=255)
                                                             tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                        <  1:1529604/100‑1 (MQ=255)
                                                                cTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgt                                     >  1:605193/1‑100 (MQ=255)
                                                                  cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                   <  2:1536270/100‑1 (MQ=255)
                                                                   gTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                  <  2:227886/100‑1 (MQ=255)
                                                                   gTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                  >  2:2790383/1‑100 (MQ=255)
                                                                      ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATACACATCGGCCTGTGCGGTg                               >  2:1118750/1‑100 (MQ=255)
                                                                      ccGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               >  1:4028990/1‑100 (MQ=25)
                                                                      ccGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               >  2:924739/1‑100 (MQ=21)
                                                                       cgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              >  1:1941481/4‑100 (MQ=25)
                                                                       cgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              >  2:3537737/4‑100 (MQ=25)
                                                                       cgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               >  2:3874330/4‑99 (MQ=21)
                                                                       cgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              >  1:2640033/4‑100 (MQ=21)
                                                                       cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              <  2:1821437/100‑1 (MQ=255)
                                                                        gccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                             >  2:829076/3‑100 (MQ=25)
                                                                        gTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              >  2:4231577/1‑99 (MQ=255)
                                                                        gTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                             <  2:411983/100‑1 (MQ=255)
                                                                          cACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                           >  1:3344003/1‑100 (MQ=37)
                                                                           acgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                          >  1:2531069/7‑100 (MQ=37)
                                                                           acgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                          >  2:3445743/7‑100 (MQ=25)
                                                                            cGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                         <  1:1391562/100‑1 (MQ=255)
                                                                             gatgcTTCAGC‑ATATTGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                         >  2:431574/5‑100 (MQ=255)
                                                                             gatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaa                        >  1:1803178/5‑100 (MQ=25)
                                                                             gAT‑CTTCAGCAATA‑TGCTGATGCATCGTGGCGGTATAACTATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTCTGCGGTGCATCaa                         <  1:3842337/99‑1 (MQ=255)
                                                                             gAT‑CTTCAGCAATA‑TGCTGAAGCATCTTGGCGGTATAACGATGATTTCTGTTTGAATAAAGATAATCCTGCATATACATCGGCCTGTGCGGTGCATCaaa                        >  1:511979/1‑100 (MQ=255)
                                                                              atgcTTCAGC‑ATATTGCTGAAG‑ATTGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       >  1:1185323/4‑100 (MQ=37)
                                                                              atgcTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       >  1:68605/4‑100 (MQ=38)
                                                                              aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       <  2:513666/100‑1 (MQ=255)
                                                                                      gCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                >  1:299690/1‑100 (MQ=255)
                                                                                      gCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                >  2:2759620/1‑100 (MQ=255)
                                                                                      gCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                <  1:995442/100‑1 (MQ=255)
                                                                                      gCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                >  1:1769907/1‑100 (MQ=255)
                                                                                       cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg               <  2:3212127/100‑1 (MQ=255)
                                                                                        aaTA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg               <  1:1345461/99‑1 (MQ=255)
                                                                                        aaTA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc              >  1:777971/1‑100 (MQ=255)
                                                                                          ta‑tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGt            >  1:4390631/1‑100 (MQ=255)
                                                                                               cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTccc        >  1:2191894/1‑100 (MQ=255)
                                                                                               cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc        >  2:1084173/1‑100 (MQ=255)
                                                                                                tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg       <  1:3253295/100‑1 (MQ=255)
                                                                                                  aaGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGa     <  1:27032/100‑1 (MQ=255)
                                                                                                   aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt    <  1:4408271/100‑1 (MQ=255)
                                                                                                    gCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAtt   >  2:2027484/1‑100 (MQ=255)
                                                                                                    gCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAtt   >  2:3262245/1‑100 (MQ=255)
                                                                                                     |                                                                                                   
TCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTC  >  NC_000913/3066220‑3066418

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: ‑