Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,376 G→A 7.8% P14L (CCC→CTC)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,3760GA 7.8% 785.4 / 38.7 255P14L (CCC→CTC) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (122/113);  new base A (10/10);  total (132/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

Only 100 of 263 total aligned reads displayed.

GGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTC  >  NC_000913/53234‑53524
                                                                                                                                              |                                                                                                                                                    
gggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGc                                                                                                                                                >  1:1700599/1‑149 (MQ=255)
gggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGc                                                                                                                                                >  1:415541/1‑149 (MQ=255)
gggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGc                                                                                                                                                >  4:419406/1‑149 (MQ=255)
     aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGGGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGa                                                                                                                                           <  2:1285383/149‑1 (MQ=255)
     aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGa                                                                                                                                           <  3:776221/149‑1 (MQ=255)
      aTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAg                                                                                                                                          >  3:2370164/1‑149 (MQ=255)
        aaGGGCGGAGTGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAAGAGAGTGGCATCTGCACAAATAACCAGGTCGACCGGCCCCTTACGCTGTGGAAACTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                        <  2:321476/149‑1 (MQ=255)
          gggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGAt                                                                                                                                      >  3:1046246/1‑149 (MQ=255)
          gggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGAt                                                                                                                                      >  1:1828181/1‑149 (MQ=255)
                aGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaac                                                                                                                                <  3:205065/149‑1 (MQ=255)
                     gAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGACCACAACACGtt                                                                                                                           <  3:1448967/149‑1 (MQ=255)
                      aGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTg                                                                                                                          <  3:318489/149‑1 (MQ=255)
                          gCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtt                                                                                                                      >  3:2389832/1‑149 (MQ=255)
                          accaaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtt                                                                                                                      >  3:767243/4‑149 (MQ=255)
                           cAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCGCAACACGTTGGGttt                                                                                                                     <  1:24737/149‑1 (MQ=255)
                             aaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGa                                                                                                                                           >  4:1145440/1‑125 (MQ=255)
                             aaCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                                                                                                                   >  2:1013928/1‑149 (MQ=255)
                               ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc                                                                                                                 <  1:908525/149‑1 (MQ=255)
                               ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc                                                                                                                 <  2:41851/149‑1 (MQ=255)
                               ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc                                                                                                                 >  4:2266480/1‑149 (MQ=255)
                                    cATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTa                                                                                                            <  1:156940/149‑1 (MQ=255)
                                     aTCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg                                                                                                           >  4:1101736/1‑149 (MQ=255)
                                     aTCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg                                                                                                           >  4:2213272/1‑149 (MQ=255)
                                       cGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTGACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGtt                                                                                                         <  4:2408930/149‑1 (MQ=255)
                                            cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCa                                                                                                    <  2:264924/149‑1 (MQ=255)
                                            cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGGGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCa                                                                                                    >  4:1061980/1‑149 (MQ=255)
                                             ccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAg                                                                                                   <  3:1702033/149‑1 (MQ=255)
                                             ccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAg                                                                                                   >  1:1639504/1‑149 (MQ=255)
                                               ggTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                                                                 <  1:2223573/149‑1 (MQ=255)
                                                  tGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAttg                                                                                              <  3:1870381/149‑2 (MQ=255)
                                                      aaGGAGAGTGGCAGCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAAc                                                                                          <  1:578517/149‑1 (MQ=255)
                                                           gagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAgg                                                                                     >  3:663612/1‑149 (MQ=255)
                                                              tGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCACTCCCGGCGAGCCCGCCGGGATTg                                                                                                                                        <  1:790906/95‑1 (MQ=38)
                                                              tGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGgcgc                                                                                  >  1:862183/1‑149 (MQ=255)
                                                                gCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTg                                                                                <  1:1032500/149‑1 (MQ=255)
                                                                  aTCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGc                                                                              >  3:1713202/1‑149 (MQ=255)
                                                                  aTCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGc                                                                              >  1:1775004/1‑149 (MQ=255)
                                                                       cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgtt                                                                         <  2:231017/149‑1 (MQ=255)
                                                                           aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTc                                                                     >  3:447284/1‑149 (MQ=255)
                                                                           aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTc                                                                     <  4:2444880/149‑1 (MQ=255)
                                                                           aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTc                                                                     >  2:769173/1‑149 (MQ=255)
                                                                            acaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                    <  2:2085812/149‑1 (MQ=255)
                                                                             caacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCt                                                                   <  1:769173/149‑1 (MQ=255)
                                                                             caacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCt                                                                   <  3:1101736/149‑1 (MQ=255)
                                                                              aacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGc                                                                                                                                                <  1:1120405/71‑1 (MQ=255)
                                                                              aacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTg                                                                  <  2:862183/149‑1 (MQ=255)
                                                                              aacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTg                                                                  >  4:1852361/1‑149 (MQ=255)
                                                                               acCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGc                                                                 <  2:1514010/149‑1 (MQ=255)
                                                                                 cAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCAt                                                               <  2:2421197/149‑1 (MQ=255)
                                                                                  aGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCCGGCACGTTGTTCCTGCATc                                                              <  4:1713202/149‑1 (MQ=255)
                                                                                       gACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTTGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCt                                                         >  2:1058967/1‑149 (MQ=255)
                                                                                       gACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                    <  2:1896756/138‑1 (MQ=255)
                                                                                        aCCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATc                                                                                                                                     >  1:1007002/1‑72 (MQ=255)
                                                                                         ccGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtg                                                       <  3:1173995/149‑1 (MQ=255)
                                                                                           ggCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTATAACGTAGGCGCTGGCACGTTGTTCCTGCACGTAGCTtgct                                                     <  2:200912/149‑1 (MQ=255)
                                                                                            gCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCGCAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgctg                                                    <  3:251798/149‑1 (MQ=255)
                                                                                            gCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgctg                                                    >  4:1255167/1‑149 (MQ=255)
                                                                                             ccACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCAGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgctgc                                                   >  3:1821626/1‑149 (MQ=255)
                                                                                               aCTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGctt                                                 >  2:956502/1‑149 (MQ=255)
                                                                                                cTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGcttc                                                <  4:566849/149‑1 (MQ=255)
                                                                                                 tCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCGCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGcttct                                               >  4:1607222/1‑149 (MQ=255)
                                                                                                 tCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                                                                                                                  <  4:750513/82‑1 (MQ=255)
                                                                                                 tCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGcttct                                               <  1:1400103/149‑1 (MQ=255)
                                                                                                 tCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGcttct                                               >  1:2033211/1‑149 (MQ=255)
                                                                                                 tCACGCTGCGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGcttct                                               >  1:362929/1‑149 (MQ=255)
                                                                                                   aCGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc                                                                                                                >  1:470857/1‑82 (MQ=255)
                                                                                                   aCGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc                                                                                                                <  2:470857/82‑1 (MQ=255)
                                                                                                    cGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTgg                                                                                                                         >  1:2294804/1‑72 (MQ=255)
                                                                                                     gCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCg                                           >  4:527542/1‑149 (MQ=255)
                                                                                                     gCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTGTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCg                                           >  4:212785/1‑149 (MQ=255)
                                                                                                             aGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAg                                                                                      >  2:113206/1‑98 (MQ=255)
                                                                                                              gCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTAc                                  >  3:1761827/1‑149 (MQ=255)
                                                                                                                tGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACgg                                <  1:114766/149‑1 (MQ=255)
                                                                                                                  gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGtt                              >  2:1805659/1‑149 (MQ=255)
                                                                                                                  gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATATTAACGTAGGCGCTGGCa                                                                             >  3:1636780/1‑102 (MQ=255)
                                                                                                                     aaCTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcat                           >  1:1670650/1‑149 (MQ=255)
                                                                                                                     aaCTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcat                           >  3:1383787/1‑149 (MQ=255)
                                                                                                                       cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatca                         <  4:889285/149‑1 (MQ=255)
                                                                                                                       cTAAGTCCAGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCtt                                                        <  3:2175582/118‑1 (MQ=255)
                                                                                                                         aaGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGc                       >  2:178852/1‑149 (MQ=255)
                                                                                                                           gTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCAt                     >  4:295077/1‑149 (MQ=255)
                                                                                                                            tCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatca                         >  1:1369993/1‑144 (MQ=255)
                                                                                                                              cGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCg                  >  1:1891251/1‑149 (MQ=255)
                                                                                                                               gggCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCgg                 >  4:801858/1‑149 (MQ=255)
                                                                                                                               gggCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCgg                 <  2:1700599/149‑1 (MQ=255)
                                                                                                                                 gCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTa               >  1:96770/1‑149 (MQ=255)
                                                                                                                                 atcactCCCGGCGAGCCCGCCGGGATTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtc                             >  4:333356/6‑135 (MQ=255)
                                                                                                                                 atcactCCCGGCGAGCCCGCCGGGATTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtc                             <  3:333356/130‑1 (MQ=255)
                                                                                                                                  tcAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTCGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTAt              >  2:670554/2‑149 (MQ=255)
                                                                                                                                    aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGc            <  2:1871186/149‑1 (MQ=255)
                                                                                                                                     aTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCt                                                                   >  4:1938400/1‑93 (MQ=255)
                                                                                                                                     aTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCt                                                                   <  3:1938400/93‑1 (MQ=255)
                                                                                                                                      tCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCAc          >  3:910431/1‑149 (MQ=255)
                                                                                                                                      tCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCAc          >  3:53409/1‑149 (MQ=255)
                                                                                                                                      tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCa           <  2:1233154/148‑1 (MQ=255)
                                                                                                                                      tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCAc          >  4:627772/1‑149 (MQ=255)
                                                                                                                                        ccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGa        <  2:1639504/149‑1 (MQ=255)
                                                                                                                                        ccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGa        <  1:242900/149‑1 (MQ=255)
                                                                                                                                         cggcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGAt       <  1:1185723/149‑1 (MQ=255)
                                                                                                                                          ggcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATc      <  1:172061/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTC  >  NC_000913/53234‑53524

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