Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 606,958 Δ1 bp 5.1% coding (426/1119 nt) ybdK ← weak gamma‑glutamyl:cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913606,9580G. 5.1% 269.8 / 7.8 79coding (426/1119 nt)ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (47/28);  new base . (2/2);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATT  >  NC_000913/606810‑607105
                                                                                                                                                    |                                                                                                                                                   
tATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTg                                                                                                                                                     >  8:215135/1‑149 (MQ=255)
    aGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGgcg                                                                                                                                                  >  6:316297/1‑148 (MQ=255)
    aGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGgcg                                                                                                                                                  >  2:110118/1‑148 (MQ=255)
             ggAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACACAGCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCCCTGGCGCAGCCAACAt                                                                                                                                        >  3:223786/1‑149 (MQ=255)
                 aaaaTATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGAc                                                                                                                                    <  6:394332/149‑1 (MQ=255)
                 aaaaTATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGAc                                                                                                                                    <  7:309185/149‑1 (MQ=255)
                                gaggagGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAAc                                                                                                                     >  3:221173/1‑149 (MQ=255)
                                gaggagGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAAc                                                                                                                     >  5:309519/1‑149 (MQ=255)
                                      gCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGc                                                                                                               <  1:62048/149‑1 (MQ=255)
                                        aaaaCGCGTATCCGTTCCCTGCATATATGGCGACGCGTCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGCCCAAAAACGGTCTCct                                                                                                             >  7:375013/1‑149 (MQ=255)
                                        aaaaCGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGCCATGCTGACCAAAAACGGTCTCct                                                                                                             >  3:26442/1‑149 (MQ=255)
                                        aaaaCGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCct                                                                                                             >  5:99372/1‑149 (MQ=255)
                                        aaaaCGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCct                                                                                                             >  2:5602/1‑149 (MQ=255)
                                        aaaaCGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCct                                                                                                             >  2:323651/1‑149 (MQ=255)
                                        aaaaCGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCct                                                                                                             >  2:103551/1‑149 (MQ=255)
                                               gTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAAt                                                                                                      >  6:229958/1‑149 (MQ=255)
                                                 aTCCGTTCCCTGCATATATGGCGCCGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAAt                                                                                                      >  2:116103/1‑147 (MQ=255)
                                                  tCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAAt                                                                                                      >  3:256487/1‑146 (MQ=255)
                                                    cGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATgaga                                                                                                  <  3:130697/148‑1 (MQ=255)
                                                    cGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGAt                                                                                                 <  2:362215/149‑1 (MQ=255)
                                                     gTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATa                                                                                                <  7:215135/149‑1 (MQ=255)
                                                      ttCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATaa                                                                                               <  2:148299/149‑1 (MQ=255)
                                                           tGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGCCATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCa                                                                                            >  3:97658/1‑147 (MQ=255)
                                                           tGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCa                                                                                            >  8:121370/1‑147 (MQ=255)
                                                            gCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCa                                                                                                                                                <  4:345868/94‑1 (MQ=255)
                                                            gCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCa                                                                                                                                                >  3:345868/1‑94 (MQ=255)
                                                            gCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGt                                                                                        >  5:365567/1‑149 (MQ=255)
                                                            gCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGt                                                                                        >  5:365570/1‑149 (MQ=255)
                                                              atatatGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCGAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGGCGCCTGCTGAATGAGATAACCAAAGtt                                                                                       >  5:19060/1‑149 (MQ=255)
                                                               tatatGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATg                                                                                                                                 >  7:331278/1‑106 (MQ=255)
                                                               tatatGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATg                                                                                                                                 <  8:331278/106‑1 (MQ=255)
                                                               tatatGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGttt                                                                                      >  5:340405/1‑149 (MQ=255)
                                                                     gCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgc                                                                                <  1:323651/149‑1 (MQ=255)
                                                                       gACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGAGATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgcgt                                                                              <  1:110118/149‑1 (MQ=255)
                                                                        aCGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgcgtg                                                                             >  2:407218/1‑149 (MQ=255)
                                                                          gcggcgGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgcgtgcg                                                                           >  3:199293/1‑149 (MQ=255)
                                                                                gAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGa                                                                       >  1:37798/1‑147 (MQ=255)
                                                                                gAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGa                                                                       >  1:301674/1‑147 (MQ=255)
                                                                                gAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATa                                                                     <  6:175547/149‑1 (MQ=255)
                                                                                 aaaGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAg                                                                    <  1:407218/149‑1 (MQ=255)
                                                                                   aGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcg                                                                  >  6:98977/1‑149 (MQ=255)
                                                                                    gggCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcgt                                                                 >  6:353366/1‑149 (MQ=255)
                                                                                      gCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcgttc                                                               <  5:98977/149‑1 (MQ=255)
                                                                                       cGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcgttcg                                                              <  1:5602/149‑1 (MQ=255)
                                                                                        gATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcgttcgt                                                             <  1:116103/149‑1 (MQ=255)
                                                                                            aaGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATc                                                         <  4:366299/149‑1 (MQ=255)
                                                                                             aGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCg                                                        >  6:105603/1‑149 (MQ=255)
                                                                                                    cacaAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGCCATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAcct                                                 >  2:269618/1‑149 (MQ=255)
                                                                                                    cacaAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGATGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCAtcact                                                 >  7:206096/1‑145 (MQ=255)
                                                                                                     acaAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATc                                                         <  3:14687/140‑1 (MQ=255)
                                                                                                     acaAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATc                                                         >  4:14687/1‑140 (MQ=255)
                                                                                                      caAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAcctcc                                               >  8:401258/1‑149 (MQ=255)
                                                                                                       aaaTCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAcctcct                                              >  1:399484/1‑149 (MQ=255)
                                                                                                        aaTCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGa                                            <  6:365570/149‑1 (MQ=255)
                                                                                                        aaTCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGa                                            <  6:365567/149‑1 (MQ=255)
                                                                                                         aTCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGa                                            <  5:316297/149‑1 (MQ=255)
                                                                                                         aTCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGa                                            >  7:85720/1‑149 (MQ=255)
                                                                                                            gTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgcg                                                                               >  8:313929/1‑111 (MQ=255)
                                                                                                            gTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgcg                                                                               <  7:313929/111‑1 (MQ=255)
                                                                                                               aCAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCATCATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGc                                      <  6:99372/149‑1 (MQ=255)
                                                                                                                 aaGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCa                                    >  4:98739/1‑149 (MQ=255)
                                                                                                                          agcaAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGaaa                           <  7:121370/149‑1 (MQ=255)
                                                                                                                             aaataaatGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACgg                        >  1:18581/1‑149 (MQ=255)
                                                                                                                               ataaatGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgt                      <  6:266185/149‑1 (MQ=255)
                                                                                                                               ataaatGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgt                      >  4:44974/1‑149 (MQ=255)
                                                                                                                               ataaatCGCGTCATCGCAACTGGCGCAGCCAACATGGAAATGCTCCCCAAAAACGGTCGCCTGCTGAATGAGATAACCAATGTTTTCCACCGTGAGTTGATAGCGTTCGTTATCGCATCCCTCCTGACGCTGCCAGTTCTGATCCGGgt                      >  6:46989/1‑149 (MQ=255)
                                                                                                                                taaatGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCTGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgtg                     <  6:340405/149‑1 (MQ=255)
                                                                                                                                  aatGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgtgc                    <  2:37798/148‑1 (MQ=255)
                                                                                                                                      gCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgtgcgtgc                <  6:309519/148‑1 (MQ=255)
                                                                                                                                         tCATCGCCACTGGCGCAGCCAACATGGCCATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATACCCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgt                      <  5:280174/139‑1 (MQ=255)
                                                                                                                                         tCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCATCGTGCGTTGATAGCGTTCGTGATCGCATACCTCCTGACGCTGCCAGTTCTGAAACGGGTGCGTGCCAcc            >  6:378786/1‑149 (MQ=255)
                                                                                                                                         tCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAcc            >  2:250495/1‑149 (MQ=255)
                                                                                                                                         tCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTGTTCCAGCGTGCGTTGCTAGCGGTCGTTATCGCAAACCTCCTGACGCTGCCAGTTCTGCAACGGgggcgctcgg              >  6:280174/1‑142 (MQ=255)
                                                                                                                                           aTCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgc          >  3:95180/1‑149 (MQ=255)
                                                                                                                                           aTCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgc          >  1:312366/1‑149 (MQ=255)
                                                                                                                                            tCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgcc         >  1:283081/1‑149 (MQ=255)
                                                                                                                                            tCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgcc         >  6:351214/1‑149 (MQ=255)
                                                                                                                                                  cTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgtgcgtgc                <  3:281991/136‑1 (MQ=255)
                                                                                                                                                  cTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgtgcgtgc                >  4:281991/1‑136 (MQ=255)
                                                                                                                                                   tGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAAtt  >  7:367528/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATT  >  NC_000913/606810‑607105

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