Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,703,252 T→G 5.1% L199L (CTA→CTC dppD ← dipeptide/heme ABC transporter ATPas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,703,2520TG 5.1% 345.3 / 10.0 98L199L (CTA→CTCdppDdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base T (44/49);  new base G (3/2);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

ACGGGTGACGCGGCGCATGGAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCGGG  >  NC_000913/3703104‑3703399
                                                                                                                                                    |                                                                                                                                                   
aCGGGTGACGCGGCGCATGGAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGt                                                                                                                                                     <  4:760033/149‑1 (MQ=255)
  gggTGACGCGGCGCATGGAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAg                                                                                                                                                   >  4:5739/1‑149 (MQ=255)
        cgcgGCGCATGGAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCca                                                                                                                                              <  1:326476/148‑1 (MQ=255)
                 tGGAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTc                                                                                                                                    >  4:536808/1‑149 (MQ=255)
                  ggAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTc                                                                                                                                    >  2:620393/1‑148 (MQ=255)
                   gAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCga                                                                                                                                  <  3:89171/149‑1 (MQ=255)
                    aaGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGTAGTAGTTCgat                                                                                                                                 <  3:720468/149‑1 (MQ=255)
                    aaGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCga                                                                                                                                  >  1:248337/1‑148 (MQ=255)
                    aaGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCga                                                                                                                                  <  2:330907/148‑1 (MQ=255)
                    aaGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCga                                                                                                                                  <  1:813799/148‑1 (MQ=255)
                       aTGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCgatga                                                                                                                               >  2:496563/1‑148 (MQ=255)
                        tGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCgatgat                                                                                                                              <  1:393828/148‑1 (MQ=255)
                        tGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGAtt                                                                                                                             <  4:829625/149‑1 (MQ=255)
                           cGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTgc                                                                                                                          <  4:4586/149‑1 (MQ=255)
                              gcgcATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTgcgc                                                                                                                        >  1:235794/1‑148 (MQ=255)
                               cgcATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCt                                                                                                                      <  4:525522/149‑1 (MQ=255)
                                gcATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTg                                                                                                                     >  3:111701/1‑149 (MQ=255)
                                 cATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTg                                                                                                                     >  1:785117/1‑148 (MQ=255)
                                  aTCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGaa                                                                                                                   <  3:500926/149‑1 (MQ=255)
                                    cACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAAt                                                                                                                  >  2:649940/1‑148 (MQ=255)
                                     aCCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATgg                                                                                                                >  4:379075/1‑149 (MQ=255)
                                      ccGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATgg                                                                                                                >  2:841182/1‑148 (MQ=255)
                                      ccGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATgg                                                                                                                >  1:180294/1‑148 (MQ=255)
                                       cGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGt                                                                                                               <  2:559898/148‑1 (MQ=255)
                                        ggTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTgc                                                                                                                          <  2:493132/136‑1 (MQ=255)
                                        ggTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTgc                                                                                                                          >  1:493132/1‑136 (MQ=255)
                                             ccaccaCCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTcc                                                                                                        <  4:855140/149‑1 (MQ=255)
                                               accaccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAg                                                                                                      >  4:244013/1‑149 (MQ=255)
                                                ccaccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCga                                                                                                                                  <  1:545868/120‑1 (MQ=255)
                                                ccaccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCga                                                                                                                                  >  2:545868/1‑120 (MQ=255)
                                                ccaccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAg                                                                                                      >  1:613384/1‑148 (MQ=255)
                                                ccaccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAg                                                                                                      >  2:35408/1‑148 (MQ=255)
                                                ccaccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAg                                                                                                      >  2:270613/1‑148 (MQ=255)
                                                  accTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAgcg                                                                                                    >  1:460613/1‑148 (MQ=255)
                                                   ccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATgg                                                                                                                <  2:101714/135‑1 (MQ=255)
                                                   ccTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATgg                                                                                                                >  1:101714/1‑135 (MQ=255)
                                                     tGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCgg                                                                                                 <  1:270667/148‑1 (MQ=255)
                                                       gCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCggtg                                                                                               <  1:496563/148‑1 (MQ=255)
                                                           gCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCg                                                                                           <  2:248337/148‑1 (MQ=255)
                                                           gCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCg                                                                                           <  1:87725/148‑1 (MQ=255)
                                                             aTACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGtt                                                                                        <  3:181070/149‑1 (MQ=255)
                                                              tACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGtt                                                                                        <  2:722274/148‑1 (MQ=255)
                                                                  tCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATc                                                                                    >  2:558060/1‑148 (MQ=255)
                                                                    aCGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGc                                                                                 >  3:300360/1‑149 (MQ=255)
                                                                     cGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGc                                                                                 <  2:460613/148‑1 (MQ=255)
                                                                         gatTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATca                                                                            <  4:756363/149‑1 (MQ=255)
                                                                            tttATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATcagc                                                                          >  1:102165/1‑148 (MQ=255)
                                                                            tttATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATcagc                                                                          >  1:757039/1‑148 (MQ=255)
                                                                            tttATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATcagc                                                                          <  1:649940/148‑1 (MQ=255)
                                                                              tATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATcagcag                                                                        <  2:613384/148‑1 (MQ=255)
                                                                                    ccGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTgg                                                                  <  2:716608/148‑1 (MQ=255)
                                                                                    ccGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGc                                                                 >  3:855439/1‑149 (MQ=255)
                                                                                        ttCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCggtg                                                                                               >  3:523149/1‑115 (MQ=255)
                                                                                        ttCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCggtg                                                                                               <  4:523149/115‑1 (MQ=255)
                                                                                          ccGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACa                                                            <  1:694494/148‑1 (MQ=255)
                                                                                           cGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAgg                                                          <  3:916623/149‑1 (MQ=255)
                                                                                             ccaccaGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGc                                                         >  2:804904/1‑148 (MQ=255)
                                                                                                ccaGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAAt                                                      <  2:455828/148‑1 (MQ=255)
                                                                                                   gcgcCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGc                                                   >  2:305142/1‑148 (MQ=255)
                                                                                                    cgcCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGcc                                                  <  2:235794/148‑1 (MQ=255)
                                                                                                       cAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATgg                                              >  3:251307/1‑149 (MQ=255)
                                                                                                       cAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATg                                               <  1:305142/148‑1 (MQ=255)
                                                                                                          atCAAGGGTAATTAACACCCGCGCCATGTTCTCATTCGACTGTAGCTCCAGCAGTAGTTCTATGTTTTGCGCCTGACTGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCGGCAGTTTTGGCCCACAAGCAATCGCCATGGCg                                            <  2:693198/147‑1 (MQ=255)
                                                                                                           tCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCAc                                                                                                            >  3:102020/1‑83 (MQ=255)
                                                                                                           tCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATc                                                                             >  2:120438/1‑114 (MQ=255)
                                                                                                           tCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATc                                                                             <  1:120438/114‑1 (MQ=255)
                                                                                                           tCATGGGTAATTAACACCAGCGCCATGTCCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCAc                                                                                                            <  4:102020/83‑1 (MQ=255)
                                                                                                             aTGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGatc                                         <  2:719174/148‑1 (MQ=255)
                                                                                                               gggTAATTCACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGatcatc                                      >  4:83696/1‑149 (MQ=255)
                                                                                                               gggTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCAc                                                                                                            >  4:864545/1‑79 (MQ=255)
                                                                                                               gggTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCAc                                                                                                            <  3:864545/79‑1 (MQ=255)
                                                                                                               gggTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGatcat                                       <  1:804904/148‑1 (MQ=255)
                                                                                                                ggTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGAtcatca                                     >  3:18645/1‑149 (MQ=255)
                                                                                                                   aaTTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCAcgc                                  <  3:480069/149‑1 (MQ=255)
                                                                                                                     ttAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACgcgc                                >  4:894350/1‑149 (MQ=255)
                                                                                                                         cacCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGc                            <  3:83696/149‑1 (MQ=255)
                                                                                                                          acCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATcagca                                                                         >  3:191159/1‑103 (MQ=255)
                                                                                                                          acCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATcagca                                                                         <  4:191159/103‑1 (MQ=255)
                                                                                                                          acCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCt                           >  3:386154/1‑149 (MQ=255)
                                                                                                                              gcgcCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCAt                        >  1:434820/1‑148 (MQ=255)
                                                                                                                              gcgcCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATg                       >  4:559139/1‑149 (MQ=255)
                                                                                                                                  cATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATgccgc                   >  3:397782/1‑149 (MQ=255)
                                                                                                                                    tGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATgccgcc                  <  1:73890/148‑1 (MQ=255)
                                                                                                                                         aactactgctgGAGCTACAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCAcgc                                  >  4:304815/6‑127 (MQ=255)
                                                                                                                                         aactactgctgGAGCTACAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCAcgc                                  <  3:304815/122‑1 (MQ=255)
                                                                                                                                         aactactgctgGAGCTACAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGaaa             >  1:642710/6‑148 (MQ=255)
                                                                                                                                          ctTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAg            >  2:879975/1‑148 (MQ=255)
                                                                                                                                          ctTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAg            >  1:375543/1‑148 (MQ=255)
                                                                                                                                          ctTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGc           <  3:244013/149‑1 (MQ=255)
                                                                                                                                           tttCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGc           >  2:399891/1‑148 (MQ=255)
                                                                                                                                             actgctgGAGCTACAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAg            <  1:252799/144‑1 (MQ=255)
                                                                                                                                             actgctgGAGCTACAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAg            >  2:252799/2‑145 (MQ=255)
                                                                                                                                                gctgTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATg      <  1:841182/148‑1 (MQ=255)
                                                                                                                                                  tgTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCgg   <  4:251307/149‑1 (MQ=255)
                                                                                                                                                  tgTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCgg   <  3:946555/149‑1 (MQ=255)
                                                                                                                                                  tgTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCgg   <  3:536808/149‑1 (MQ=255)
                                                                                                                                                  tgTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCg    <  1:483571/148‑1 (MQ=255)
                                                                                                                                                  tgTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCg    <  1:715208/148‑1 (MQ=255)
                                                                                                                                                   gTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCggg  >  3:916819/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGGGTGACGCGGCGCATGGAAGATGGCGTGCGCATCACCGGTTTCCACCACCTGGCCTGCATACATCACGATGATTTTATGTGCCGCTTCCGCCACCAGCGCCAGGTCATGGGTAATTAACACCAGCGCCATGTTCTCTTTCTGCTGTAGCTCCAGCAGTAGTTCGATGATTTGCGCCTGAATGGTCACGTCCAGCGCGGTGGTCGGTTCATCGGCAATCAGCAGTTTTGGCCGACAGGCAATCGCCATGGCGATCATCACGCGCTGGCTCATGCCGCCGGAAAGCTGATGCGGG  >  NC_000913/3703104‑3703399

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