Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,739 T→C 8.9% intergenic (+28/+71) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,243,7390TC 8.9% 184.1 / 7.9 67intergenic (+28/+71)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (3/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

CCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTT  >  NC_002947/1243597‑1243882
                                                                                                                                              |                                                                                                                                                
ccACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGc                                                                                                                                            >  4:160873/1‑149 (MQ=255)
 cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGcc                                                                                                                                           >  7:46896/1‑149 (MQ=255)
     aaGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtt                                                                                                                                       <  6:156388/149‑1 (MQ=255)
      agagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttt                                                                                                                                      >  7:116490/1‑149 (MQ=255)
        agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttt                                                                                                                                    <  8:116490/149‑1 (MQ=255)
            cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                <  1:23118/149‑1 (MQ=255)
             ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATg                                                                                                                               >  6:66890/1‑149 (MQ=255)
               cGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGcc                                                                                                                             >  3:19242/1‑149 (MQ=255)
                  aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTc                                                                                                                          <  7:109461/149‑1 (MQ=255)
                      tgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCt                                                                                                                            <  2:142957/143‑1 (MQ=255)
                      tgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCt                                                                                                                            >  1:142957/1‑143 (MQ=255)
                        cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGc                                                                                                                    <  2:96040/149‑1 (MQ=255)
                        cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGc                                                                                                                    >  4:277412/1‑149 (MQ=255)
                           tGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCgg                                                                                                                 >  5:128410/1‑149 (MQ=255)
                             ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGc                                                                                                                <  8:46896/148‑1 (MQ=255)
                               aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCctc                                                                                                             >  7:84780/1‑149 (MQ=255)
                                    gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgc                                                                                                        >  2:84836/1‑149 (MQ=255)
                                    gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgc                                                                                                        >  6:51555/1‑149 (MQ=255)
                                    gAAAAAAGGCGACAAGGGCAGGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcg                                                                                                         <  4:19242/148‑1 (MQ=255)
                                         aaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcg                                                                                                   >  6:112197/1‑149 (MQ=255)
                                            gCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcgcgc                                                                                                <  6:279647/149‑1 (MQ=255)
                                              gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTGTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                              >  3:35554/1‑149 (MQ=255)
                                                  aGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGa                                                                                         <  5:51555/149‑1 (MQ=255)
                                                       aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGcc                                                                                    >  4:196759/1‑149 (MQ=255)
                                                             gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTCCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                               >  1:336341/1‑148 (MQ=255)
                                                                 cGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAAGCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                           <  6:207777/148‑1 (MQ=255)
                                                                 cGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTCCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                           <  2:336341/148‑1 (MQ=255)
                                                                  tccTGAAGGCGTTGCTGCCTCCCTGCGCGCGTACGCGGGCATGCTGTACTTGCTGCGTGCATGAAACAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATc                                                                         <  4:43573/148‑1 (MQ=255)
                                                                  gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                                       >  7:200767/1‑135 (MQ=255)
                                                                  gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                                       <  8:200767/135‑1 (MQ=255)
                                                                   ccTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATc                                                                         >  8:171849/1‑148 (MQ=255)
                                                                    cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                               >  7:239476/1‑141 (MQ=255)
                                                                    cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGAt                                                                          <  8:239476/146‑1 (MQ=255)
                                                                    cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCaa                                                                       >  7:300530/1‑149 (MQ=255)
                                                                    cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCa                                                                        >  4:315478/1‑148 (MQ=255)
                                                                       aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                     >  6:209882/1‑148 (MQ=255)
                                                                       aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                     <  8:96958/148‑1 (MQ=255)
                                                                       aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                     <  5:209882/148‑1 (MQ=255)
                                                                       aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                     <  3:315478/148‑1 (MQ=255)
                                                                         ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                  <  8:84780/149‑1 (MQ=255)
                                                                         ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                  <  3:196759/149‑1 (MQ=255)
                                                                         ggCGCTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                  <  3:155850/149‑1 (MQ=255)
                                                                          gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                                 <  5:112197/149‑1 (MQ=255)
                                                                             tttctgcCTACCTGATCTAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATTAAAAAACCGGCCTTCGTGGCCGTTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGaaa                                                              <  3:185189/146‑1 (MQ=255)
                                                                                 tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                          >  4:70597/1‑149 (MQ=255)
                                                                                  gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGAt                                                                          <  4:79112/132‑1 (MQ=255)
                                                                                   ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCgcgc                                                        <  3:160873/149‑1 (MQ=255)
                                                                                        cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc                                                   <  5:66890/149‑1 (MQ=255)
                                                                                                  aCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttct                                         >  6:26044/1‑149 (MQ=255)
                                                                                                  aCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttct                                         >  5:19616/1‑149 (MQ=255)
                                                                                                    gCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                                        <  1:84836/149‑1 (MQ=255)
                                                                                                        gCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                    <  5:26044/148‑1 (MQ=255)
                                                                                                         cAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttcttt                                       <  4:4696/144‑1 (MQ=255)
                                                                                                         cAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttcttt                                       >  3:4696/1‑144 (MQ=255)
                                                                                                              tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTATTTCCacatcat                              >  1:287619/1‑142 (MQ=255)
                                                                                                                ttaTTGCTGCGTGCATGAAAAAACCGGGGTTCGTGGCCGGTTTTTTTATGCATTCCTGTGCGGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCATTCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                            <  8:30016/146‑1 (MQ=255)
                                                                                                                tAATTGCTGCTTGCATGAAAAAACCGGCCTTCGTTGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGttta                            <  6:128410/148‑3 (MQ=255)
                                                                                                                 aaTTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTAAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGAc                           >  6:293659/1‑148 (MQ=255)
                                                                                                                  aTTGCTGCGTTCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACg                          <  4:202456/148‑1 (MQ=255)
                                                                                                                  aTTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCTGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGa                         <  3:277412/149‑1 (MQ=255)
                                                                                                                              ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                             <  2:42987/114‑1 (MQ=37)
                                                                                                                              ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                             >  1:42987/4‑117 (MQ=37)
                                                                                                                                aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCa                                                      <  3:88401/105‑1 (MQ=38)
                                                                                                                                aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCa                                                      >  4:88401/2‑106 (MQ=38)
                                                                                                                                 aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAGCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGCACAGCt          >  5:137770/1‑149 (MQ=255)
                                                                                                                                   aaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGCTGCCGCAGAACAGCt          >  7:146786/1‑147 (MQ=255)
                                                                                                                                         ggCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtt  >  7:30365/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTT  >  NC_002947/1243597‑1243882

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: