Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,580,798 T→G 5.2% I388S (ATC→AGC) ‡ tktB → transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,580,7980TG 5.2% 277.0 / 4.2 89I388S (ATC→AGC) ‡tktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (44/40);  new base G (3/2);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCAGAAATATATCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTG  >  NC_000913/2580653‑2580944
                                                                                                                                                 |                                                                                                                                                  
aCTCAGAAATATATCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTg                                                                                                                                                 <  5:220050/149‑1 (MQ=255)
  tCAGAAATATATCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTgg                                                                                                                                                >  4:323969/1‑148 (MQ=255)
  tCAGAAATATATCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCAGCTgg                                                                                                                                                >  8:178484/1‑148 (MQ=255)
    aGAAATATATCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGa                                                                                                                                               >  3:268015/1‑147 (MQ=255)
    aGAAATATATCAATGAGTTACAGGCAAATCCGGAGAAAATCGCTACCCGTAAGGCTTCGCAAAATACTCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGa                                                                                                                                               >  8:190494/1‑147 (MQ=255)
      aaaTATATCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAagg                                                                                                                                           >  8:176692‑M2/1‑146 (MQ=255)
           tatCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCt                                                                                                                                       <  5:308569/148‑1 (MQ=255)
           tatCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTg                                                                                                                                      >  8:379339/1‑149 (MQ=255)
             tCAATGAGTTCCAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGtt                                                                                                                                    <  7:178484/149‑1 (MQ=255)
               aaTGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGttt                                                                                                                                   <  7:379339/148‑1 (MQ=255)
               aaTGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTc                                                                                                                                  >  7:360494/1‑149 (MQ=255)
               aaTGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTc                                                                                                                                  <  7:56068/149‑1 (MQ=255)
                aTGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCg                                                                                                                                 <  8:227941/149‑1 (MQ=255)
                         aGGCAAATCCGGCTAAAATCTCAACCCGTAAGGCTTCGCAAAATACTCTTAACGCTTACGTGCCGATGCTGCCTGAGTTGCTCGGCGGTTAGGCGGATCTGGCTCCCAGCAACATGAGCATATGGAAAGGTTATGTTTCGATgaagg                                                                                                                          >  4:338352/1‑147 (MQ=255)
                           gCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTACCGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCACCCTGCCCATCTGGAAAGGTTCTGTTTCGCTGAAGGACGa                                                                                                                      >  1:382485/1‑149 (MQ=255)
                           gCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTgaagg                                                                                                                          >  7:28679/1‑145 (MQ=255)
                            cAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGAt                                                                                                                     <  3:323969/149‑1 (MQ=255)
                                  cGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGCTGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCACCCTGACCATCTGGAAAGGGTCTGTTTCGCTGAAGGAAGCTCCAGCg                                                                                                               >  7:80655/1‑149 (MQ=255)
                                   ggCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATc                                                                                                                    >  6:101558/1‑143 (MQ=255)
                                   ggCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATc                                                                                                                    <  5:101558/143‑1 (MQ=255)
                                     cGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTCCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTGCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGc                                                                                                            <  2:369191/149‑1 (MQ=255)
                                      gAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGc                                                                                                            >  6:367888/1‑148 (MQ=255)
                                         aaTCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTTTCCAGATGGTCaggttgctgggagc                                                                                                                                >  8:10804‑M2/1‑111 (MQ=255)
                                         aaTCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTTTCCAGATGGTCaggttgctgggagc                                                                                                                                <  7:10804‑M2/125‑15 (MQ=255)
                                         aaTCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACt                                                                                                        <  3:296961/149‑1 (MQ=255)
                                             gCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACa                                                                                                     <  7:380333/148‑1 (MQ=255)
                                                aCCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCa                                                                                                 >  6:176204/1‑149 (MQ=255)
                                                   cGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTa                                                                                              <  7:190494/149‑1 (MQ=255)
                                                   cGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTa                                                                                              <  1:394665/149‑1 (MQ=255)
                                                     tAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACg                                                                                            >  4:309011/1‑149 (MQ=255)
                                                       aGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACggg                                                                                          <  8:226136/149‑1 (MQ=255)
                                                           ttCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGACTGAGTTGCTCGGCGGTTCGGCGGAACTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGCAGGAAGATCCAGCGGGCACCTACATTCACTACgtgg                                                                                         >  4:48377/1‑146 (MQ=255)
                                                           ttCGCAAAATACGCGTCACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTGCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACgggg                                                                                         <  3:48377/146‑1 (MQ=255)
                                                             cGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGt                                                                                    >  7:237385/1‑149 (MQ=255)
                                                                aaaaTACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTTTCCAGATGGTCaggttgctgggagccagatccgc                                                                                                                       >  5:370837‑M2/1‑88 (MQ=255)
                                                                aaaaTACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTTTCCAGATGGTCaggttgctgggagccagatccgc                                                                                                                       <  6:370837‑M2/111‑24 (MQ=255)
                                                                   aTACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGTTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAAttt                                                                              <  6:239870/149‑1 (MQ=255)
                                                                    tACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTg                                                                             <  5:176204/149‑1 (MQ=255)
                                                                      cgcttaccgcttacGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCACAGAAACCTGAACATCTGGAAATGTTCTGTTTCGCTGAAGGAAGATCCAGCGTGAAACTACTTTCACTACGGTGTTCGTGAATTTgtc                                                                           <  4:385069/149‑3 (MQ=255)
                                                                      cgcttaacgcttaCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTgg                                                                            <  8:331537/148‑1 (MQ=255)
                                                                          taacgcttaCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGa                                                                       >  4:128998/1‑149 (MQ=255)
                                                                           aacgcgtaCGGGCCGATGCGGCCTGAGGTGCTCGGCGGTGCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGa                                                                       <  8:360494/148‑1 (MQ=255)
                                                                              gcttaCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCgg                                                                                                              <  8:43694/106‑1 (MQ=255)
                                                                              gcttaCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCgg                                                                                                              >  7:43694/1‑106 (MQ=255)
                                                                              gcttaCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTg                                                                             >  2:350290/1‑139 (MQ=255)
                                                                              gcttaCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTg                                                                             <  1:350290/139‑1 (MQ=255)
                                                                               cttaCGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCACCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGc                                                                   >  7:129882/1‑148 (MQ=255)
                                                                                           tgctgcCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACgg                                                       >  7:45156/1‑148 (MQ=255)
                                                                                              tgcCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTa                                                                                              <  5:27076/106‑1 (MQ=255)
                                                                                              tgcCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTa                                                                                              >  6:27076/1‑106 (MQ=255)
                                                                                               gcCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATcg                                                  >  3:284837/1‑149 (MQ=255)
                                                                                                  tGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATcgc                                                 >  6:243962/1‑147 (MQ=255)
                                                                                                  tGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATcgc                                                 <  5:243962/147‑1 (MQ=255)
                                                                                                   gAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCATCGGGCAACTACATTCACTACTCGGTGCGTGAATTTGGCATGACAGCTATCGCCAACGGCATCGCGca                                              >  5:162382/1‑149 (MQ=255)
                                                                                                        gCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAcg                                         >  4:279898/1‑149 (MQ=255)
                                                                                                            ggcggTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAcggcgg                                     <  8:45156/149‑1 (MQ=255)
                                                                                                              cggTTCGGCGGATCTGGCTCCCAGCAACATTACCATCTGGAAAGGTTCTGTTTCGCTGAAGTAATATCCAGCGGGCAACTACATTCACTACGGGGTGCGTTAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAcggcg                                      <  3:338352/146‑1 (MQ=255)
                                                                                                                agctggcggATCTGGCTCCCAGCAAGCTGACCATCTGGAAAGGGTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGGGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCttt                                 <  5:133179/145‑1 (MQ=255)
                                                                                                                     gcggATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGcc                            >  6:213478/1‑149 (MQ=255)
                                                                                                                     gcggATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGcc                            <  8:251648/149‑1 (MQ=255)
                                                                                                                     gcggATCTGGCTCCCAGCAACCTGACCATCAGGAAAGGTTCTGTTACTCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGCAGCCTGAATTTGGCATGACCACTATCGCCAAACGCATCGCGCACCACCGCGGCTTtgaga                             >  6:133179/1‑145 (MQ=255)
                                                                                                                      cggATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCg                           <  8:237385/149‑1 (MQ=255)
                                                                                                                       ggATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCg                           >  7:115542/1‑148 (MQ=255)
                                                                                                                            tGGCTCCCAGCAACCTGCCCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGAACCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCCCCACGGCGGCTTTGTGCCGAATAcc                     >  1:107850/1‑149 (MQ=255)
                                                                                                                            tGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATAcc                     >  4:317012/1‑149 (MQ=255)
                                                                                                                             ggCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAAttt                                                                              >  1:73229/1‑91 (MQ=255)
                                                                                                                             ggCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAAttt                                                                              <  2:73229/91‑1 (MQ=255)
                                                                                                                             ggCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACcg                    >  1:280835/1‑149 (MQ=255)
                                                                                                                              gCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACcgc                   <  3:317012/149‑1 (MQ=255)
                                                                                                                                  ccAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACg               >  5:187774/1‑149 (MQ=255)
                                                                                                                                   cAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGt              >  4:140369/1‑149 (MQ=255)
                                                                                                                                   cAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGt              >  7:137226/1‑149 (MQ=255)
                                                                                                                                    aGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGtt             >  6:292833/1‑149 (MQ=255)
                                                                                                                                     gCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTc            <  5:292833/149‑1 (MQ=255)
                                                                                                                                        aCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTGTCGCCAACGGCa                                                     <  4:174236/105‑1 (MQ=255)
                                                                                                                                        aCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTGTCGCCAACGGCa                                                     >  3:174236/1‑105 (MQ=255)
                                                                                                                                          cTGACCATCTGGAACGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAt       <  2:377147/149‑1 (MQ=255)
                                                                                                                                           tGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGt                                                                                    <  3:137098/71‑1 (MQ=255)
                                                                                                                                           tGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGt                                                                                    >  4:137098/1‑71 (MQ=255)
                                                                                                                                           tGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtg      >  4:264192/1‑149 (MQ=255)
                                                                                                                                            gACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgt     >  3:75427/1‑149 (MQ=255)
                                                                                                                                             aCCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgtt    <  3:279898/149‑1 (MQ=255)
                                                                                                                                              ccATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgttt   <  3:264192/149‑1 (MQ=255)
                                                                                                                                              ccATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgttt   <  7:134454/149‑1 (MQ=255)
                                                                                                                                               cATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCttt                                 >  5:407341/1‑118 (MQ=255)
                                                                                                                                               cATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCttt                                 <  6:407341/118‑1 (MQ=255)
                                                                                                                                               cATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgtttg  >  5:233143/1‑149 (MQ=255)
                                                                                                                                                aTCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGa        >  7:382208/1‑142 (MQ=255)
                                                                                                                                                aTCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGa        <  8:382208/142‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACTCAGAAATATATCAATGAGTTACAGGCAAATCCGGCGAAAATCGCTACCCGTAAGGCTTCGCAAAATACGCTTAACGCTTACGGGCCGATGCTGCCTGAGTTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTG  >  NC_000913/2580653‑2580944

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