Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 299,652 T→A 6.1% E428V (GAG→GTG)  paoC ← PaoABC aldehyde oxidoreductase, Moco‑containing subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913299,6520TA 6.1% 207.0 / 2.0 64E428V (GAG→GTG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base A (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAGTGTCATTCAGGATG  >  NC_000913/299508‑299795
                                                                                                                                                |                                                                                                                                               
aGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTcccc                                                                                                                                             <  8:142421/149‑1 (MQ=255)
  aaTGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGt                                                                                                                                           >  1:32677/1‑149 (MQ=255)
   aTGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTc                                                                                                                                          <  8:60480/149‑1 (MQ=255)
      gtgtAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTcgcg                                                                                                                                       >  1:190610/1‑149 (MQ=255)
                         tGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCg                                                                                                                    >  1:156487/1‑149 (MQ=255)
                         tGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCg                                                                                                                    >  4:312564/1‑149 (MQ=255)
                               tCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTgcgc                                                                                                              <  8:85204/149‑1 (MQ=255)
                               tCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTgcgc                                                                                                              <  2:156487/149‑1 (MQ=255)
                                 aTGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCgctg                                                                                                            >  7:216392/1‑149 (MQ=255)
                                    tCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTc                                                                                                                               >  2:162454/1‑127 (MQ=255)
                                    tCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTc                                                                                                                               <  1:162454/127‑1 (MQ=255)
                                    tCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTCCtt                                                                                                         >  1:45420/1‑149 (MQ=255)
                                       gTTTATACGGTAACGGTTCCGTTTTGTTCGAGGTGAACCCTAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTATCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCa                                                                                                      <  3:79370/149‑1 (MQ=255)
                                          tCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCg                                                                                                                                            >  1:277504/1‑108 (MQ=255)
                                          tCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCg                                                                                                                                            <  2:277504/108‑1 (MQ=255)
                                             aCGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTGGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCaaa                                                                                                >  4:273088/1‑149 (MQ=255)
                                                 tAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTa                                                                                            >  1:212597/1‑149 (MQ=255)
                                                   aCGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTAt                                                                                           <  2:32677/148‑1 (MQ=255)
                                                    cGGTGCCGGTTTCTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTCGCCGACTACCCACTCCCCGTCGCGCACCTGGCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATcc                                                                                         <  4:9462/149‑1 (MQ=255)
                                                     ggTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCg                                                                                        <  1:241397/149‑1 (MQ=255)
                                                           gTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCATATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCgg                                                                                  >  7:184396/1‑149 (MQ=255)
                                                            ttttGTTCGAGGTGAAACCGATCACCCGATCTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAAAACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCGGCTGCCAGCCAAATTTATCCGCTCCGGt                                                                                 <  8:258440/149‑1 (MQ=255)
                                                                   cGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGACGACTAGCCACTCCCCTTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAACCAAATTTATCCGCTCCGGTGCGCaa                                                                           <  2:80412/148‑1 (MQ=255)
                                                                    gAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAg                                                                          >  2:229885/1‑148 (MQ=255)
                                                                      ggTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCAc                                                                       >  1:365429/1‑149 (MQ=255)
                                                                       gTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACt                                                                      <  2:212597/149‑1 (MQ=255)
                                                                          aaCCCGAGCACCCGATTTTTCCAGCATATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGa                                                                   <  8:184396/149‑1 (MQ=255)
                                                                          aaCCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGa                                                                   <  2:365429/149‑1 (MQ=255)
                                                                             ccGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTCGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATaa                                                                <  3:273088/149‑1 (MQ=255)
                                                                             ccGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGAt                                                                  >  8:133574/1‑147 (MQ=255)
                                                                             ccGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGAt                                                                  >  4:195259/1‑147 (MQ=255)
                                                                             ccGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATaa                                                                <  7:202837/149‑1 (MQ=255)
                                                                                     ccGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATc                                                                                          <  5:153599/115‑1 (MQ=255)
                                                                                     ccGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATc                                                                                          >  6:153599/1‑115 (MQ=255)
                                                                                      cGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGc                                                       <  6:208843/149‑1 (MQ=255)
                                                                                       gATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCg                                                      <  8:216392/149‑1 (MQ=255)
                                                                                        aTTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCg                                                                                   <  5:53608/119‑1 (MQ=255)
                                                                                        aTTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCg                                                                                   >  6:53608/1‑119 (MQ=255)
                                                                                         tttttCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGc                                                       >  7:166580/1‑146 (MQ=255)
                                                                                               cagcagATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCAACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCAAGACAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAag                                               >  4:187993/1‑148 (MQ=255)
                                                                                                agcagATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAag                                               >  3:1623/1‑147 (MQ=255)
                                                                                                    gATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAgcgc                                         <  2:190610/149‑1 (MQ=255)
                                                                                                      ttattGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCgt                                       >  5:190073/1‑149 (MQ=255)
                                                                                                         ttGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCgtcgg                                    >  3:329223/1‑149 (MQ=255)
                                                                                                         ttGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCgtcgg                                    <  7:133574/149‑1 (MQ=255)
                                                                                                         ttGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCgtcgg                                    <  3:195259/149‑1 (MQ=255)
                                                                                                           gCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCgtcggg                                   <  6:190073/148‑1 (MQ=255)
                                                                                                             gAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCgtcgggtcg                                >  7:319279/1‑149 (MQ=255)
                                                                                                                 gCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGCCGGCGAGAGAAGCAGCGCGTCGGGGgggcgg                            >  6:225878/1‑149 (MQ=255)
                                                                                                                 gCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTcggcgg                            >  8:368235/1‑149 (MQ=255)
                                                                                                                 gCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTCCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTcggcgg                            <  2:45420/149‑1 (MQ=255)
                                                                                                                         cAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAAcc                    >  7:156024/1‑149 (MQ=255)
                                                                                                                             acCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCTCGTCGGGTCGGCGGGGTCAACCTGAg                >  8:163457/1‑149 (MQ=255)
                                                                                                                               cGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAgtg              >  2:183905/1‑149 (MQ=255)
                                                                                                                               cGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAgtg              <  3:312564/149‑1 (MQ=255)
                                                                                                                                 tGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAgtgt             >  6:181641/1‑148 (MQ=255)
                                                                                                                                    ccGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAGTGTCAtt         >  8:228422/1‑149 (MQ=255)
                                                                                                                                    ccGACTAGCCACACCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGa                 >  2:167256/1‑141 (MQ=255)
                                                                                                                                    ccGACTAGCCACACCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGa                 >  2:162149/1‑141 (MQ=255)
                                                                                                                                    ccGACTAGCCACACCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGa                 <  1:167256/141‑1 (MQ=255)
                                                                                                                                    ccGACTAGCCACACCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGa                 <  1:162149/141‑1 (MQ=255)
                                                                                                                                       aCTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAGTGTCATTCAg      >  8:41470/1‑149 (MQ=255)
                                                                                                                                         tAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAGTGTCATTCAGGa    >  2:115170/1‑149 (MQ=255)
                                                                                                                                           gCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAGTGTCATTCAGGATg  >  4:88157/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAGTGTCATTCAGGATG  >  NC_000913/299508‑299795

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