Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,475,772 G→T 5.3% R199R (CGC→CGA)  dusA ← tRNA dihydrouridine synthase DusA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,475,7720GT 5.3% 285.8 / 3.7 88R199R (CGC→CGA) dusAtRNA dihydrouridine synthase DusA
Reads supporting (aligned to +/- strand):  ref base G (44/39);  new base T (2/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GATAAGGGTTGTGATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGCTGTCACGGCCGTTGATAC  >  NC_002947/2475624‑2475920
                                                                                                                                                    |                                                                                                                                                    
gATAAGGGTTGTGATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAg                                                                                                                                                      >  3:126537/1‑149 (MQ=255)
  tAAGGGTTGTGATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAgcg                                                                                                                                                    <  7:250243/149‑1 (MQ=255)
     gggTTGTGATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAg                                                                                                                                                 >  3:25373/1‑149 (MQ=255)
      ggTTGTGATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAgc                                                                                                                                                >  3:157471/1‑149 (MQ=255)
        ttGTGATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAgcgg                                                                                                                                              <  7:111472/149‑1 (MQ=255)
          gtgATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAgcggcg                                                                                                                                            >  4:14016/1‑149 (MQ=255)
           tgATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAgcggcg                                                                                                                                            <  3:345581/148‑1 (MQ=255)
           tgATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGATGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGAGATGTCATAGCGCAgcggcg                                                                                                                                            >  2:309926/1‑148 (MQ=255)
               aCGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAAtc                                                                                                                                       >  8:328190/1‑149 (MQ=255)
                cGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAAtct                                                                                                                                      >  7:194362/1‑149 (MQ=255)
                          cccAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTtctc                                                                                                                            <  3:14016/149‑1 (MQ=255)
                          cccAGCATCACGCCATCGAAGGTTTCAAGGAGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTtctc                                                                                                                            >  1:31787/1‑149 (MQ=255)
                           ccAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTctct                                                                                                                           <  4:48572/149‑1 (MQ=255)
                            caccATGACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCGCTTCAACAGCGCTGCGATGTCATAGCGCAGCGGCGGAATATGGCGGTAATCtt                                                                                                                          >  4:85320/4‑149 (MQ=255)
                              gCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGTATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTc                                                                                                                        <  4:198619/149‑1 (MQ=255)
                              gCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTTTTCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTc                                                                                                                        <  1:285663/149‑1 (MQ=255)
                                    cGCCATCGAAGGTTTCAAGGTGTGCCCGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGa                                                                                                                   >  1:169722/1‑148 (MQ=255)
                                     gCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACa                                                                                                                 <  4:25373/149‑1 (MQ=255)
                                        aTCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGcc                                                                                                              >  3:86551/1‑149 (MQ=255)
                                        aTCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGcc                                                                                                              >  8:45869/1‑149 (MQ=255)
                                         tCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGccc                                                                                                             >  3:136970/1‑149 (MQ=255)
                                         tCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGccc                                                                                                             <  4:157471/149‑1 (MQ=255)
                                         tCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGc                                                                                                               <  2:283807/147‑1 (MQ=255)
                                         tCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGc                                                                                                               >  1:283807/1‑147 (MQ=255)
                                           gAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCtt                                                                                                           >  4:292394/1‑149 (MQ=255)
                                            aaGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTc                                                                                                          <  4:86551/149‑1 (MQ=255)
                                            aaGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTc                                                                                                          >  2:103018/1‑149 (MQ=255)
                                            aaGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTc                                                                                                          <  6:3392/149‑1 (MQ=255)
                                            aaGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTc                                                                                                          <  7:215737/149‑1 (MQ=255)
                                                       gtgtGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGa                                                                                                                                          <  5:273733/106‑1 (MQ=255)
                                                       gtgtGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGa                                                                                                                                          >  6:273733/1‑106 (MQ=255)
                                                          tgCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGc                                                                                            >  2:347691/1‑149 (MQ=255)
                                                            ccTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGt                                                                                          <  1:309926/149‑1 (MQ=255)
                                                            ccTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGt                                                                                          <  1:347691/149‑1 (MQ=255)
                                                            ccTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCCGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGt                                                                                          <  5:355995/149‑1 (MQ=255)
                                                                gCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCTGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGtgcat                                                                                      <  3:85320/149‑1 (MQ=255)
                                                                gCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGtgca                                                                                       <  2:234979/148‑1 (MQ=255)
                                                                  aTTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGtgcatg                                                                                     >  3:244627/1‑148 (MQ=255)
                                                                  aTTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGtgcatg                                                                                     >  6:15922/1‑148 (MQ=255)
                                                                   ttCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGtgcatgca                                                                                   <  4:136970/149‑1 (MQ=255)
                                                                   ttCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGtgcatgca                                                                                   <  7:45869/149‑1 (MQ=255)
                                                                       tCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTtct                                                                                                                             >  7:188678/1‑103 (MQ=255)
                                                                       tCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTtct                                                                                                                             <  8:188678/103‑1 (MQ=255)
                                                                       tCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACgg                                                                                >  3:293817/1‑148 (MQ=255)
                                                                         cAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTa                                                                              >  8:123111/1‑148 (MQ=255)
                                                                         cAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTa                                                                              >  1:177699/1‑148 (MQ=255)
                                                                         cAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTa                                                                              <  1:103018/148‑1 (MQ=255)
                                                                                ttGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCg                                                                      >  8:204647/1‑149 (MQ=255)
                                                                                  gATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGt                                                                               >  8:183132/1‑138 (MQ=255)
                                                                                  gATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGt                                                                               <  7:183132/138‑1 (MQ=255)
                                                                                  gATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCgg                                                                     <  4:293817/148‑1 (MQ=255)
                                                                                  gATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCg                                                                      <  4:244627/147‑1 (MQ=255)
                                                                                    tCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTcc                                                                       <  2:340052/144‑1 (MQ=255)
                                                                                    tCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTcc                                                                       >  1:340052/1‑144 (MQ=255)
                                                                                    tCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAg                                                                  >  1:324605/1‑149 (MQ=255)
                                                                                         ccgTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGc                                                             <  7:204647/149‑1 (MQ=255)
                                                                                                 cACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACgg                                                                                >  2:235725/1‑122 (MQ=255)
                                                                                                  aCGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCAc                                                    <  2:324605/149‑1 (MQ=255)
                                                                                                         ccAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTgg                                                >  5:20805/1‑146 (MQ=255)
                                                                                                         ccAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTgg                                                <  6:20805/146‑1 (MQ=255)
                                                                                                          cAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCcga                                            >  6:30348/1‑149 (MQ=255)
                                                                                                          cAGGTCCGGGAAGTCCCTCATCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGcc                                              <  4:84696/147‑1 (MQ=255)
                                                                                                           aGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATcgc                                                                                                 >  4:302998/1‑95 (MQ=255)
                                                                                                           aGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATcgc                                                                                                 <  3:302998/95‑1 (MQ=255)
                                                                                                              tCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCcgacg                                          >  6:310129/1‑147 (MQ=255)
                                                                                                              tCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGaa                                        <  7:328190/149‑1 (MQ=255)
                                                                                                              tCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGaa                                        <  2:169722/149‑1 (MQ=255)
                                                                                                                cGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAg                                       <  5:310129/148‑1 (MQ=255)
                                                                                                                    aaGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCAc                                   >  8:173069/1‑148 (MQ=255)
                                                                                                                        cccTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGt                              >  3:222967/1‑149 (MQ=255)
                                                                                                                         ccTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGtt                             >  2:1551/1‑149 (MQ=255)
                                                                                                                         ccTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGt                              >  5:223017/1‑148 (MQ=255)
                                                                                                                         ccTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGt                              <  5:15922/148‑1 (MQ=255)
                                                                                                                          ctctTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTc                            <  7:123111/149‑1 (MQ=255)
                                                                                                                          ctctTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTc                            >  1:210176/1‑149 (MQ=255)
                                                                                                                             ttCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGc                         >  4:11761/1‑149 (MQ=255)
                                                                                                                               ttccgccgctgcgCTATGACATCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCg                           <  8:275624/139‑1 (MQ=255)
                                                                                                                               ttccgccgctgcgCTATGACATCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCg                           >  7:275624/7‑145 (MQ=255)
                                                                                                                                 aCTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAg                     >  4:56371/1‑149 (MQ=255)
                                                                                                                                  ctgcgctgcgATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGc                    >  3:291119/1‑149 (MQ=255)
                                                                                                                                   gccgctgcgCTATGACATCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTc                                                          <  3:280903/108‑1 (MQ=255)
                                                                                                                                   gccgctgcgCTATGACATCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTc                                                          >  4:280903/3‑110 (MQ=255)
                                                                                                                                     cgctgcgATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGCt                   <  5:30348/147‑1 (MQ=255)
                                                                                                                                        tgcgATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGCCGAAGTCACTCAGTTCGGCGTAGCTGTCAc              >  8:103809/1‑149 (MQ=255)
                                                                                                                                        tgcgATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGCTGTCAc              >  7:215388/1‑149 (MQ=255)
                                                                                                                                           gATCTCATATCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTgg                                                <  5:210003/112‑1 (MQ=255)
                                                                                                                                           gATCTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTgg                                                >  6:210003/1‑112 (MQ=255)
                                                                                                                                            aTGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCgg                                                                     >  3:226617/1‑90 (MQ=255)
                                                                                                                                            aTGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCgg                                                                     <  4:226617/90‑1 (MQ=255)
                                                                                                                                               tCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGCTGTCACGGCCGt        <  3:292394/148‑1 (MQ=255)
                                                                                                                                                    gcgcAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAg                                       >  8:13054/1‑112 (MQ=255)
                                                                                                                                                    gcgcAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAg                                       <  7:13054/112‑1 (MQ=255)
                                                                                                                                                    gcgcAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGCt                   >  2:119558/1‑132 (MQ=255)
                                                                                                                                                    gcgcAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGCTGTCACGGCCGTTGATAc  >  8:245581/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATAAGGGTTGTGATACGCCTCGCGCCCCAGCATCACGCCATCGAAGGTTTCAAGGTGTGCCTGGCATTCGTCCAGGGTCTTGATCCCGCCGTTGAGCACGATCTCCAGGTCCGGGAAGTCCCTCTTCAACTGCGCTGCGATGTCATAGCGCAGCGGCGGAATCTCGCGGTTCTCTTTCGGCGACAGCCCTTCCAGAATCGCGATGCGTGCATGCACGGTAAAACTCCGGCAGCCGGCCTCACGCACCTGGCCGACGAAGTCACTCAGTTCGGCGTAGCTGTCACGGCCGTTGATAC  >  NC_002947/2475624‑2475920

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