Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,725 G→A 7.2% intergenic (+14/+85) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,243,7250GA 7.2% 250.9 / 4.4 84intergenic (+14/+85)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base G (36/41);  new base A (3/3);  total (39/45)
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

TGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGC  >  NC_002947/1243590‑1243873
                                                                                                                                       |                                                                                                                                                     
tGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTTCCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGc                                                                                                                                            <  1:312375/147‑1 (MQ=255)
     ttCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTg                                                                                                                                     >  5:424372/1‑149 (MQ=255)
        cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGcc                                                                                                                                  <  3:254744/149‑1 (MQ=255)
               agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttt                                                                                                                           <  6:424372/149‑1 (MQ=255)
                gCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtttttt                                                                                                                          >  7:213474/1‑149 (MQ=255)
                  aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTa                                                                                                                        <  1:109116/149‑1 (MQ=255)
                  aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGttttttca                                                                                                                        <  1:210282/149‑3 (MQ=255)
                   cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                       <  6:156043/149‑1 (MQ=255)
                   cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                       >  1:56989/1‑149 (MQ=255)
                   cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                       <  8:106882/149‑1 (MQ=255)
                   cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTCGCCGGTTTTTTTAt                                                                                                                       <  7:351708/149‑1 (MQ=255)
                    ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATg                                                                                                                      >  5:298512/1‑149 (MQ=255)
                    ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATg                                                                                                                      >  2:328185/1‑149 (MQ=255)
                      cGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGcc                                                                                                                    >  5:409807/1‑149 (MQ=255)
                      cGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAACGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGACCCAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGca                                                                                                                    <  5:267898/149‑2 (MQ=38)
                          gtgtGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGAATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTcc                                                                                                                >  6:268094/1‑149 (MQ=255)
                           tgtgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑t                                                                                                              <  8:213474/149‑1 (MQ=255)
                               cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                       >  1:137746/1‑137 (MQ=38)
                               cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                       <  2:142442/137‑1 (MQ=38)
                               cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                       <  2:137746/137‑1 (MQ=38)
                               cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                       >  1:142442/1‑137 (MQ=38)
                                 cTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCg                                                                                                        <  7:111809/149‑1 (MQ=255)
                                     gaaCTGGAAAAAAGGCGTCAAGGGCATGATCGCTTCGCCGGATGGCGTTGCTGCCTACCTGTGTGAGAACTCTGGCTAGCTGTATTTGGTGCGTGCTTCATTAACCCGTCCTTCTTGGCCGGTTTTTTTATGCCTTCG‑TGTGCCGGCct                                                                                                    <  4:258930/148‑1 (MQ=255)
                                       aCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCctct                                                                                                  >  6:94987/1‑149 (MQ=255)
                                         tGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTc                                                                                                >  2:54080/1‑149 (MQ=255)
                                            aaaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTcgcg                                                                                             >  3:284148/1‑149 (MQ=255)
                                             aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCgg                                                                                            >  1:74431/1‑149 (MQ=255)
                                              aaaaGGCGACAAGGGCATGAAAGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCggg                                                                                           >  4:266018/1‑149 (MQ=255)
                                                 aGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGgcgc                                                                                        <  5:268094/149‑1 (MQ=255)
                                                   gCGACAACGGCATGAACGCTTCGCCCGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGgcgcgc                                                                                      <  2:18784/149‑1 (MQ=255)
                                                     gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                    >  5:416653/1‑149 (MQ=255)
                                                       cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGgcgcgc                                                                                      <  3:172321/145‑1 (MQ=255)
                                                       cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGATGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGgcgcgc                                                                                      >  4:172321/1‑145 (MQ=255)
                                                          gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGaa                                                                               <  5:94987/149‑1 (MQ=255)
                                                            gCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATg                                                                             >  4:293268/1‑149 (MQ=255)
                                                            gCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATg                                                                             >  5:351783/1‑149 (MQ=255)
                                                             cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGc                                                                            <  2:56989/149‑1 (MQ=255)
                                                             cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTAACTGAGCGACAACGCGGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATg                                                                             <  1:54080/149‑1 (MQ=255)
                                                                gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                         >  6:320222/1‑149 (MQ=255)
                                                                 aacccgtCGCCTGAAGGCGTTGCTGCCTACCTGAGCGACAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCttt                                                                        <  7:254784/143‑1 (MQ=255)
                                                                 aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCAAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTc                                                                                                                                        >  6:336770/1‑86 (MQ=255)
                                                                    gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGc                                                                     <  8:35172/149‑1 (MQ=255)
                                                                         gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATc                                                                <  3:293268/149‑1 (MQ=255)
                                                                         gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATc                                                                <  6:409807/149‑1 (MQ=255)
                                                                            tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGtttt                                                                                                                            <  3:47270/87‑1 (MQ=25)
                                                                            tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGtttt                                                                                                                            >  4:47270/1‑87 (MQ=25)
                                                                                 gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCGTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                         >  5:424375/1‑148 (MQ=255)
                                                                                   gTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                              <  6:172553/125‑1 (MQ=255)
                                                                                   gTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                              >  5:172553/1‑125 (MQ=255)
                                                                                       ctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATc                                                   <  8:387028/148‑1 (MQ=255)
                                                                                          ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCgcgc                                               >  1:208055/1‑149 (MQ=255)
                                                                                          ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCgcgc                                               <  2:74431/149‑1 (MQ=255)
                                                                                            tACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCgcgc                                               >  3:295148/1‑147 (MQ=255)
                                                                                               cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc                                          >  2:234462/1‑149 (MQ=255)
                                                                                                tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                    <  2:19451/122‑1 (MQ=255)
                                                                                                tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccg                                         >  3:113976/1‑149 (MQ=255)
                                                                                                    cGAGAACGCTGGCAAGCTGTAATTGCTGTGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCa                                     >  5:138761/1‑149 (MQ=255)
                                                                                                    cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCa                                     <  1:234462/149‑1 (MQ=255)
                                                                                                     gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                    <  6:416653/149‑1 (MQ=255)
                                                                                                     gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                    >  8:232038/1‑149 (MQ=255)
                                                                                                     gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCACCCAt                                    >  6:316934/1‑149 (MQ=255)
                                                                                                      agaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTATAAGATAAATCGTTGAAATCGCTCCAGCCCGCCAtt                                   >  8:304971/1‑149 (MQ=255)
                                                                                                      agaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtt                                   >  7:97910/1‑149 (MQ=255)
                                                                                                      agaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtt                                   <  6:424375/149‑1 (MQ=255)
                                                                                                       gaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAttt                                  <  7:232038/149‑1 (MQ=255)
                                                                                                         aCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttct                                <  4:295148/149‑1 (MQ=255)
                                                                                                            cTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                             <  4:113976/149‑1 (MQ=255)
                                                                                                            cTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCACCCAtttctttc                             <  5:316934/149‑1 (MQ=255)
                                                                                                             tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGc                                                                            >  2:352423/1‑101 (MQ=255)
                                                                                                             tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGc                                                                            <  1:352423/101‑1 (MQ=255)
                                                                                                               gCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                           <  5:320222/148‑1 (MQ=255)
                                                                                                                    cTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGAGCCAGCCGCCCATTTCTTTCCAGCGGtt                     <  6:351783/149‑1 (MQ=255)
                                                                                                                     tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                     >  2:192479/1‑148 (MQ=255)
                                                                                                                      gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                   >  6:7308/1‑149 (MQ=255)
                                                                                                                       tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGAc                  >  6:172256/1‑149 (MQ=255)
                                                                                                                         aTTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                    <  3:172363/81‑1 (MQ=255)
                                                                                                                         aTTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                    >  4:172363/1‑81 (MQ=255)
                                                                                                                             ctgcTTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc             <  7:304971/148‑1 (MQ=255)
                                                                                                                             ctgcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                             <  7:163364/100‑1 (MQ=255)
                                                                                                                             ctgcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                             >  8:163364/1‑100 (MQ=255)
                                                                                                                                cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCg           <  5:172256/147‑1 (MQ=255)
                                                                                                                                cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCa         <  2:208055/149‑1 (MQ=255)
                                                                                                                                    cATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAAc     <  6:138761/149‑1 (MQ=255)
                                                                                                                                     aTGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                               >  1:161028/1‑122 (MQ=255)
                                                                                                                                       gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                               <  5:20454/120‑1 (MQ=255)
                                                                                                                                       gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                               >  6:20454/1‑120 (MQ=255)
                                                                                                                                       gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGc  >  5:306627/1‑149 (MQ=255)
                                                                                                                                       gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCG‑ACGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCAGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGc  >  2:81321/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                                     
TGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCC‑TGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGC  >  NC_002947/1243590‑1243873

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