Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,921,101 T→C 14.4% M334V (ATG→GTG)  gudP ← putative D‑glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,921,1010TC14.4% 2291.3 / 314.7 875M334V (ATG→GTG) gudPputative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (448/300);  new base C (80/46);  total (528/347)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

Only 100 of 882 total aligned reads displayed.

CCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATC  >  NC_000913/2920953‑2921249
                                                                                                                                                    |                                                                                                                                                    
ccGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACAt                                                                                                                                                      <  1:753415/149‑1 (MQ=255)
    tGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGccc                                                                                                                                                  >  4:7079573/1‑149 (MQ=255)
    tGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGccc                                                                                                                                                  <  3:5120003/149‑1 (MQ=255)
      aTCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAATCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCAt                                                                                                                                                >  3:3473997/1‑149 (MQ=255)
      aTCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCAt                                                                                                                                                >  1:2177414/1‑149 (MQ=255)
      aTCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCAt                                                                                                                                                >  1:4054048/1‑149 (MQ=255)
      aTCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACACGCCCAt                                                                                                                                                >  2:1373668/1‑149 (MQ=255)
       tctcTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCAtt                                                                                                                                               <  4:5616816/149‑1 (MQ=255)
         tctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTAc                                                                                                                                             >  4:336032/1‑149 (MQ=255)
         tctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTAc                                                                                                                                             >  1:3666265/1‑149 (MQ=255)
         tctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTAc                                                                                                                                             >  4:5253201/1‑149 (MQ=255)
          ctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACACGCCCATTACg                                                                                                                                            >  1:238905/1‑149 (MQ=255)
           ttttGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGa                                                                                                                                           <  2:3967973/149‑1 (MQ=255)
           ttttGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGa                                                                                                                                           >  4:1371343/1‑149 (MQ=255)
           ttttGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACACGCCCATTACGa                                                                                                                                           <  4:2313471/149‑1 (MQ=255)
             ttGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATc                                                                                                                                         <  2:458208/149‑1 (MQ=255)
               ggCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCgg                                                                                                                                       >  3:5797426/1‑149 (MQ=255)
                    ggTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGttt                                                                                                                                  <  1:4441341/149‑1 (MQ=255)
                    ggTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGttt                                                                                                                                  <  3:891792/149‑1 (MQ=255)
                           gCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACACGCCCATTACGATCGGTGTTTTACgcgc                                                                                                                           >  2:1578946/1‑149 (MQ=255)
                            ccATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCTGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCa                                                                                                                          <  3:4353682/149‑1 (MQ=255)
                              aTTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCGAt                                                                                                                        >  3:3179075/1‑149 (MQ=255)
                                 aCTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGtt                                                                                                                     <  4:3080828/149‑1 (MQ=255)
                                  cTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTc                                                                                                                    >  2:2409763/1‑149 (MQ=255)
                                  cTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTc                                                                                                                    >  1:981609/1‑149 (MQ=255)
                                   tGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCa                                                                                                                   >  2:3805255/1‑149 (MQ=255)
                                    gCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAg                                                                                                                  <  1:1479482/149‑1 (MQ=255)
                                        aaCCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATc                                                                                                                                         >  3:1977939/1‑122 (MQ=255)
                                        aaCCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATc                                                                                                                                         <  4:1977939/122‑1 (MQ=255)
                                        aaCCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGAt                                                                                                              >  4:4690418/1‑149 (MQ=255)
                                            cAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCg                                                                                                          <  3:5612097/149‑1 (MQ=255)
                                                  cccGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTgcggc                                                                                                    <  1:1879066/149‑1 (MQ=255)
                                                  cccGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTgcggc                                                                                                    >  1:1007062/1‑149 (MQ=255)
                                                      aTGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTAc                                                                                                                               <  3:1699598/118‑1 (MQ=255)
                                                          ccTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTgc                                                                                                       <  3:6337449/138‑1 (MQ=255)
                                                          ccTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGc                                                                                            >  3:600621/1‑149 (MQ=255)
                                                            ttACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGttt                                                                                                                                  >  4:3868757/1‑109 (MQ=255)
                                                              aCCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGt                                                                                        <  3:727113/149‑1 (MQ=255)
                                                                 gaagaaGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCgg                                                                                     >  1:3872115/1‑149 (MQ=255)
                                                                   agaagGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGaa                                                                                   >  2:6914052/1‑149 (MQ=255)
                                                                          cAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcc                                                                            <  1:6644220/149‑1 (MQ=255)
                                                                          cAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcc                                                                            <  2:3548308/149‑1 (MQ=255)
                                                                              gcCATAAAGCCTATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccacc                                                                        >  3:4163470/1‑149 (MQ=255)
                                                                              gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccacc                                                                        <  4:5758762/149‑1 (MQ=255)
                                                                                 aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAg                                                                     >  3:2505670/1‑149 (MQ=255)
                                                                                  tAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc                                                                    >  2:3261792/1‑149 (MQ=255)
                                                                                   aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGca                                                                   <  1:302515/149‑1 (MQ=255)
                                                                                   aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGca                                                                   >  2:937732/1‑149 (MQ=255)
                                                                                        cGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATccc                                                                                                           >  3:48788/1‑104 (MQ=255)
                                                                                         gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgc                                                             >  4:1773207/1‑149 (MQ=255)
                                                                                         gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgc                                                             >  2:662892/1‑149 (MQ=255)
                                                                                         gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgc                                                             <  4:4163470/149‑1 (MQ=255)
                                                                                              tcatcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATa                                                        >  3:2954157/1‑149 (MQ=255)
                                                                                                atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaaa                                                      >  2:5177558/1‑149 (MQ=255)
                                                                                                atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaaa                                                      >  1:2587867/1‑149 (MQ=255)
                                                                                                atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGAGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaaa                                                      >  2:6419616/1‑149 (MQ=255)
                                                                                                 tcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaaaa                                                     >  3:6924604/1‑149 (MQ=255)
                                                                                                 tcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaaaa                                                     >  3:4635102/1‑149 (MQ=255)
                                                                                                  ccACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaac                                                    >  4:4350759/1‑149 (MQ=255)
                                                                                                     cTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                 >  1:4600947/1‑149 (MQ=255)
                                                                                                      tCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAaataa                                                                               <  1:5302201/118‑1 (MQ=255)
                                                                                                      tCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgca                                                >  2:1823303/1‑149 (MQ=255)
                                                                                                        aaCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaa                                              <  4:1641145/149‑1 (MQ=255)
                                                                                                         aCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaac                                             <  1:6205842/149‑1 (MQ=255)
                                                                                                           gTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCg                                                                                      >  2:5384055/1‑106 (MQ=255)
                                                                                                           gTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATc                                                                             <  2:6033578/115‑1 (MQ=255)
                                                                                                              gACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAAccgccg                                        <  3:3968809/149‑1 (MQ=255)
                                                                                                                 gTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaa                                     >  2:1261460/1‑149 (MQ=255)
                                                                                                                 gTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaa                                     >  2:2228547/1‑149 (MQ=255)
                                                                                                                    gTTGCAGAACACCATCACCATCGACAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgg                                  >  1:5682758/1‑149 (MQ=255)
                                                                                                                      tGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacggaa                                <  4:6732966/149‑1 (MQ=255)
                                                                                                                      tGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacggaa                                <  1:6453139/149‑1 (MQ=255)
                                                                                                                       gCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAg                               >  3:5312508/1‑149 (MQ=255)
                                                                                                                       gCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAg                               <  3:1361472/149‑1 (MQ=255)
                                                                                                                          gAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCca                            >  2:7197981/1‑149 (MQ=255)
                                                                                                                          gAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCca                            >  3:181873/1‑149 (MQ=255)
                                                                                                                            acacCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcaca                          >  1:3311453/1‑149 (MQ=255)
                                                                                                                                catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGc                      >  2:5866921/1‑149 (MQ=255)
                                                                                                                                catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGc                      >  4:2738539/1‑149 (MQ=255)
                                                                                                                                catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGc                      >  1:4966719/1‑149 (MQ=255)
                                                                                                                                   caccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGTCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCg                   >  1:5930744/1‑149 (MQ=255)
                                                                                                                                   caccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCg                   >  4:6066438/1‑149 (MQ=255)
                                                                                                                                    accatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGc                  <  1:280119/149‑1 (MQ=255)
                                                                                                                                      catcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCGGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCtt                >  1:1269444/1‑149 (MQ=255)
                                                                                                                                      catcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCtt                >  4:4985093/1‑149 (MQ=255)
                                                                                                                                        tcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc              >  3:3687358/1‑149 (MQ=255)
                                                                                                                                        tcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc              <  2:2587867/149‑1 (MQ=255)
                                                                                                                                        tcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc              <  2:289229/149‑1 (MQ=255)
                                                                                                                                        tcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc              >  4:291173/1‑149 (MQ=255)
                                                                                                                                        tcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc              >  2:5676951/1‑149 (MQ=255)
                                                                                                                                        tcGACAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATGATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc              >  1:6915501/1‑149 (MQ=255)
                                                                                                                                         cGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCa             <  4:6550878/149‑1 (MQ=255)
                                                                                                                                          gACAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAg            >  2:4208204/1‑149 (MQ=255)
                                                                                                                                           aCAGCAACATGCCCATTACGATGGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa           >  1:34397/1‑149 (MQ=255)
                                                                                                                                           aCAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                       >  4:3592347/1‑89 (MQ=255)
                                                                                                                                           aCAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAg                       <  2:4851559/137‑1 (MQ=255)
                                                                                                                                           aCAGCAACACGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAg                       >  1:4851559/1‑137 (MQ=255)
                                                                                                                                             aGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAAt         >  3:6629247/1‑149 (MQ=255)
                                                                                                                                              gCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc        >  3:1166792/1‑149 (MQ=255)
                                                                                                                                              gCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc        <  4:1106286/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATC  >  NC_000913/2920953‑2921249

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