Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,319 Δ1 bp 19.3% 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.19.3% 544.3 / 137.2 192coding (855/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (83/71);  new base . (18/19);  total (102/90)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

TCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTC  >  NC_000913/3066220‑3066418
                                                                                                     |                                                                                                   
 cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCa                                                                                                    >  1:1027039/1‑100 (MQ=255)
    aGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCg                                                                                                 <  2:760567/100‑1 (MQ=255)
    aGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGC‑ATATTGCTGAAG‑ATCGt                                                                                                <  1:3764666/100‑1 (MQ=255)
     gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGt                                                                                                <  2:668686/100‑1 (MQ=25)
      ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTg                                                                                               <  1:854942/100‑1 (MQ=255)
       cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTgg                                                                                              >  1:4008215/1‑100 (MQ=255)
        aCGTCTTCGACTTACTCCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGc                                                                                             <  2:3683951/100‑1 (MQ=255)
        aCGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGc                                                                                             <  1:1727084/100‑1 (MQ=255)
        aCGTCTGCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTgg                                                                                              <  2:3271963/99‑1 (MQ=255)
         cGTCTTCGACTTACATCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                            <  1:4078922/100‑1 (MQ=21)
         cGTCTTCGACTTACACCGTACTTTAATTTCAACCGGCTTGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACg                                                                                            <  2:2730161/100‑1 (MQ=18)
         cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                            <  1:1884120/100‑1 (MQ=21)
         cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                            <  1:2097601/100‑1 (MQ=21)
         cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                            <  1:2152544/100‑1 (MQ=21)
          ttataccaCTTACACCATAATTTACTGGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCAACGATGCTACTGCATAT‑TGCTGAAG‑ATCGTGACgg                                                                                           <  2:4063576/93‑1 (MQ=12)
          gTCTTCGACTTACACCGTAATTTAATTGCAGCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACg                                                                                            <  2:543589/99‑1 (MQ=21)
          gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                           <  2:2937083/100‑1 (MQ=255)
          gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                           >  1:294614/1‑100 (MQ=255)
          gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                           <  1:131721/100‑1 (MQ=255)
          gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACGCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACgg                                                                                           <  2:1264173/100‑1 (MQ=21)
           tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGt                                                                                          >  1:3283264/1‑100 (MQ=255)
            cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         <  1:2025460/100‑1 (MQ=255)
            cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         <  2:1898727/100‑1 (MQ=255)
            cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         <  1:3080500/100‑1 (MQ=255)
            cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                         >  2:2614456/1‑100 (MQ=255)
             ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGtat                                                                                        <  1:2344680/100‑1 (MQ=255)
              tCGACTTACACCGTAATTTAATTGCAACCGGCTCGGTGTTGGGTTAATGACTCGTTACGTCAAGAT‑CTTCAGCAATA‑TTGTGAAGACTCGTGTGGGtat                                                                                        >  2:1920790/1‑99 (MQ=255)
                gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                     >  2:554838/1‑100 (MQ=255)
                gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                     >  2:3185112/1‑100 (MQ=255)
                gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                     >  2:1297846/1‑100 (MQ=255)
                gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                     <  2:2079165/100‑1 (MQ=255)
                  cTTACACCGTGATTTCATTGCACCCCGCTCGATGCAGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACggg                                                                                          <  2:223612/93‑2 (MQ=9)
                  cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                     >  2:858951/1‑98 (MQ=255)
                  cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACga                                                                                   >  1:606856/1‑100 (MQ=255)
                    tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGATGAAGCATCGTGGCGGTATAACgatg                                                                                 >  1:1914640/1‑100 (MQ=255)
                        ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                             >  1:2866943/1‑100 (MQ=255)
                         cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                             >  2:3269239/1‑99 (MQ=255)
                          gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                           <  2:1034246/100‑1 (MQ=255)
                          gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                           <  2:3286348/100‑1 (MQ=255)
                              tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                       >  2:1341103/1‑100 (MQ=255)
                               ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGTTGATTTCTGATTg                                                                      <  1:2808052/100‑1 (MQ=255)
                                  attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTGTAACGATGATTTCTGATTGAAt                                                                   >  2:3071404/1‑100 (MQ=255)
                                  attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                   >  1:2258076/1‑100 (MQ=255)
                                     gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                >  1:3140840/1‑100 (MQ=255)
                                      cAACCGGCTCGATGCTGGGTTAATGACTCGTCGCGTCACGAT‑CTTCAGCAATA‑TGGTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                               >  1:3927616/1‑100 (MQ=255)
                                       aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                              >  1:3576182/1‑100 (MQ=255)
                                            gCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                         >  1:1716766/1‑100 (MQ=255)
                                                 aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc                                                     <  2:902965/99‑1 (MQ=255)
                                                 aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                    <  1:2597810/100‑1 (MQ=255)
                                                 aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                    >  1:3527278/1‑100 (MQ=255)
                                                    cTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATGTCTGATTGAATAAAGATATTCCTGCatat                                                 >  2:1739484/1‑100 (MQ=255)
                                                          tAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                            >  1:478580/1‑99 (MQ=255)
                                                             tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                        <  1:1579925/100‑1 (MQ=255)
                                                               aCTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtg                                      <  2:3708649/100‑1 (MQ=255)
                                                                  cGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                   >  1:1061413/1‑100 (MQ=255)
                                                                    tCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                 <  1:333174/100‑1 (MQ=255)
                                                                     cccGTCACGAT‑CTTCAGCAATA‑TGATGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAGATCGGCCtgtg                                    >  2:19607/1‑96 (MQ=255)
                                                                      ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                                >  1:77580/1‑99 (MQ=255)
                                                                      ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               <  1:1021072/100‑1 (MQ=255)
                                                                      ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               >  1:2417157/1‑100 (MQ=255)
                                                                      ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               >  1:2208635/1‑100 (MQ=255)
                                                                      ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                               >  2:1024689/1‑100 (MQ=255)
                                                                      ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTCCGGTg                               >  2:3499516/1‑100 (MQ=255)
                                                                       cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                              <  1:3462728/100‑1 (MQ=255)
                                                                        gccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                             >  1:3811469/3‑100 (MQ=25)
                                                                         tCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCAt                            >  1:438568/1‑100 (MQ=255)
                                                                          cACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                           >  2:4041543/1‑100 (MQ=255)
                                                                            cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                         >  1:2636645/6‑100 (MQ=25)
                                                                            cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                         >  1:1173122/6‑100 (MQ=25)
                                                                            cGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                         >  2:2115342/1‑100 (MQ=255)
                                                                            cGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                          >  2:43143/1‑99 (MQ=255)
                                                                              aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       >  1:1136233/1‑100 (MQ=255)
                                                                              aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       <  1:3370521/100‑1 (MQ=255)
                                                                              aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       <  1:2749913/100‑1 (MQ=255)
                                                                              aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       <  2:2835550/100‑1 (MQ=255)
                                                                              aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                       <  1:3681138/100‑1 (MQ=255)
                                                                                  ttCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                    >  2:3952308/1‑100 (MQ=255)
                                                                                  ttCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                    >  2:1988111/1‑100 (MQ=255)
                                                                                   tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                   <  1:3175544/100‑1 (MQ=255)
                                                                                   tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                   >  2:1121422/1‑100 (MQ=255)
                                                                                   tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                   >  1:944020/1‑100 (MQ=255)
                                                                                     cagcatattGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttag                 >  2:3767152/9‑98 (MQ=38)
                                                                                      gCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                >  1:3654058/1‑100 (MQ=255)
                                                                                       cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg               <  1:677681/100‑1 (MQ=255)
                                                                                       cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg               <  1:372427/100‑1 (MQ=255)
                                                                                         atattGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATCCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg             >  1:2698281/5‑100 (MQ=255)
                                                                                           a‑tGCTGCAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCCGCGt            <  2:138944/99‑1 (MQ=255)
                                                                                             tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg          <  2:3561580/100‑1 (MQ=255)
                                                                                             tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg          >  1:1448616/1‑100 (MQ=255)
                                                                                             tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg          >  1:1472557/1‑100 (MQ=255)
                                                                                              gCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGc         >  1:3349511/1‑100 (MQ=255)
                                                                                               cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc        >  2:2186158/1‑100 (MQ=255)
                                                                                                tGAAGCATCGTGGCGGTGTAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg       <  2:2529457/100‑1 (MQ=255)
                                                                                                 gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg      >  1:2477566/1‑100 (MQ=255)
                                                                                                 gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg      <  2:3842633/100‑1 (MQ=255)
                                                                                                   aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGa     >  1:903651/1‑99 (MQ=255)
                                                                                                   aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt    <  2:474561/100‑1 (MQ=255)
                                                                                                   aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt    <  1:543972/100‑1 (MQ=255)
                                                                                                   aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt    <  1:514731/100‑1 (MQ=255)
                                                                                                     cATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTc  <  1:2045378/100‑1 (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: ‑