Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,298 A→T 19.4% H293L (CAT→CTT)  ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,2980AT19.4% 434.6 / 85.9 188H293L (CAT→CTT) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base A (90/61);  new base T (19/18);  total (109/79)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

GAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGG  >  NC_000913/3181200‑3181397
                                                                                                  |                                                                                                   
  agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                                  <  1:643958/100‑1 (MQ=255)
   gagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcatc                                                                                                   <  2:4196423/98‑1 (MQ=255)
   gagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCg                                                                                                 >  2:3079627/1‑100 (MQ=255)
       aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                             <  2:2654344/100‑1 (MQ=255)
        ttATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaa                                                                                            <  1:2677189/100‑2 (MQ=39)
        ttATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAAt                                                                                            <  1:3495374/100‑1 (MQ=255)
         tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                           <  2:244288/100‑3 (MQ=39)
         tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                           <  1:2113020/100‑1 (MQ=255)
          aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                          <  1:4164305/100‑4 (MQ=38)
            ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataag                                                                                        <  1:1141828/100‑6 (MQ=38)
             tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                                       <  1:2595743/100‑7 (MQ=38)
             tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                                       <  1:494879/100‑1 (MQ=255)
             tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGAAAATGGaaa                                                                                       >  1:1302063/1‑100 (MQ=255)
              ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                                      <  1:514245/100‑8 (MQ=37)
               caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaaa                                                                                     >  1:4021266/1‑100 (MQ=255)
               caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaaa                                                                                     >  1:2508904/1‑100 (MQ=255)
                aacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcggg                                                                                    <  2:3589354/100‑10 (MQ=37)
                 acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                                   <  2:4032433/100‑1 (MQ=255)
                 acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTCATTTTGCGGATGATCATCTTCCGCaaaataagcgggc                                                                                   <  2:2241059/100‑11 (MQ=25)
                  caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                                  >  2:682547/1‑100 (MQ=255)
                  caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                                  >  2:2325466/1‑100 (MQ=255)
                   aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                                 <  2:1951185/100‑1 (MQ=255)
                    aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                                >  2:2161894/1‑100 (MQ=255)
                    aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                                >  1:3258075/1‑100 (MQ=255)
                     ggCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgtg                                                                               <  1:2873873/100‑1 (MQ=255)
                     ggCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgtg                                                                               >  1:183436/1‑100 (MQ=255)
                         cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtca                                                                           <  1:2846283/100‑3 (MQ=255)
                            cTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                                        <  1:2633900/100‑1 (MQ=255)
                             tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                       >  2:3931102/1‑100 (MQ=255)
                              aCCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCg                                                                      <  2:2688152/100‑1 (MQ=255)
                               ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                     >  2:3168741/1‑100 (MQ=255)
                               ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTACCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                     >  2:4166284/1‑100 (MQ=255)
                                cGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAt                                                                    >  1:395695/1‑100 (MQ=255)
                                 gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                                   <  1:2974350/100‑1 (MQ=255)
                                 gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                                   <  2:846799/100‑1 (MQ=255)
                                 gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCGGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                                   <  2:2489563/100‑1 (MQ=255)
                                     cTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACGGAtcttc                                                                >  2:2909849/1‑99 (MQ=255)
                                      tgtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCACATGGAAAAATATGTGGTTAAACCGATCTtctc                                                              >  2:1842711/1‑100 (MQ=255)
                                       gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                             >  1:3594296/1‑100 (MQ=255)
                                       gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                             <  2:1207723/100‑1 (MQ=255)
                                         gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                           >  1:2220756/1‑100 (MQ=255)
                                           gagaTGTTCCCGAATCACCCGAACCTGCTCCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGTAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                         <  2:2523523/100‑1 (MQ=255)
                                               tgtcgccGAAGCACGCGAACCTGCTGCCTGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAAGATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                     <  2:3635023/95‑1 (MQ=255)
                                               tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                     >  2:3117692/1‑100 (MQ=255)
                                               tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                     <  1:3617558/100‑1 (MQ=255)
                                               tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                     >  1:542183/1‑100 (MQ=255)
                                               tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                     <  1:3523652/100‑1 (MQ=255)
                                               tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                     <  2:4125934/100‑1 (MQ=255)
                                               tGTTCCCGAATCACCCGAACCTGCAGCCCGCTTATTTTGCGGAAGATGATCACCCGCAAATGGGAAAATACGTGGATAAACCGATCTTCTCCCGTGAAgg                                                     >  2:2787335/1‑100 (MQ=255)
                                                gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                                    >  1:2157820/1‑100 (MQ=255)
                                                   cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                                 >  1:2522162/1‑100 (MQ=255)
                                                    ccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaa                                                >  2:3816130/1‑100 (MQ=255)
                                                         tCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtg                                           <  1:2487136/100‑1 (MQ=255)
                                                          cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                          >  1:1416867/1‑100 (MQ=255)
                                                          cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                          >  2:2977646/1‑100 (MQ=255)
                                                           aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                         <  1:2744666/100‑1 (MQ=255)
                                                            cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTCAACCGATCTTCTCCCGTGAACGCGCAAACGTGTCg                                        >  1:1572074/1‑100 (MQ=255)
                                                            cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                        >  1:3392314/1‑100 (MQ=255)
                                                             ccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGa                                       >  1:827532/1‑100 (MQ=255)
                                                              cGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGAt                                      >  2:3266886/1‑100 (MQ=255)
                                                                 aCCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAt                                   <  2:2173295/100‑1 (MQ=255)
                                                                   ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                                 >  2:620356/1‑100 (MQ=255)
                                                                   ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                                 <  2:1698777/100‑1 (MQ=255)
                                                                    tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                                >  2:1649250/1‑100 (MQ=255)
                                                                    tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                                <  1:917212/100‑1 (MQ=255)
                                                                    tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                                >  1:1218215/1‑100 (MQ=255)
                                                                       tgaccGCTTATTTTGCGGATGATGTTCATCCACAAATGGAAAAATATGAGGTTAAACCGATCTTCTCCCGTGATGGCGGAAACGTGTCGATCATTg                                 <  2:3602711/93‑1 (MQ=255)
                                                                        gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                            >  1:1288356/1‑100 (MQ=255)
                                                                           cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGACAAATATGTGGCTAAACCGCTCTTCTCCCGTGAAGGCGCAACCCTGTCGATCATTGAGAACg                           >  2:1099159/1‑98 (MQ=255)
                                                                           cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                         >  1:3141350/1‑100 (MQ=255)
                                                                           cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                         >  1:3905368/1‑100 (MQ=255)
                                                                           cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                         >  1:3467010/1‑100 (MQ=255)
                                                                            tttccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCa                        >  2:159339/7‑100 (MQ=37)
                                                                            gCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCa                        >  2:113195/1‑100 (MQ=255)
                                                                              tccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa                      >  1:7144/5‑100 (MQ=38)
                                                                               ccatttGCGGATGATCATCTTCCGCGAATGGAAAAATATGTGGTTAAACCGATCGTCTCCCGTGAAGGCGCAAAGGTGACAATCATTGAGAACGGCaga                      >  2:1017302/4‑97 (MQ=37)
                                                                                catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                    >  1:3329323/3‑100 (MQ=38)
                                                                                catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                    >  1:1506575/3‑100 (MQ=38)
                                                                                catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAAACGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                    >  2:3796302/3‑100 (MQ=38)
                                                                                catttGCGGATGATCATCTTCCGCAAATGGAAAAAGATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTAGATCATTGAGAACGGCAAAAc                    >  2:1771656/3‑100 (MQ=38)
                                                                                aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                    >  1:1096019/1‑100 (MQ=255)
                                                                                 atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc                   >  2:91612/2‑100 (MQ=39)
                                                                                   ttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAt                 <  2:1506697/100‑1 (MQ=255)
                                                                                    tGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt                >  2:3449743/1‑100 (MQ=255)
                                                                                    tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt                <  1:1909930/100‑1 (MQ=255)
                                                                                     gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg               <  1:37775/100‑1 (MQ=255)
                                                                                     gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg               <  2:669346/100‑1 (MQ=255)
                                                                                       ggAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGaa             >  2:3395234/1‑100 (MQ=255)
                                                                                        gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag            >  1:88285/1‑100 (MQ=255)
                                                                                        gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag            >  1:372303/1‑100 (MQ=255)
                                                                                        gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag            <  1:3074243/100‑1 (MQ=255)
                                                                                        gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag            >  1:2933947/1‑100 (MQ=255)
                                                                                         aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagc           <  2:1190429/100‑1 (MQ=255)
                                                                                         aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagc           <  1:113195/100‑1 (MQ=255)
                                                                                             tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg       <  2:1722949/100‑1 (MQ=255)
                                                                                             tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg       >  2:1398050/1‑100 (MQ=255)
                                                                                             tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg       <  1:2544465/100‑1 (MQ=255)
                                                                                             tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg       >  1:1791148/1‑100 (MQ=255)
                                                                                                 catcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAg   >  1:2926565/1‑100 (MQ=255)
                                                                                                  atcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAgg  >  2:3456444/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGG  >  NC_000913/3181200‑3181397

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