Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,909,949 C→G 9.9% intergenic (+342/+137) PP_4318 → / ← PP_4319 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,909,9490CG 9.9% 174.8 / 10.6 70intergenic (+342/+137)PP_4318/PP_4319transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base G (5/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

AACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACTGTGC  >  NC_002947/4909807‑4910094
                                                                                                                                               |                                                                                                                                                 
aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATACATGGGGCCTCAGCGCTTGACGGGCATCCATGcc                                                                                                                                             >  7:78226/1‑149 (MQ=17)
    cgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGt                                                                                                                                          >  8:150897/1‑148 (MQ=32)
    cgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                         <  2:302422/149‑1 (MQ=32)
      aCAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                       <  7:14762/149‑2 (MQ=17)
       cAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTACGGGGTTACCTCGCTACGCTCTAGCTACCGGCAGTACGCGCTGCAGGCCCCATCCATGCGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCaac                                                                                                                                      >  5:362721/1‑147 (MQ=11)
        aaaaaGTAACCAAAAAACGCTGGATCCTATCATCCGGCCCCTGCGCTCCGCTCCGGGGTTACCTCGCTACGGAGATTCACCCCGGAGTACGCGCTGCACGCCCCAGCCATGGGGCCTCAGCGCAAGCCGGGCCAACATGCCCGTCACCt                                                                                                                                     >  5:201436/1‑149 (MQ=11)
         aaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTc                                                                                                                                    >  1:32997/1‑149 (MQ=32)
            aGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCt                                                                                                                                 <  2:32997/149‑1 (MQ=32)
               aaCCAAAAAACGCTGGCTACTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATCGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCg                                                                                                                              >  2:77988/1‑149 (MQ=21)
                  cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCaa                                                                                                                           <  3:6257/149‑1 (MQ=34)
                     aaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGcc                                                                                                                        >  6:193025/1‑149 (MQ=34)
                                   tATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCaa                                                                                                                           <  3:151756/132‑1 (MQ=33)
                                   tATCATCCGGCCCCTGCGCTAAGCTCCGGGGTTACCTCGCTACGGTCTTGCTACCGGGAGTACTCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCAGTCACCTCCATACGCaa                                                                                                                           >  4:151756/1‑132 (MQ=11)
                                    atcatcCGGCCCCTGCGCTACGCTCCGGGGTCTCCCTCGCTACGGTCTTGATCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCtt                                                                                                          >  1:232956/1‑149 (MQ=255)
                                    atcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTc                                                                                                         >  8:68006/1‑149 (MQ=35)
                                    atcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTc                                                                                                         >  8:267705/1‑149 (MQ=35)
                                     tcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCtt                                                                                                          <  7:150897/147‑1 (MQ=35)
                                       atcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTg                                                                                                       <  6:189366/148‑1 (MQ=37)
                                        tcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGaa                                                                                                     <  5:193025/149‑1 (MQ=37)
                                           ggCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTc                                                                                                  <  2:23057/149‑1 (MQ=37)
                                              cccTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTc                                                                                                  <  7:143358/146‑1 (MQ=37)
                                              cccTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTAAGCGCTGCAGGCCCCCTCCAAGGGGCCTCAGCGCTTGACGGGCATCCCTGCCCCTCACCACCCGTCGCAAGCCCTGCGCGCCGCCATCGGAAGTc                                                                                                  >  8:143358/1‑146 (MQ=11)
                                                cTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCg                                                                                                               <  2:206728/131‑1 (MQ=34)
                                                cTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCg                                                                                                               >  1:206728/1‑131 (MQ=34)
                                                cTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAAt                                                                                             <  7:267705/149‑1 (MQ=37)
                                                   cgcTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCt                                                                                           <  7:68006/148‑1 (MQ=37)
                                                    gcTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                                                         <  3:224043/149‑1 (MQ=37)
                                                      tACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGt                                                                                       >  3:193316/1‑149 (MQ=37)
                                                      tACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGt                                                                                       >  8:184672/1‑149 (MQ=17)
                                                      tACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGCGGCCTCAGCGCTTGACGGGCATACATGCCCGTCACCTCCCTTCGCAAGCCATgcgcacc                                                                                                               >  6:305760/1‑122 (MQ=255)
                                                      tACGCTCCGGGGTTCCCTCGCTAAGGTCTTGCTACCGGGAGTACGCGCTGCAGGCCCCCTCCATGGGGCCTCAGCGCTTGACGGGCATACATGGCCGTGACCTCCCTTCGCAAGCCCTGCGCTCGGGCTTCTGAAGTAGCAATATGCGt                                                                                       >  5:33245/1‑149 (MQ=11)
                                                             cGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTg                                                                                 >  7:102148/1‑148 (MQ=37)
                                                              ggggTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGaa                                                                               <  1:217949/149‑1 (MQ=37)
                                                                     cccTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTa                                                                            <  2:60715/146‑1 (MQ=37)
                                                                     cccTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTa                                                                            >  1:60715/1‑146 (MQ=37)
                                                                      ccTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGc                                                                                    >  4:112972/1‑137 (MQ=37)
                                                                      ccTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGc                                                                                    <  3:112972/137‑1 (MQ=37)
                                                                       cTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCAAGGGGCCTCAGCGCTTGACAGGCATCCATGCCCGTCACCTCCCATCCCAAGCCCTGAGCTCCGCCTTCTGAAGTCGCAAACTACGTCCCCTGAACAACcgca                                                                       >  6:330996/1‑148 (MQ=14)
                                                                           cTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCt                                                                   <  6:33245/149‑1 (MQ=255)
                                                                             aCGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTAGCAATCTGCGTCGCCTGAACTACCGCGCGCTTa                                                                 >  6:60379/1‑149 (MQ=255)
                                                                              cGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCtt                                                                  <  8:102148/147‑1 (MQ=39)
                                                                                      cTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCaaa                                                         <  7:184672/148‑1 (MQ=25)
                                                                                        cccGGGAGTACGCGCTGCAGGCCCCTTCCTTGGGGCCTCAGCGCTTGTCGGTCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGGTCGGCCTTCTGAAGTCGCAATCTGCGTGGCCTGAACTACCGCGCGCTTACAAGGAAaagg                                                      <  6:201436/149‑2 (MQ=14)
                                                                                            ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaa                                                  <  8:112139/149‑1 (MQ=255)
                                                                                                  cgcgCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGaa                                                                               <  1:197888/114‑1 (MQ=35)
                                                                                                  cgcgCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGaa                                                                               >  2:197888/1‑114 (MQ=35)
                                                                                                    cgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                                           >  6:28294/1‑148 (MQ=255)
                                                                                                        gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaa                                      >  7:257759/1‑149 (MQ=255)
                                                                                                            gCCCCATCCATGGGGACTCAGCGCTTGACGGGCATCCATGCCCGTAACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag                                  >  2:31859/1‑149 (MQ=255)
                                                                                                             ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag                                  >  7:30217/1‑148 (MQ=255)
                                                                                                             ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagc                                 <  8:257759/149‑1 (MQ=255)
                                                                                                                     atGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcagaag                            >  6:146123/1‑146 (MQ=255)
                                                                                                                      tGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCagcagaagcaga                        <  5:28294/149‑1 (MQ=255)
                                                                                                                       ggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAACAATAGCAGAAGCTGAAGCAGAAGCagcagaagcaga                        >  7:174901/1‑148 (MQ=255)
                                                                                                                        gggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATa                      <  1:370874/149‑1 (MQ=255)
                                                                                                                           ccTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAg                   <  5:146123/149‑1 (MQ=255)
                                                                                                                            cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGACGCAGAGACAg                   >  8:147726/1‑148 (MQ=255)
                                                                                                                             tCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAg                   <  4:276368/147‑1 (MQ=255)
                                                                                                                             tCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAg                   >  3:276368/1‑147 (MQ=255)
                                                                                                                             tCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGa                  >  8:264282/1‑148 (MQ=255)
                                                                                                                                 gaggtGACCGGCATGGATGCCCCTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTTAAGTCGCAAAATGCGTCGCCGGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCACAAGCAGAACCAGCAGAAGCAGATAAACAAACAg             >  6:233285/5‑149 (MQ=18)
                                                                                                                                   cTTGACGGGCATCCATGCCCGTCACCTCCGTTCGCAAGCCCTGCGCTCTGCCTTCTGAAGTCGGAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAgtg           <  8:78226/149‑1 (MQ=255)
                                                                                                                                    gtGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGCACTACCGCGAGCTTACAAGCAAAAGCAGAAGCCGAAGCAGAAGCAGCAGAAGCAGATACTGAAACAgtgt          >  8:135999/2‑149 (MQ=18)
                                                                                                                                     tGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAgtgt          >  8:220730/1‑148 (MQ=255)
                                                                                                                                     tGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAgtgt          >  8:226448/1‑148 (MQ=255)
                                                                                                                                      gACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagc                                 <  8:311114/124‑1 (MQ=255)
                                                                                                                                      gACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagc                                 >  7:311114/1‑124 (MQ=255)
                                                                                                                                      gACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAgtgt          <  4:193316/147‑1 (MQ=255)
                                                                                                                                         gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATAc                     >  8:185141/1‑133 (MQ=255)
                                                                                                                                         gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATAc                     <  7:185141/133‑1 (MQ=255)
                                                                                                                                            cagccATGCCCGTCACCTCCCTTCGCTAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAACGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAACCAGTGGTACTGTGc  <  5:220566/146‑1 (MQ=21)
                                                                                                                                               |                                                                                                                                                 
AACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACTGTGC  >  NC_002947/4909807‑4910094

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