Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,874 T→C 5.1% Y346H (TAC→CAC)  ttuD → hydroxypyruvate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,8740TC 5.1% 235.3 / 2.2 78Y346H (TAC→CAC) ttuDhydroxypyruvate reductase
Reads supporting (aligned to +/- strand):  ref base T (32/42);  new base C (2/2);  total (34/44)
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.88e-01

TTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCG  >  NC_002947/4891728‑4892021
                                                                                                                                                  |                                                                                                                                                     
ttCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTAc                                                                                                                                                     >  7:240136/1‑149 (MQ=255)
     cacGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCt                                                                                                                                                >  1:303539/1‑149 (MQ=255)
     cacGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCt                                                                                                                                                >  2:362696/1‑149 (MQ=255)
          ccAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCg                                                                                                                                           <  2:303539/149‑1 (MQ=255)
           cAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCgg                                                                                                                                          <  3:220979/149‑1 (MQ=255)
           cAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCgg                                                                                                                                          <  7:13637/149‑1 (MQ=255)
            aGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCGCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGt                                                                                                                                         <  1:320988/149‑1 (MQ=255)
               cGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGa                                                                                                                                       >  7:380089/1‑148 (MQ=255)
               cGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGa                                                                                                                                       >  7:322248/1‑148 (MQ=255)
                 cTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGacac                                                                                                                                    >  7:85356/1‑149 (MQ=255)
                    aaaGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGa                                                                                                                                 <  2:291683/149‑1 (MQ=255)
                    aaaGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTACTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGa                                                                                                                                 >  3:163647/1‑149 (MQ=255)
                     aaGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTgg                                                                                                                                              >  3:344481/1‑135 (MQ=255)
                     aaGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTgg                                                                                                                                              <  4:344481/135‑1 (MQ=255)
                      aGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACg                                                                                                                               <  7:348826/149‑1 (MQ=255)
                       gcgcCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACgg                                                                                                                              <  6:108597/149‑1 (MQ=255)
                       gcgcCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACgg                                                                                                                              <  4:70707/149‑1 (MQ=255)
                          cccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTACACGcccggcag                                                                                                                                            >  7:81529/1‑127 (MQ=255)
                          cccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTACACGcccggcag                                                                                                                                            <  8:81529/132‑6 (MQ=255)
                          cccTGCGTGATCCTGTCGGGTGGCGAAACAACCGTGACCGTGCGCGGCAATGGACGTGGCGGGCGTAAAGCCTAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCAt                                                                                                                           >  8:6677/1‑149 (MQ=255)
                               cGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGa                                                                                                                        <  2:24096/147‑1 (MQ=255)
                               cGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGa                                                                                                                        <  8:113904/147‑1 (MQ=255)
                               cGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGAt                                                                                                                       <  2:283923/148‑1 (MQ=255)
                               cGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGAt                                                                                                                       <  1:362696/148‑1 (MQ=255)
                                 tGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGAt                                                                                                                       >  7:63986/1‑146 (MQ=255)
                                 tGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGAt                                                                                                                       <  8:63986/146‑1 (MQ=255)
                                 tGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGc                                                                                                                    >  5:94777/1‑149 (MQ=255)
                                         tCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGgaa                                                                                                              <  2:402996/147‑1 (MQ=255)
                                             gTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAAc                                                                                                        >  3:371924/1‑149 (MQ=255)
                                              tGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACg                                                                                                       >  4:234226/1‑149 (MQ=255)
                                                 cGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCg                                                                                                    <  2:118959/149‑1 (MQ=255)
                                                       caccGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACg                                                                                                                                                    <  8:119075/95‑1 (MQ=255)
                                                       caccGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACg                                                                                                                                                    >  7:119075/1‑95 (MQ=255)
                                                           gtgaccgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCAt                                                                                          <  2:90172/149‑1 (MQ=255)
                                                           gtgaccgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCAt                                                                                          <  5:214666/149‑1 (MQ=255)
                                                                cgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGAccc                                                                                     <  8:380089/149‑1 (MQ=255)
                                                                cgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGAccc                                                                                     >  6:204814/1‑149 (MQ=255)
                                                                cgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGAccc                                                                                     >  7:375399/1‑149 (MQ=255)
                                                                 gtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGAcccc                                                                                    <  5:18606/149‑1 (MQ=255)
                                                                    cgcgGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCgg                                                                                  <  8:85356/148‑1 (MQ=255)
                                                                    cgcgGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGc                                                                                 <  4:260591/149‑1 (MQ=255)
                                                                       ggCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAg                                                                             >  1:380511/1‑149 (MQ=255)
                                                                         cAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGc                                                                            >  4:120943/1‑148 (MQ=255)
                                                                         cAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCt                                                                           <  1:59625/149‑1 (MQ=255)
                                                                            tGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGga                                                                                                               <  3:316071/111‑1 (MQ=255)
                                                                            tGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGga                                                                                                               >  4:316071/1‑111 (MQ=255)
                                                                              gCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGc                                                                       <  5:299205/148‑1 (MQ=255)
                                                                               ccGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCa                                                                     <  1:174920/149‑1 (MQ=255)
                                                                               ccGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCa                                                                     <  8:375399/149‑1 (MQ=255)
                                                                                    gcgggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAgcgc                                                                <  4:371924/148‑1 (MQ=255)
                                                                                      gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAgcgcggc                                                             >  5:339566/1‑146 (MQ=255)
                                                                                      gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGcc                                                               >  5:135956/1‑147 (MQ=255)
                                                                                      gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGcc                                                               <  6:135956/147‑1 (MQ=255)
                                                                                        gcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAg                                                           >  6:173072/1‑149 (MQ=255)
                                                                                           tAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCt                                                        >  1:117558/1‑149 (MQ=255)
                                                                                              cGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTgg                                                     <  7:110350/149‑1 (MQ=255)
                                                                                                ccGAGTTCCTGCTCAGCCTCCCCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                                   <  8:322248/149‑1 (MQ=255)
                                                                                                ccGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                                   >  3:288242/1‑149 (MQ=255)
                                                                                                  gAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCTGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCt                                                 >  4:323154/1‑149 (MQ=255)
                                                                                                       ccTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCg                                             <  2:380511/148‑1 (MQ=255)
                                                                                                        cTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                                   <  6:251015/141‑1 (MQ=255)
                                                                                                        cTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                                   >  5:251015/1‑141 (MQ=255)
                                                                                                           cTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCa                                         <  7:6677/148‑1 (MQ=255)
                                                                                                            tCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGc                                       >  5:83954/1‑149 (MQ=255)
                                                                                                               gCCTCACCGAACGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgac                                    <  1:189175/149‑1 (MQ=255)
                                                                                                                    aCCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCt                               >  6:201809/1‑149 (MQ=255)
                                                                                                                     ccGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                              >  6:270696/1‑149 (MQ=255)
                                                                                                                        aaaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGac                           <  2:117558/149‑1 (MQ=255)
                                                                                                                         aaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACTTTTTCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGac                           <  1:334262/148‑1 (MQ=255)
                                                                                                                          aGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCTGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGacaa                         <  3:323154/149‑1 (MQ=255)
                                                                                                                             cTGAAAGGCCTGCCGGGCGTGTTCGCCCTGGCCGGGGACACCGACGGCATCGATGGCGCGGAAGAAAACGCCGGTGCCTTCATTACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCg                                             <  1:160759/126‑1 (MQ=255)
                                                                                                                              tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCAACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGacaacaa                      >  5:15975/1‑148 (MQ=255)
                                                                                                                               gAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGacaacaac                     >  3:221686/1‑148 (MQ=255)
                                                                                                                                 aaGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACaacaacaac                  <  8:240136/149‑1 (MQ=255)
                                                                                                                                      ggccagggCGTACACGCCCGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggctac            >  5:118754/6‑149 (MQ=255)
                                                                                                                                         ccGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACaacaacaac                  <  4:371084/141‑1 (MQ=255)
                                                                                                                                         ccGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACaacaacaac                  >  3:371084/1‑141 (MQ=255)
                                                                                                                                         cagggCGTACACGCCCGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggctacggc         <  6:118754/147‑1 (MQ=255)
                                                                                                                                             tcttgtACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACggctacggctac      <  7:108248/146‑1 (MQ=255)
                                                                                                                                               gtgtACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGATGCCTTCATGACCCCGGCGAGGCTACGCCATCGGGCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACggctacggctac      >  1:125302/1‑149 (MQ=255)
                                                                                                                                                tgtACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTc   >  5:94734/1‑149 (MQ=255)
                                                                                                                                                  tACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTcg  <  5:173072/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCG  >  NC_002947/4891728‑4892021

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: