Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,217,392 T→C V95A (GTA→GCA)  aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,217,3920TC98.3% 390.4 / NA 118V95A (GTA→GCA) aceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (53/63);  minor base A (1/1);  total (54/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 121 total aligned reads displayed.

GCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCG  >  NC_000913/4217244‑4217538
                                                                                                                                                    |                                                                                                                                                  
 cgcAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGCGTCGAAAAAAGGCTACATCCACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAAGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAgc                                                                                                                                                 >  2:1067338/1‑151 (MQ=255)
 cgcAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGgcaggcag                                                                                                                                                  >  2:1337065/1‑150 (MQ=255)
  gcAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAgcg                                                                                                                                                >  2:1176885/1‑151 (MQ=255)
   cAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGCAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAgcgg                                                                                                                                               >  2:965260/1‑151 (MQ=255)
    aaCTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAgcggc                                                                                                                                              <  3:47956/151‑1 (MQ=255)
    aaCTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAgcggc                                                                                                                                              <  2:1235077/151‑1 (MQ=255)
        gggCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGa                                                                                                                                           <  2:420748/150‑1 (MQ=255)
         ggCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACg                                                                                                                                         <  1:846362/151‑1 (MQ=255)
          gcgcAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGc                                                                                                                                        <  1:881510/151‑1 (MQ=255)
           cgcAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCt                                                                                                                                       <  1:435311/151‑1 (MQ=255)
            gcAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTa                                                                                                                                      <  1:104069/151‑1 (MQ=255)
             cAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTaa                                                                                                                                     <  4:510608/151‑1 (MQ=255)
             cAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTaa                                                                                                                                     <  2:11955/151‑1 (MQ=255)
              aGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGTCGCTaa                                                                                                                                     >  2:1026197/1‑150 (MQ=255)
               gCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAAcc                                                                                                                                   >  3:127399/1‑151 (MQ=255)
                cGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCt                                                                                                                                  <  2:1009462/151‑1 (MQ=255)
                 gAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTg                                                                                                                                 <  1:379216/151‑1 (MQ=255)
                 gAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGGCTATCTGTCGGGATGGCAGGCAGCGGCGGGCGCTAACCt                                                                                                                                  >  1:282737/1‑150 (MQ=255)
                  aaaaTGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTgg                                                                                                                                <  2:464045/151‑1 (MQ=255)
                   aaaTGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTggc                                                                                                                               <  2:1737607/151‑1 (MQ=255)
                     aTGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTggcgg                                                                                                                             <  1:971309/151‑1 (MQ=255)
                     aTGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTggcgg                                                                                                                             >  4:118876/1‑151 (MQ=255)
                         ggCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCCGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAAACTGGCGGcc                                                                                                                           >  2:439061/1‑149 (MQ=255)
                         ggCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGGAAGCGGGTATTGAAGCAGACTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGcc                                                                                                                           >  2:347482/1‑149 (MQ=255)
                         ggCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCa                                                                                                                          <  1:539047/150‑1 (MQ=255)
                           cGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCa                                                                                                                       <  2:995793/151‑1 (MQ=255)
                             tCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGAGGGTATTCAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCAt                                                                                                                      >  1:1161967/1‑150 (MQ=255)
                             tCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAAGCGCTGCAACCCGCCAAAGCGGGGATTGAAGCAGTCTATTTGTCGGGGTGGCAGGCAGCGGCCGACGCTAACCTGGCGGCCCGCGTg                                                                                                                     >  1:78854/1‑151 (MQ=255)
                              ctgctgCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATg                                                                                                                     >  2:1451971/1‑150 (MQ=255)
                              ctgctgCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCAATGACTGGCGGTCAGGCGCTGCCACCCGGGAAAGCGGGGATTAAAGCAGTCCATCTGTCGGGAGGGACGGCAGCGGCGGACCCTAACCATGTGGCCAGcgt                                                                                                                      >  1:19900/1‑147 (MQ=255)
                                   gcACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCg                                                                                                               <  1:1062821/151‑1 (MQ=255)
                                   gcACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCg                                                                                                               >  3:218331/1‑151 (MQ=255)
                                   gcACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCg                                                                                                               <  4:396979/151‑1 (MQ=255)
                                    cACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCgg                                                                                                              <  2:537548/151‑1 (MQ=255)
                                     aCGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCgg                                                                                                              >  2:1667935/1‑150 (MQ=255)
                                     aCGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCgg                                                                                                              >  2:1050590/1‑150 (MQ=255)
                                     aCGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGACGCTGCAACAGGAGAAAGCGGGAATTGAAGCAGTCTATCTGTCGGGATGGCAAGACGCGGCGGCAGCTAAACTGGCGGGCAGCCAGTATCCg                                                                                                               >  1:1461106/1‑149 (MQ=255)
                                       ggTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCGGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGAt                                                                                                            <  1:262459/150‑1 (MQ=255)
                                        gTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCa                                                                                                          <  1:1604138/151‑1 (MQ=255)
                                        gTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCa                                                                                                          <  1:1628650/151‑1 (MQ=255)
                                           aGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGGGAAAGCGGGTATTGAAGCCGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTc                                                                                                       >  1:505303/1‑151 (MQ=255)
                                               gAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGGATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGctc                                                                                                   >  2:1187272/1‑151 (MQ=255)
                                                aaaaaaGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAAGCGCTGCAACAGGAGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCCGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATAAGTCGct                                                                                                    >  1:301038/1‑149 (MQ=255)
                                                 aaaaaGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGctct                                                                                                  >  2:453595/1‑150 (MQ=255)
                                                  aaaaGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTAt                                                                                                >  2:1673769/1‑151 (MQ=255)
                                                     aGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCg                                                                                             <  1:1469212/151‑1 (MQ=255)
                                                     aGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCg                                                                                             <  2:213186/151‑1 (MQ=255)
                                                      ggCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCgg                                                                                            <  1:1300120/151‑1 (MQ=255)
                                                      ggCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGGATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCgg                                                                                            >  1:807404/1‑151 (MQ=255)
                                                      ggCTACATCAACAGCCTCGGCGCACTGACTGGCGGACAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCgg                                                                                            >  3:376053/1‑151 (MQ=255)
                                                          aCATCATCAGCCTCGGCGTGCTGACTGGCGGTCTGTGTCTGCAACAGGGGAAAGCGGATATTGAAGCAGGCAAGCTGTCGGTAGGGCAGGCCGCGGCCGGCGCTAACCTGGCGGCCAGCATGTATACGTATCAGACCCGCTTTCCGTCaaa                                                                                        <  1:897918/151‑1 (MQ=255)
                                                           cATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTTCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCCGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAAc                                                                                       >  4:122438/1‑151 (MQ=255)
                                                            aTCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCCAAc                                                                                       >  1:1544626/1‑150 (MQ=255)
                                                            aTCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACt                                                                                      <  3:38072/151‑1 (MQ=255)
                                                              cAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTc                                                                                     >  1:1709928/1‑150 (MQ=255)
                                                                aCAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGt                                                                                  <  4:448209/151‑1 (MQ=255)
                                                                    ccTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGcc                                                                               <  1:1475596/150‑1 (MQ=255)
                                                                    ccTCGGCGCACTGACTGACGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCa                                                                              <  4:313372/151‑1 (MQ=255)
                                                                      tCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTTGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGc                                                                            >  1:1221591/1‑151 (MQ=255)
                                                                      tCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGc                                                                            >  2:1539916/1‑151 (MQ=255)
                                                                      tCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGc                                                                            >  2:652537/1‑151 (MQ=255)
                                                                       cGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCt                                                                           <  3:446204/151‑1 (MQ=255)
                                                                         gcgcACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCtgt                                                                         <  3:307502/151‑1 (MQ=255)
                                                                          cgcACTGACTGGCGGTCGGGCGCTGCAACAGGCGAAAGCGCGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGGCGCTAACCTGGCGGCCAGTATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCa                                                                              <  1:1018407/145‑1 (MQ=255)
                                                                          cgcACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTCGCGGCCAGCATGTATCCGGATCCGTCGCTCTATCCGGCAAACTCGGTGCCAGCtgtg                                                                        >  1:253303/1‑151 (MQ=255)
                                                                            cACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGGATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCACGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTgtggg                                                                      >  1:1677930/1‑150 (MQ=255)
                                                                             actgactgGCGGTCAGGCGCTGCAACAGGCGAAAACGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTgtggtg                                                                     <  2:777558/151‑1 (MQ=255)
                                                                                    ggCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCgg                                                               <  2:1434828/150‑1 (MQ=255)
                                                                                     gCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGAt                                                             <  2:1732399/151‑1 (MQ=255)
                                                                                      cGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATc                                                            <  4:290858/151‑1 (MQ=255)
                                                                                              cgcTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAacac                                                    >  1:1648310/1‑151 (MQ=255)
                                                                                               gcTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGAAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAacac                                                    <  1:957272/150‑1 (MQ=255)
                                                                                                     aCAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCcgt                                             <  2:1115208/151‑1 (MQ=255)
                                                                                                      cAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCcgtc                                            >  2:405099/1‑151 (MQ=255)
                                                                                                       aGGCGAAAGCGGGTATTGGAGCAGTCTATGTGGCGGGATGTGTGGAAGCGGCGGAGGCTAGCCTGGCGGGCAGCATGTATCCGGTTCAGTCGCTCTATCCGGCAAACGCGTTGCCAGCTGTGGTGGAGCGGatcaccatcaccagc                                                <  1:560166/146‑4 (MQ=255)
                                                                                                       aGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCcgtcg                                           >  2:470462/1‑151 (MQ=255)
                                                                                                       aGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCcgtcg                                           >  3:179575/1‑151 (MQ=255)
                                                                                                           gAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGcc                                       >  4:346455/1‑151 (MQ=255)
                                                                                                             aaGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGa                                     <  2:790065/151‑1 (MQ=255)
                                                                                                               gCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCg                                      >  1:248984/1‑148 (MQ=255)
                                                                                                               gCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGAt                                    <  4:388950/150‑1 (MQ=255)
                                                                                                               gCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATc                                   >  1:520203/1‑151 (MQ=255)
                                                                                                                 gggTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGTGCGGATCAACAACCCCTTCCGTCGTGCCGATCAg                                 <  1:1736429/151‑1 (MQ=255)
                                                                                                                   gTATTGAAGCAGTCTATCTGTCGGGATGGCCGGCCGCGGCGGACGCTAACCTGGCTGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCTGCTGTGGTGGAGCGGATCAACAACCCCTTCCGTCGTGCCGATCAGAt                               <  1:1220408/151‑1 (MQ=255)
                                                                                                                   gTATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATc                                   >  2:72551/1‑147 (MQ=255)
                                                                                                                    tATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATc                                   >  2:307762/1‑146 (MQ=255)
                                                                                                                    tATTGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATc                                   >  3:450006/1‑146 (MQ=255)
                                                                                                                       tGAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCa                            <  2:121594/150‑1 (MQ=255)
                                                                                                                        gAAGCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACTCCCTCCGTCGTGCCGATCAGATCCa                            <  1:1295276/149‑1 (MQ=255)
                                                                                                                           gCAGTCTATCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGt                       <  1:534661/151‑1 (MQ=255)
                                                                                                                                  aTCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcgc                   <  2:1392545/148‑1 (MQ=255)
                                                                                                                                  aTCTGTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCgg                 <  2:1067403/150‑1 (MQ=255)
                                                                                                                                      gTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACGCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCa              <  3:248366/149‑1 (MQ=255)
                                                                                                                                      gTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGc               <  2:1475589/148‑1 (MQ=255)
                                                                                                                                      gTCGGGATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCAt             <  3:127864/150‑1 (MQ=255)
                                                                                                                                         gggATGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAg         <  1:1654661/151‑1 (MQ=255)
                                                                                                                                            aTGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCg      <  1:1360436/151‑1 (MQ=255)
                                                                                                                                             tGGCAGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCg      >  2:1047218/1‑150 (MQ=255)
                                                                                                                                             tGGCAGGCAGCGGCGGAAGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCgg     >  1:1053436/1‑151 (MQ=255)
                                                                                                                                                 aGGCAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCg  >  3:289199/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCG  >  NC_000913/4217244‑4217538

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