Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,410,990 T→A 7.2% intergenic (+106/+23) dbpA → / ← ttcA ATP‑dependent RNA helicase, specific for 23S rRNA/tRNA s(2)C32 thioltransferase, iron sulfur cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,410,9900TA 7.2% 272.9 / 8.2 80intergenic (+106/+23)dbpA/ttcAATP‑dependent RNA helicase, specific for 23S rRNA/tRNA s(2)C32 thioltransferase, iron sulfur cluster protein
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base A (3/3);  total (43/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 = 1.00e+00

GAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCA  >  NC_000913/1410842‑1411137
                                                                                                                                                    |                                                                                                                                                   
gAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCt                                                                                                                                                     >  1:507007/1‑149 (MQ=255)
 aaGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCtt                                                                                                                                                    <  2:540099/149‑1 (MQ=255)
        aaGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATc                                                                                                                                             <  4:16142/149‑1 (MQ=255)
          ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCgg                                                                                                                                           <  3:618464/149‑1 (MQ=255)
          ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCgg                                                                                                                                           <  2:339524/149‑1 (MQ=255)
          ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCgg                                                                                                                                           <  1:79313/149‑1 (MQ=255)
          ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCgg                                                                                                                                           <  2:801560/149‑1 (MQ=255)
           gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                                                          >  2:96374/1‑149 (MQ=255)
           gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                                                          >  2:516614/1‑149 (MQ=255)
           gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                                                          >  4:754810/1‑149 (MQ=255)
           gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                                                          <  4:546352/149‑1 (MQ=255)
                     ccGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTg                                                                                                                                <  1:197820/149‑1 (MQ=255)
                        ggTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTg                                                                                                                                <  4:898933/146‑1 (MQ=255)
                        ggTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTg                                                                                                                                >  3:898933/1‑146 (MQ=255)
                             ggTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTc                                                                                                                        >  4:387024/1‑149 (MQ=255)
                              gTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCa                                                                                                                       <  1:719995/149‑1 (MQ=255)
                              gTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGAATTGACGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCa                                                                                                                       <  4:3638/149‑1 (MQ=255)
                               ttattaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCAc                                                                                                                      <  3:44785/149‑1 (MQ=255)
                                tattaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACt                                                                                                                     <  2:507007/149‑1 (MQ=255)
                                tattaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACt                                                                                                                     <  2:353972/149‑1 (MQ=255)
                                 attaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACtt                                                                                                                    >  3:112713/1‑149 (MQ=255)
                                   taAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCa                                                                                                                  >  3:714451/1‑149 (MQ=255)
                                     aaaTAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTAt                                                                                                                           <  2:358260/138‑1 (MQ=255)
                                     aaaTAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTAt                                                                                                                           >  1:358260/1‑138 (MQ=255)
                                     aaaTAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAAc                                                                                                                >  3:538816/1‑149 (MQ=255)
                                     aaaTAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAAc                                                                                                                >  4:359876/1‑149 (MQ=255)
                                       ataatGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACca                                                                                                              >  1:386610/1‑149 (MQ=255)
                                        taatGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACcac                                                                                                             >  4:856733/1‑149 (MQ=255)
                                           tGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACAtt                                                                                                          >  1:816878/1‑149 (MQ=255)
                                            gAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTc                                                                                                         >  4:802739/1‑149 (MQ=255)
                                              aaTGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAg                                                                                                       <  3:754810/149‑1 (MQ=255)
                                                tGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGcc                                                                                                     >  1:780115/1‑149 (MQ=255)
                                                       ttGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTc                                                                                              <  1:96374/149‑1 (MQ=255)
                                                          ccGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATc                                                                                           <  3:802285/149‑1 (MQ=255)
                                                             ggTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCaa                                                                                        <  2:386610/149‑1 (MQ=255)
                                                                  aaGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCAc                                                                                                                      >  3:237047/1‑114 (MQ=255)
                                                                  aaGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCAc                                                                                                                      <  4:237047/114‑1 (MQ=255)
                                                                  aaGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAt                                                                                   >  3:78283/1‑149 (MQ=255)
                                                                  aaGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAt                                                                                   <  3:387024/149‑1 (MQ=255)
                                                                  aaGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGAAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAt                                                                                   <  2:567330/149‑1 (MQ=255)
                                                                     agTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAtttt                                                                                <  1:751763/149‑1 (MQ=255)
                                                                     agTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAtttt                                                                                >  2:489698/1‑149 (MQ=255)
                                                                       tAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCa                                                                              <  1:369731/149‑1 (MQ=255)
                                                                            aTCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtcttc                                                                         >  4:443605/1‑149 (MQ=255)
                                                                             tCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtcttct                                                                        >  2:263981/1‑149 (MQ=255)
                                                                             tCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtcttct                                                                        >  3:389758/1‑149 (MQ=255)
                                                                                 aTTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTg                                                                    <  4:562013/149‑1 (MQ=255)
                                                                                  ttCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCAc                                                                                                                      >  2:586155/1‑98 (MQ=255)
                                                                                  ttCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCAc                                                                                                                      <  1:586155/98‑1 (MQ=255)
                                                                                     gggATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTg                                                                <  3:480455/149‑1 (MQ=255)
                                                                                       gATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCt                                                                 >  3:549252/1‑146 (MQ=255)
                                                                                       gATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGcc                                                              >  4:514500/1‑149 (MQ=255)
                                                                                       gATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCACATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCt                                                                 <  4:549252/146‑1 (MQ=255)
                                                                                            ccGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCAc                                                         >  4:484107/1‑149 (MQ=255)
                                                                                                    gACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg                                                 >  3:265443/1‑149 (MQ=255)
                                                                                                    gACGCGGCCAGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg                                                 >  1:49252/1‑149 (MQ=255)
                                                                                                        cTGGCCGCGTCTTATCCGGCCTCTATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGt                                             <  4:643875/149‑1 (MQ=255)
                                                                                                          ggCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTgg                                           >  1:695733/1‑149 (MQ=255)
                                                                                                           gCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTggg                                          >  3:507956/1‑149 (MQ=255)
                                                                                                            ccGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGa                                         >  4:879134/1‑149 (MQ=255)
                                                                                                              gcgTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGAtc                                       <  4:112713/149‑1 (MQ=255)
                                                                                                                   ttATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTc                                  >  4:144998/1‑149 (MQ=255)
                                                                                                                      tCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTcgcg                               >  2:1069/1‑149 (MQ=255)
                                                                                                                      tCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGa                              >  1:676779/1‑149 (MQ=255)
                                                                                                                       ccGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGa                              <  2:250435/149‑1 (MQ=255)
                                                                                                                         ggCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATc                            <  2:816878/149‑1 (MQ=255)
                                                                                                                               cATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGc                      >  2:416477/1‑149 (MQ=255)
                                                                                                                                   agagTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGa                  <  3:802739/149‑1 (MQ=255)
                                                                                                                                   agagTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGa                  >  3:680966/1‑149 (MQ=255)
                                                                                                                                    gagTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAg                   >  1:714034/1‑147 (MQ=255)
                                                                                                                                    gagTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAg                   <  2:714034/147‑1 (MQ=255)
                                                                                                                                      gTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATc                            <  2:478201/136‑1 (MQ=255)
                                                                                                                                      gTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATc                            >  1:478201/1‑136 (MQ=255)
                                                                                                                                      gTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCa               <  4:730365/149‑1 (MQ=255)
                                                                                                                                         gCCCGATGCGCAAGCGTATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCAccc            >  4:667779/1‑149 (MQ=255)
                                                                                                                                           ccGATGCGCAAGCGTATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCGCGCCGGCTGTAGTGGGATCTCTTCGCGAGCCAACGCCAGATCAgcccc          >  2:128552/1‑149 (MQ=255)
                                                                                                                                           ccGATGCGCAAGCGTATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCAccccc          >  3:327596/1‑149 (MQ=255)
                                                                                                                                                aCGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTaa     <  1:312756/149‑1 (MQ=255)
                                                                                                                                                 cGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc   <  2:676779/149‑1 (MQ=255)
                                                                                                                                                  gCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc   >  1:925754/1‑149 (MQ=255)
                                                                                                                                                  gCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc   >  4:531392/1‑149 (MQ=255)
                                                                                                                                                   cTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCa  >  3:49769/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCA  >  NC_000913/1410842‑1411137

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