Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,157,277 C→A 5.3% R219S (CGC→AGC)  dkgA → 2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,157,2770CA 5.3% 243.9 / 2.7 75R219S (CGC→AGC) dkgA2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base A (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTC  >  NC_000913/3157132‑3157424
                                                                                                                                                 |                                                                                                                                                   
ccAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCt                                                                                                                                                  <  1:17303/149‑1 (MQ=255)
ccAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCt                                                                                                                                                  >  5:257777/1‑149 (MQ=255)
 cAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTg                                                                                                                                                 >  4:71906/1‑149 (MQ=255)
       cacGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATc                                                                                                                                            <  3:104272/148‑1 (MQ=255)
        acGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTg                                                                                                                                          <  6:257777/149‑1 (MQ=255)
        acGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTg                                                                                                                                          <  6:241491/149‑1 (MQ=255)
            cTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACACCGGCGCAGATTGTTATCAGCTGGCATCTGGa                                                                                                                                        >  6:157244/1‑147 (MQ=255)
             tGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGAATGTTATCCGCTGGCATCTGGATAg                                                                                                                                     >  7:235057/1‑149 (MQ=255)
             tGGAAAGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGa                                                                                                                                        >  6:237659/1‑146 (MQ=255)
              ggAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGAAGGGAAAGGAGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCAGGCGCAGATTGTTATCCGCTGCCATCTGTATAGc                                                                                                                                    >  2:220747/1‑149 (MQ=255)
                 actcgACACACAAAATGCAGACCGAATCCTGGAGCCCATTAGGGCAAGGAGGTAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATTCGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGc                                                                                                                                 <  3:280296/146‑1 (MQ=255)
                     gACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCtgg                                                                                                                             >  8:33629/1‑149 (MQ=255)
                       cacacaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCtggtg                                                                                                                           <  5:237659/149‑1 (MQ=255)
                        acacaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCtgg                                                                                                                             >  1:383585/1‑146 (MQ=255)
                        acacaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCtgg                                                                                                                             <  2:383585/146‑1 (MQ=255)
                        acacaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTAACCGCTGGCTACTGCATAGCTGACtggtgg                                                                                                                          >  8:178643/1‑149 (MQ=255)
                         cacaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTg                                                                                                                                                 >  8:186338/1‑125 (MQ=255)
                         cacaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTg                                                                                                                                                 <  7:186338/125‑1 (MQ=255)
                          acaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGAt                                                                                                                                       >  5:15094/1‑134 (MQ=255)
                          acaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGAt                                                                                                                                       <  6:15094/134‑1 (MQ=255)
                          acaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTggtggt                                                                                                                         >  7:303639/1‑148 (MQ=255)
                          acaAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCAGCTGGCATCTGGATAGCGGCCTggtggt                                                                                                                         >  4:329657/1‑148 (MQ=255)
                             aaaTCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATc                                                                                                                     <  3:71906/149‑1 (MQ=255)
                             aaaTCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATc                                                                                                                     <  5:290051/149‑1 (MQ=255)
                              aaTCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATcc                                                                                                                    <  3:362966/149‑1 (MQ=255)
                                  cAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCt                                                                                                                                           <  8:106618/122‑1 (MQ=255)
                                  cAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCt                                                                                                                                           >  7:106618/1‑122 (MQ=255)
                                  cAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGaa                                                                                                                <  5:341003/149‑1 (MQ=255)
                                       cGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCgg                                                                                                           >  2:242191/1‑149 (MQ=255)
                                             cTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGCTCCCGAAATCGGTCacac                                                                                                     >  3:360384/1‑149 (MQ=255)
                                              tGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCg                                                                                                            <  6:49409/141‑1 (MQ=255)
                                              tGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCg                                                                                                            >  5:49409/1‑141 (MQ=255)
                                                    ccATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACg                                                                                              >  1:286912/1‑149 (MQ=255)
                                                                   gggAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGTTGAACCCGAAAACGTACACACGTACACGATTTGCCGAAAAGtt                                                                               >  5:197648/1‑149 (MQ=255)
                                                                    ggAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACttt                                                                              >  6:251220/1‑149 (MQ=255)
                                                                      aaaGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGa                                                                            >  2:14369/1‑149 (MQ=255)
                                                                      aaaGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGa                                                                            >  7:32152/1‑149 (MQ=255)
                                                                        aGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATg                                                                          >  4:290923/1‑149 (MQ=255)
                                                                          gCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCAGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTc                                                                        <  3:329657/149‑1 (MQ=255)
                                                                            gTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGAt                                                                           <  2:66219/144‑1 (MQ=255)
                                                                            gTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGAt                                                                           >  1:66219/1‑144 (MQ=255)
                                                                            gTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTg                                                                      <  1:132378/149‑1 (MQ=255)
                                                                            gTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTg                                                                      >  2:132378/1‑149 (MQ=255)
                                                                             ttttCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTgg                                                                     <  1:242191/149‑1 (MQ=255)
                                                                                tCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGa                                                                   <  8:303639/148‑1 (MQ=255)
                                                                                tCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGAt                                                                  >  7:40339/1‑149 (MQ=255)
                                                                                  gATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAACTCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGAttt                                                                >  5:244926/1‑149 (MQ=255)
                                                                                   aTCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGcc                                                                                       <  1:315651/125‑1 (MQ=255)
                                                                                   aTCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGcc                                                                                       >  2:315651/1‑125 (MQ=255)
                                                                                    tCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACttt                                                                              <  2:36888/133‑1 (MQ=255)
                                                                                    tCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACttt                                                                              >  1:36888/1‑133 (MQ=255)
                                                                                      aGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGTAAACTTTGATGTCTGGGATTTCCGt                                                            >  5:313310/1‑149 (MQ=255)
                                                                                      aGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGt                                                            <  3:290923/149‑1 (MQ=255)
                                                                                        aaaGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGtct                                                          <  6:313310/149‑1 (MQ=255)
                                                                                         taGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGtct                                                          <  7:178643/147‑1 (MQ=255)
                                                                                          aGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGc                                                                                                                                    >  5:273840/1‑73 (MQ=255)
                                                                                          aGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGc                                                                                                                                    <  6:273840/73‑1 (MQ=255)
                                                                                            tCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGa                                                       >  7:136683/1‑148 (MQ=255)
                                                                                                   cgATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACg                                               <  2:286912/149‑1 (MQ=255)
                                                                                                    gATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACg                                               >  8:346741/1‑148 (MQ=255)
                                                                                                        tGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAGATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTc                                          >  8:290135/1‑149 (MQ=255)
                                                                                                          gCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCgg                                        >  5:184516/1‑149 (MQ=255)
                                                                                                           cAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGc                                       <  4:360384/149‑1 (MQ=255)
                                                                                                                      ggCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCaaa                            <  7:33629/149‑1 (MQ=255)
                                                                                                                         aaaaCCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACt                         <  5:251220/149‑1 (MQ=255)
                                                                                                                         aaaaCCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGACGAAAACTTTGATGTCTGGGATTTCCGTCTCGAAAAAGACGAACTCGGCGAAATTGCAAAACt                         >  7:225321/1‑149 (MQ=255)
                                                                                                                          aaaCCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTc                        <  1:14369/149‑1 (MQ=255)
                                                                                                                           aaCCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTc                        >  3:165230/1‑148 (MQ=255)
                                                                                                                               ccGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCg                       <  8:215061/145‑1 (MQ=255)
                                                                                                                               ccGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCg                       >  7:215061/1‑145 (MQ=255)
                                                                                                                               ccGGCGCAGATTGTTATCAGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCa                   <  5:157244/149‑1 (MQ=255)
                                                                                                                                 ggCGAAGATTGTTATACGCTGGCATCTGGATAGCTGCCTGTTTGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGAAGAACTCGGCGAAATTGCAAAACTCGATCAgg                 <  1:220747/149‑1 (MQ=255)
                                                                                                                                  gcgcAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAggg                <  8:32152/149‑1 (MQ=255)
                                                                                                                                     cAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAGATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCat             <  7:290135/149‑2 (MQ=255)
                                                                                                                                       gATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGc           <  8:40339/149‑1 (MQ=255)
                                                                                                                                                 cGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTc  <  8:136683/148‑1 (MQ=255)
                                                                                                                                                 cGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTc  <  7:346741/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTAGCGCAAGGAGGGAAAGGCGTTTTCGATCAGAAAGTCATTCGCGATCTGGCAGATAAATACGGCAAAACCCCGGCGCAGATTGTTATCCGCTGGCATCTGGATAGCGGCCTGGTGGTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTC  >  NC_000913/3157132‑3157424

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: ‑