Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,294 G→C 19.8% D292H (GAT→CAT)  ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,2940GC19.8% 449.7 / 101.9 187D292H (GAT→CAT) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (90/60);  new base C (19/18);  total (109/78)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 196 total aligned reads displayed.

GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCG  >  NC_000913/3181195‑3181392
                                                                                                   |                                                                                                  
  gTGGAATAGCATTATCTACAACAAGGCACTTCTACCGCTACTGTGGGATATGTTCCTGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                  <  2:1823475/100‑1 (MQ=255)
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                  <  1:1023703/100‑1 (MQ=255)
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                  >  1:4049415/1‑100 (MQ=255)
     gAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcat                                                                                               >  1:669347/1‑100 (MQ=255)
       agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                             >  2:1662739/1‑100 (MQ=255)
        gagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcatc                                                                                              <  2:4196423/98‑1 (MQ=255)
        gagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCg                                                                                            >  2:3079627/1‑100 (MQ=255)
            aTTATCTCCAACAAGTCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaa                                                                                        <  1:4133876/100‑1 (MQ=255)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                      <  1:1737571/100‑3 (MQ=38)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                      <  1:722605/100‑1 (MQ=255)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                      <  1:2113020/100‑1 (MQ=255)
               aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                     <  1:4001941/100‑4 (MQ=39)
                tctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGa                                                                                    >  2:2220582/1‑100 (MQ=255)
                tctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGa                                                                                    >  2:1716799/1‑100 (MQ=255)
                 ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataag                                                                                   <  1:1141828/100‑6 (MQ=38)
                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                                  >  2:1567907/1‑100 (MQ=255)
                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGAAAATGGaaa                                                                                  >  1:1302063/1‑100 (MQ=255)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                                 <  1:1904127/100‑8 (MQ=37)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaa                                                                                 >  2:1207596/1‑100 (MQ=255)
                    caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaaa                                                                                >  1:2508904/1‑100 (MQ=255)
                    caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaaa                                                                                >  1:2613309/1‑100 (MQ=255)
                     aacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcggg                                                                               <  2:3589354/100‑10 (MQ=37)
                     aacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAACTGGaaaaa                                                                                >  1:1794299/1‑99 (MQ=255)
                      acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                              <  2:4032433/100‑1 (MQ=255)
                       caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                             >  2:682547/1‑100 (MQ=255)
                       caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                             >  1:424866/1‑100 (MQ=255)
                       caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                             >  1:1929784/1‑100 (MQ=255)
                        aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                            <  2:1951185/100‑1 (MQ=255)
                         aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                           >  1:3258075/1‑100 (MQ=255)
                         aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                           >  2:2161894/1‑100 (MQ=255)
                           gCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTGTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgtg                                                                          <  2:1021240/99‑1 (MQ=255)
                             aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                       <  2:1729790/100‑1 (MQ=255)
                              cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTa                                                                      <  2:101311/100‑1 (MQ=255)
                              cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCACCCGCAAATGGAAAAATATGTGGTTa                                                                      >  2:1310602/1‑100 (MQ=255)
                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                  >  2:1145701/1‑100 (MQ=255)
                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                  >  1:497924/1‑100 (MQ=255)
                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                  >  2:3931102/1‑100 (MQ=255)
                                  tACCGCCACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                  <  2:249319/100‑1 (MQ=255)
                                    ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                >  2:3168741/1‑100 (MQ=255)
                                    ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTACCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                >  2:4166284/1‑100 (MQ=255)
                                      gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                              <  2:846799/100‑1 (MQ=255)
                                      gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                              <  1:1511846/100‑1 (MQ=255)
                                      gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                              <  1:2974350/100‑1 (MQ=255)
                                          cTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACGGAtcttc                                                           >  2:2909849/1‑99 (MQ=255)
                                           tgtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCACATGGAAAAATATGTGGTTAAACCGATCTtctc                                                         >  2:1842711/1‑100 (MQ=255)
                                              gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                      >  1:2220756/1‑100 (MQ=255)
                                              gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                      >  1:3949502/1‑100 (MQ=255)
                                              gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                      >  1:538793/1‑100 (MQ=255)
                                                    tgtcgccGAAGCACGCGAACCTGCTGCCTGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAAGATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                <  2:3635023/95‑1 (MQ=255)
                                                    tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                >  2:3117692/1‑100 (MQ=255)
                                                    tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                <  2:4125934/100‑1 (MQ=255)
                                                    tGTTCCCGAATCACCCGAACCTGCAGCCCGCTTATTTTGCGGAAGATGATCACCCGCAAATGGGAAAATACGTGGATAAACCGATCTTCTCCCGTGAAgg                                                >  2:2787335/1‑100 (MQ=255)
                                                     gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                               >  1:2157820/1‑100 (MQ=255)
                                                     gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                               >  1:2560452/1‑100 (MQ=255)
                                                       tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                             >  1:63244/1‑100 (MQ=255)
                                                       tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                             <  1:2644817/100‑1 (MQ=255)
                                                       tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                             >  1:2933821/1‑100 (MQ=255)
                                                        cccGAATCACCCGAACCTTCGGCCCGCTTATTTTGCGGAATATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                            >  2:2443888/1‑100 (MQ=255)
                                                         ccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaa                                           >  2:3816130/1‑100 (MQ=255)
                                                            aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                        <  1:3665504/100‑1 (MQ=255)
                                                              tCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtg                                      <  1:2487136/100‑1 (MQ=255)
                                                               cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                     >  2:1432450/1‑100 (MQ=255)
                                                               cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                     >  2:2977646/1‑100 (MQ=255)
                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                   >  1:3392314/1‑100 (MQ=255)
                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCAACCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                   <  2:1613377/100‑1 (MQ=255)
                                                                   cGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGAt                                 >  2:3266886/1‑100 (MQ=255)
                                                                   cGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGAt                                 >  2:2582437/1‑100 (MQ=255)
                                                                        ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                            >  2:620356/1‑100 (MQ=255)
                                                                         tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                           >  1:1218215/1‑100 (MQ=255)
                                                                            tgaccGCTTATTTTGCGGATGATGTTCATCCACAAATGGAAAAATATGAGGTTAAACCGATCTTCTCCCGTGATGGCGGAAACGTGTCGATCATTg                            <  2:3602711/93‑1 (MQ=255)
                                                                             gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                       >  1:1288356/1‑100 (MQ=255)
                                                                              cccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                      >  1:3733937/1‑100 (MQ=255)
                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                    >  1:3467010/1‑100 (MQ=255)
                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                    <  1:2729958/100‑1 (MQ=255)
                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                    >  1:3509414/1‑100 (MQ=255)
                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                    >  1:3141350/1‑100 (MQ=255)
                                                                                   ttATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa                 >  1:1472772/1‑100 (MQ=255)
                                                                                    tATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                <  1:2801996/100‑1 (MQ=255)
                                                                                    ccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                >  1:354159/4‑100 (MQ=38)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCCAAGGTGTCGATCATTAAGAACGGCtaa                 >  2:1107277/3‑98 (MQ=37)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc               >  1:3844561/3‑100 (MQ=38)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc               >  1:3329323/3‑100 (MQ=38)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAAACGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc               >  2:3796302/3‑100 (MQ=38)
                                                                                     catttGCGGATGATCATCTTCCGCAAAGGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAGAc               >  1:3608752/3‑100 (MQ=38)
                                                                                     catttGCGGATGAACATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc               >  1:2846758/3‑100 (MQ=255)
                                                                                     aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc               >  1:1096019/1‑100 (MQ=255)
                                                                                     aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCGCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc               >  2:2682720/1‑100 (MQ=255)
                                                                                      atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc              >  2:91612/2‑100 (MQ=39)
                                                                                         tGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt           >  2:3449743/1‑100 (MQ=255)
                                                                                          gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg          <  1:841851/100‑1 (MQ=255)
                                                                                          gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg          <  2:669346/100‑1 (MQ=255)
                                                                                           cGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTAAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGa         >  2:2245928/1‑100 (MQ=255)
                                                                                            gggaGATGATCATCCGCAAATTGAAAAATATGGGGTTAATCCGCTCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGaa        <  2:2380980/97‑1 (MQ=255)
                                                                                            ggAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGaa        >  2:3395234/1‑100 (MQ=255)
                                                                                            ggAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGaa        >  2:2742946/1‑100 (MQ=255)
                                                                                            ggAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGaa        >  1:2528166/1‑100 (MQ=255)
                                                                                             gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag       >  1:88285/1‑100 (MQ=255)
                                                                                             gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag       <  1:3074243/100‑1 (MQ=255)
                                                                                             gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag       >  1:2933947/1‑100 (MQ=255)
                                                                                              aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagc      <  1:113195/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCG  >  NC_000913/3181195‑3181392

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