Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,534 G→C 8.1% A217G (GCC→GGC)  yjhP ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,5340GC 8.1% 431.2 / 23.4 124A217G (GCC→GGC) yjhPputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (54/60);  new base C (5/5);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 136 total aligned reads displayed.

GCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGG  >  NC_000913/4532387‑4532682
                                                                                                                                                   |                                                                                                                                                    
gcaaagGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAAGCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggc                                                                                                                                                     >  4:1019548/1‑149 (MQ=255)
   aagGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCggcggc                                                                                                                                                  >  3:634319/1‑149 (MQ=255)
     ggTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                                <  1:1034770/149‑1 (MQ=255)
      gTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGaa                                                                                                                                               <  4:319778/149‑1 (MQ=255)
      gTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGaa                                                                                                                                               <  4:211415/149‑1 (MQ=255)
        aaaGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGt                                                                                                                                             >  4:1262026/1‑149 (MQ=255)
            ccTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                         <  2:202954/149‑1 (MQ=255)
             cTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                         >  2:669626/1‑148 (MQ=255)
             cTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                         <  1:669626/148‑1 (MQ=255)
             cTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCg                                                                                                                                        >  3:695479/1‑149 (MQ=255)
              tCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                       >  3:1079527/1‑149 (MQ=255)
              tCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                       >  4:880337/1‑149 (MQ=255)
                gATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                                                                     <  1:1308240/149‑1 (MQ=255)
                gATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                                                                     <  1:300351/149‑1 (MQ=255)
                      tCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTc                                                                                                                               <  4:1109624/149‑1 (MQ=255)
                          gccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                                                                           >  4:1103936/1‑149 (MQ=255)
                          gccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCa                                                                                                                                           >  2:107061/1‑133 (MQ=255)
                           ccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                                                                          >  2:382413/1‑149 (MQ=255)
                           ccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                                                                          >  1:1389202/1‑149 (MQ=255)
                           ccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCAATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                                                                          <  1:33617/149‑1 (MQ=255)
                                  tACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGc                                                                                                                             <  2:626096/139‑1 (MQ=255)
                                  tACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGc                                                                                                                             >  1:626096/1‑139 (MQ=255)
                                  tACCGGGAACCTGTTGTTACCGGACGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagccag                                                                                                                   <  3:102036/149‑1 (MQ=255)
                                      gggAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcggc                                                                                                               >  4:61629/1‑149 (MQ=255)
                                      gggAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcggc                                                                                                               <  3:1041351/149‑1 (MQ=255)
                                       ggAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                                                                      <  3:566341/125‑1 (MQ=255)
                                       ggAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                                                                      >  4:566341/1‑125 (MQ=255)
                                            cTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGt                                                                                                                                             <  2:1153256/113‑1 (MQ=255)
                                            cTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGt                                                                                                                                             >  1:1153256/1‑113 (MQ=255)
                                            cTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                                                         <  4:191481/133‑1 (MQ=255)
                                            cTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                                                         >  3:191481/1‑133 (MQ=255)
                                            cTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATgg                                                                                                         >  1:921876/1‑149 (MQ=255)
                                                 gttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                                                                    >  3:499274/1‑149 (MQ=255)
                                                     ccGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtt                                                                                                <  1:1088037/149‑1 (MQ=255)
                                                     ccGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtt                                                                                                >  3:42077/1‑149 (MQ=255)
                                                     ccGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtt                                                                                                >  1:464091/1‑149 (MQ=255)
                                                         gCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                                                            <  3:800839/149‑1 (MQ=255)
                                                             tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCt                                                                                        <  4:5263/149‑1 (MQ=255)
                                                              gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                                       <  1:382413/149‑1 (MQ=255)
                                                                gTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                                       >  1:575008/1‑147 (MQ=255)
                                                                gTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                                       <  2:575008/147‑1 (MQ=255)
                                                                 tGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCAt                                                                                    >  3:927201/1‑149 (MQ=255)
                                                                    gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                                                                          <  3:1378573/108‑1 (MQ=255)
                                                                    gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGgcgc                                                                                                             <  4:689272/121‑1 (MQ=255)
                                                                    gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAt                                                                                                 >  2:248446/1‑133 (MQ=255)
                                                                    gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAt                                                                                                 <  1:248446/133‑1 (MQ=255)
                                                                    gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAcc                                                                                 <  4:499274/149‑1 (MQ=255)
                                                                    gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGATTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                                                                          >  4:1378573/1‑108 (MQ=255)
                                                                            cccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggc                                                                                                                                                     <  2:32693/73‑1 (MQ=255)
                                                                                gCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                         >  4:1291781/1‑81 (MQ=255)
                                                                                gCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                         <  3:1291781/81‑1 (MQ=255)
                                                                                  cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTcc                                                                   >  3:1104848/1‑149 (MQ=255)
                                                                                   aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                                  <  2:526384/149‑1 (MQ=255)
                                                                                   aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                                  >  1:709624/1‑149 (MQ=255)
                                                                                   aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                                  >  3:1189883/1‑149 (MQ=255)
                                                                                    aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                                                                 <  4:1063117/149‑1 (MQ=255)
                                                                                     aGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtttt                                                                                              <  3:1152505/119‑1 (MQ=255)
                                                                                      gCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt                                                               <  3:880337/149‑1 (MQ=255)
                                                                                       cATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                                                              <  4:927201/149‑1 (MQ=255)
                                                                                       cATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                                                              >  1:550271/1‑149 (MQ=255)
                                                                                        aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                                                                     >  1:1060753/1‑77 (MQ=255)
                                                                                        aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                                                                     <  2:1060753/77‑1 (MQ=255)
                                                                                        aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                                                             <  2:1389202/149‑1 (MQ=255)
                                                                                            ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCa                                                         <  4:695479/149‑1 (MQ=255)
                                                                                             cgcgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                                                         >  3:61087/1‑84 (MQ=255)
                                                                                             cgcgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                                                         <  4:61087/84‑1 (MQ=255)
                                                                                              gcgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGt                                                       <  3:1262026/149‑1 (MQ=255)
                                                                                               cgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTa                                                      <  4:1189883/149‑1 (MQ=255)
                                                                                                   tAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAt                                                  <  2:550271/149‑1 (MQ=255)
                                                                                                       tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTcc                                              <  1:663359/149‑1 (MQ=255)
                                                                                                       tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTcc                                              >  1:568957/1‑149 (MQ=255)
                                                                                                           gTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCg                                                                                <  4:845326/111‑1 (MQ=255)
                                                                                                            tAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCcaccac                                         >  4:1153137/1‑149 (MQ=255)
                                                                                                                       cgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCga                              <  2:921876/149‑1 (MQ=255)
                                                                                                                          gATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTc                                                                                                                               >  1:126541/1‑49 (MQ=255)
                                                                                                                           aTGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtc                          <  3:1103936/149‑1 (MQ=255)
                                                                                                                              tttAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCa                                                                                                      <  1:235609/70‑1 (MQ=255)
                                                                                                                                tAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                                                            <  2:996804/78‑1 (MQ=255)
                                                                                                                                       gCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                                                         >  2:1240254/1‑42 (MQ=255)
                                                                                                                                       gCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTcc              <  4:1104848/149‑1 (MQ=255)
                                                                                                                                        cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCg             >  4:275586/1‑149 (MQ=255)
                                                                                                                                         gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGaa                    >  2:166970/3‑141 (MQ=255)
                                                                                                                                         gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGa            >  3:82195/3‑149 (MQ=255)
                                                                                                                                         gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTACCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCga                              >  1:103062/3‑131 (MQ=255)
                                                                                                                                         gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTACCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCga                              <  2:103062/129‑1 (MQ=255)
                                                                                                                                         ccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTcc                                                                   >  2:578790/1‑94 (MQ=255)
                                                                                                                                          cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCTATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCAt                                    >  3:850469/1‑124 (MQ=255)
                                                                                                                                           tGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAttt                                                                                               <  4:1172568/64‑1 (MQ=255)
                                                                                                                                           tGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAttt                                                                                               >  3:1172568/1‑64 (MQ=255)
                                                                                                                                           tGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                                                  <  3:195495/77‑1 (MQ=255)
                                                                                                                                           tGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                                                  >  4:195495/1‑77 (MQ=255)
                                                                                                                                           tGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACa          >  3:239100/1‑149 (MQ=255)
                                                                                                                                             aCTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                                                                 <  1:519056/92‑1 (MQ=255)
                                                                                                                                             aCTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                                                                 >  2:519056/1‑92 (MQ=255)
                                                                                                                                                  ggcggcGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                                                              >  2:422733/1‑90 (MQ=255)
                                                                                                                                                  ggcggcGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                                                              <  1:422733/90‑1 (MQ=255)
                                                                                                                                                   gcggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                                                             >  1:110305/1‑90 (MQ=255)
                                                                                                                                                   gcggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTgg  <  4:42077/149‑1 (MQ=255)
                                                                                                                                                   gcggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTgg  >  2:406789/1‑149 (MQ=255)
                                                                                                                                                   gcggcgAAGTCATCGTCAGGATTCGCTACCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                                                             <  2:110305/90‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGG  >  NC_000913/4532387‑4532682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑