Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,566 T→A 5.8% K333N (AAA→AAT fbaB ← fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,5660TA 5.8% 301.1 / 6.4 98K333N (AAA→AATfbaBfructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base T (44/48);  new base A (3/3);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

TAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAA  >  NC_000913/2177418‑2177709
                                                                                                                                                    |                                                                                                                                               
tAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGt                                                                                                                                                 <  6:195989/149‑1 (MQ=255)
   gACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGttt                                                                                                                                               >  2:113994/1‑148 (MQ=255)
   gACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTc                                                                                                                                              <  2:302188/149‑1 (MQ=255)
   gACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTAAATCAGTTTc                                                                                                                                              <  4:392817/149‑1 (MQ=255)
    aCGCTTCAGCGTCGCATCATGGATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGTCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCa                                                                                                                                             <  8:190335/149‑1 (MQ=255)
    aCGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCa                                                                                                                                             >  5:263042/1‑149 (MQ=255)
    aCGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCa                                                                                                                                             <  4:43090/149‑1 (MQ=255)
        ttCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGcc                                                                                                                                         >  1:94689/1‑149 (MQ=255)
            gCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTGAAACTGATTAACGCCGtgc                                                                                                                                     <  3:118128/149‑3 (MQ=38)
               tCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGc                                                                                                                                   >  3:319329/1‑148 (MQ=255)
                   aTCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATc                                                                                                                               >  2:121369/1‑148 (MQ=255)
                    tCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGa                                                                                                                             <  5:260817/149‑1 (MQ=255)
                      aGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGAt                                                                                                                            <  7:255778/148‑1 (MQ=255)
                       ggCATAGGTTGCCGGATGCAGCGTAAACCCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTAAATCAGTTTCACGCCGTCAGCCATCGAttt                                                                                                                          >  6:316344/1‑149 (MQ=255)
                        gCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTc                                                                                                                         >  7:403020/1‑149 (MQ=255)
                        gCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTc                                                                                                                         <  7:75452/149‑1 (MQ=255)
                         cATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCt                                                                                                                        <  1:398075/149‑1 (MQ=255)
                           tAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCtt                                                                                                                       >  8:322326/1‑148 (MQ=255)
                              gTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGaa                                                                                                                    <  7:254178/148‑1 (MQ=255)
                                  ccGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTGAAACTGATTAACGCCGtgca                                                                                                                                    >  8:46886/1‑125 (MQ=39)
                                  ccGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTGAAACTGATTAACGCCGtgca                                                                                                                                    <  7:46886/128‑4 (MQ=39)
                                   cGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCttt                                                                                                              <  2:94689/149‑1 (MQ=255)
                                         cAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTACACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTcc                                                                                                        <  8:212724/149‑1 (MQ=255)
                                           gCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTcc                                                                                                        >  1:51141/1‑147 (MQ=255)
                                           gCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCaa                                                                                                      <  3:262555/149‑1 (MQ=255)
                                            tgTATACGCCTTATCCTTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCaa                                                                                                      <  1:119806/147‑1 (MQ=255)
                                              tAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGAGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGaa                                                                                                   >  6:180010/1‑149 (MQ=255)
                                                aaCGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATATTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATc                                                                                                 <  8:201576/149‑1 (MQ=255)
                                                  cGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTa                                                                                                             <  5:336803/135‑1 (MQ=255)
                                                  cGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTa                                                                                                             >  6:336803/1‑135 (MQ=255)
                                                   gCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGc                                                                                              <  7:428798/149‑1 (MQ=255)
                                                     cTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGccc                                                                                            <  1:121369/149‑1 (MQ=255)
                                                           cGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTcc                                                                                      >  7:118141/1‑149 (MQ=255)
                                                           cGTCCTACGGGGTGGTGGTGTTCGTATAAGGCGTATAAGGAGATAGTAATTTTGCTATCGAGATAAACTTCATGCACGGCGTAAATCAGATTCACGCCGTAATCCATCGATTTCTAGAACGCTTTACGTCCAACAAACAGCCGCATTcc                                                                                      >  1:266436/1‑149 (MQ=255)
                                                            gTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTc                                                                                                         <  6:330399/129‑1 (MQ=255)
                                                            gTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTc                                                                                                         >  5:330399/1‑129 (MQ=255)
                                                            gTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCg                                                                                     >  1:180429/1‑149 (MQ=255)
                                                                  cGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAAAATCAGCCCCATTCCGCCTgc                                                                                <  8:403020/148‑1 (MQ=255)
                                                                      gtggtgTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGtt                                                                           <  5:210897/149‑1 (MQ=255)
                                                                        ggtgTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTttg                                                                         <  6:263042/149‑1 (MQ=255)
                                                                         gtgtTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTttgt                                                                        <  2:51141/149‑1 (MQ=255)
                                                                          tgttgtTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACTTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGATCGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTttgtt                                                                       <  1:113994/149‑1 (MQ=255)
                                                                          tgttgtTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTttgtt                                                                       <  8:153675/149‑1 (MQ=255)
                                                                          tgttgtTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTttgtt                                                                       <  1:72967/149‑1 (MQ=255)
                                                                          tgttgtTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTttgtt                                                                       <  8:256227/149‑1 (MQ=255)
                                                                          tgttgtTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTACAAGAATCAGACCCATTCCGCCTTCGCGgttgtt                                                                       >  4:104914/1‑149 (MQ=255)
                                                                           gttgttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGcc                                                                                                                                         >  1:433092/1‑82 (MQ=255)
                                                                           gttgttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGcc                                                                                                                                         <  2:433092/82‑1 (MQ=255)
                                                                           gttgttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGc                                                                                    <  1:257247/135‑1 (MQ=255)
                                                                           gttgttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGc                                                                                    >  2:257247/1‑135 (MQ=255)
                                                                            ttgttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGa                                                                     >  3:285563/1‑149 (MQ=255)
                                                                              gttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGt                                                                            <  6:372916/140‑1 (MQ=255)
                                                                              gttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGt                                                                            >  5:372916/1‑140 (MQ=255)
                                                                               ttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATaa                                                                  <  4:319329/149‑1 (MQ=255)
                                                                               ttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATaa                                                                  >  8:348910/1‑149 (MQ=255)
                                                                               ttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATaa                                                                  <  6:421142/149‑1 (MQ=255)
                                                                               ttCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTATAGAACGCTTTACGTCCACGAATCAGCCACATTCCGCCTGCGCGTTTGTTGATaa                                                                  >  8:178132/1‑149 (MQ=255)
                                                                                 cGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAGCGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAAcc                                                                >  7:35773/1‑149 (MQ=255)
                                                                                  gTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAAcc                                                                <  1:115391/148‑1 (MQ=255)
                                                                                  gTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCg                                                               >  4:197166/1‑149 (MQ=255)
                                                                                  gTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGCCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAAcc                                                                <  1:49498/148‑1 (MQ=255)
                                                                                    ataAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCa                                                                                                                                             >  7:273719/1‑69 (MQ=255)
                                                                                    ataAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCa                                                                                                                                             <  8:273719/69‑1 (MQ=255)
                                                                                     taAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAg                                                            >  4:368606/1‑149 (MQ=255)
                                                                                       aGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTa                                                          >  5:160057/1‑149 (MQ=255)
                                                                                         gagTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCAGGGCGTTAATCAGTTTCACGCCGTGAGCCATCGATTTCTTGAACGCTTTACGTCCAATAATCAGGCCCATTCCGCCTGCGCGTTTGTTGATAACCGGAGTAc                                                         <  4:146003/146‑1 (MQ=255)
                                                                                         gCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACg                                                        >  2:157534/1‑149 (MQ=255)
                                                                                           gTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCa                                                      <  3:197166/149‑1 (MQ=255)
                                                                                                aGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGc                                                  <  7:348910/148‑1 (MQ=255)
                                                                                                aGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTAAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCa                                                 <  6:366989/149‑1 (MQ=255)
                                                                                                    gATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGc                                             >  3:35229/1‑149 (MQ=255)
                                                                                                    gATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGc                                             >  8:105664/1‑149 (MQ=255)
                                                                                                       aGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGaa                                                                                                                    >  6:18089/1‑75 (MQ=255)
                                                                                                       aGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGaa                                                                                                                    <  5:18089/75‑1 (MQ=255)
                                                                                                        gTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAg                                         <  7:256676/149‑1 (MQ=255)
                                                                                                         tAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCCGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGCTAACCGCAGTACGCACTGCATCGCTGAgg                                        >  3:54785/1‑149 (MQ=255)
                                                                                                         tAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAgg                                        <  8:35773/149‑1 (MQ=255)
                                                                                                         tAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAgg                                        >  1:148700/1‑149 (MQ=255)
                                                                                                           aTTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTc                                      <  1:157534/149‑1 (MQ=255)
                                                                                                                    tCGAGATAAACGTCCTTCACGTCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTTAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCTCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCa                               <  7:178132/147‑1 (MQ=255)
                                                                                                                    tCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCtt                                                                                                                       <  4:394039/59‑1 (MQ=255)
                                                                                                                    tCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCtt                                                                                                                       >  3:394039/1‑59 (MQ=255)
                                                                                                                    tCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCAcc                             <  2:180429/149‑1 (MQ=255)
                                                                                                                    tCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCAcc                             >  8:66333/1‑149 (MQ=255)
                                                                                                                    tCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACCTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCAcc                             >  6:172102/1‑149 (MQ=255)
                                                                                                                      gagaTAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCAccgc                           >  4:371819/1‑149 (MQ=255)
                                                                                                                       agaTAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCAccgcc                          >  3:196453/1‑149 (MQ=255)
                                                                                                                           aaaCGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGACAACGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTCTTGATAACCGCAGTACGCACTGCATCTCTGAGTTCAGTTTCAccgcccg                        >  6:280831/1‑147 (MQ=255)
                                                                                                                             aCGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCcgccgc                    >  2:273996/1‑149 (MQ=255)
                                                                                                                             aCGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCgcagc                    <  6:22905/149‑1 (MQ=255)
                                                                                                                             aCGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCgcag                     >  5:22905/1‑148 (MQ=255)
                                                                                                                                  cTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGACCGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAATACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCAccgc               >  5:230276/1‑149 (MQ=255)
                                                                                                                                    tgaCGGCGTGAAACTGATTAACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCAccgccg             >  2:279957/3‑149 (MQ=255)
                                                                                                                                    gCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCAccgccg             >  3:267486/1‑149 (MQ=255)
                                                                                                                                     cACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCAccgccg             >  3:109705/1‑148 (MQ=255)
                                                                                                                                      aCGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCAccgccg             >  4:112488/1‑147 (MQ=255)
                                                                                                                                      aCGGCGTGAAACTGATTAACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGa           <  1:279957/149‑1 (MQ=255)
                                                                                                                                          cGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGtt        >  4:291886/1‑148 (MQ=255)
                                                                                                                                              aaTCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATc    <  4:35229/148‑1 (MQ=255)
                                                                                                                                               aTCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCaa  <  4:285563/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAA  >  NC_000913/2177418‑2177709

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: