Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,279 C→G 6.2% K569N (AAG→AAC)  PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,2790CG 6.2% 257.3 / 6.5 81K569N (AAG→AAC) PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base G (3/2);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

GGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATT  >  NC_002947/5539134‑5539427
                                                                                                                                                 |                                                                                                                                                      
ggTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtc                                                                                                                                                     <  8:188152/149‑1 (MQ=255)
  tAGAACTCCTGTACCTGGTCGGCGCGGACGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATGATGTCATCGTCGGTGGTACCGGGGTGCGCGGCTATGAAGTACGGGAGCAGGGACTGCTCCTTGCCCGCGTCTTTCGAGTACTtctc                                                                                                                                                   >  2:155317/1‑149 (MQ=255)
  tAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                                   <  6:113875/149‑1 (MQ=255)
  tAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                                   <  7:61292/149‑1 (MQ=255)
     aaCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCGGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGaa                                                                                                                                                >  6:109223/1‑149 (MQ=255)
      aCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGaa                                                                                                                                                >  2:355045/1‑148 (MQ=255)
          cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaagaag                                                                                                                                           <  7:79522/149‑6 (MQ=255)
          cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgc                                                                                                                                           >  3:120069/1‑149 (MQ=255)
             tACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgcgc                                                                                                                                         <  8:279934/148‑1 (MQ=255)
             tACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCt                                                                                                                                        <  7:315494/149‑1 (MQ=255)
              aCCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTACGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCtt                                                                                                                                       >  1:100858/1‑149 (MQ=255)
               ccTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCtt                                                                                                                                       >  5:207092/1‑148 (MQ=255)
                   gTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGGGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGaa                                                                                                                                    <  2:276818/147‑1 (MQ=255)
                    tCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCg                                                                                                                                 >  1:151855/1‑149 (MQ=255)
                       gcgcGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTc                                                                                                                              <  2:343795/149‑1 (MQ=255)
                           ggAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAg                                                                                                                          <  7:57327/149‑1 (MQ=255)
                           ggAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAg                                                                                                                          <  6:207092/149‑1 (MQ=255)
                            gAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAgg                                                                                                                         <  5:132423/149‑1 (MQ=255)
                             aaGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGt                                                                                                                        >  6:165275/1‑149 (MQ=255)
                              aGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTg                                                                                                                       >  4:26342/1‑149 (MQ=255)
                                  gTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGa                                                                                                                   >  4:181765/1‑149 (MQ=255)
                                      cccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCGGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATAcc                                                                                                               <  5:109223/149‑1 (MQ=255)
                                       ccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATAccc                                                                                                              <  4:19038/149‑1 (MQ=255)
                                       ccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATAc                                                                                                                <  6:124575/147‑1 (MQ=255)
                                       ccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATAc                                                                                                                >  5:124575/1‑147 (MQ=255)
                                                acaGGGCCAGGTTCATCAGGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtcatc                                                                                                     <  4:301330/149‑1 (MQ=255)
                                                     gCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCtt                                                                                                >  2:168811/1‑149 (MQ=255)
                                                     gCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCtt                                                                                                <  7:78766/149‑1 (MQ=255)
                                                     gCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCtt                                                                                                >  3:81420/1‑149 (MQ=255)
                                                          gTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCa                                                                                           <  3:26342/149‑1 (MQ=255)
                                                              atcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCg                                                                                                             >  4:343934/1‑127 (MQ=255)
                                                              atcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCg                                                                                                             <  3:343934/127‑1 (MQ=255)
                                                              atcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑ggg                                                                                       <  4:120069/148‑1 (MQ=255)
                                                               tcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTg                                                                                                                       >  2:361590/1‑116 (MQ=255)
                                                               tcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTg                                                                                                                       <  1:361590/116‑1 (MQ=255)
                                                                 atGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGcc                                                                                   >  6:236294/1‑149 (MQ=255)
                                                                   gTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCAc                                                                                 >  8:155271/1‑149 (MQ=255)
                                                                        cgtcgGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGc                                                                               <  2:318853/146‑1 (MQ=255)
                                                                        cgtcgGGGGTACCGGCGTGCGCGGCAATGAAGTACGGGATCAGGCACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCg                                                                            <  8:72065/149‑1 (MQ=255)
                                                                                 tACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTcc                                                                   >  6:339007/1‑149 (MQ=255)
                                                                                   ccGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCgg                                                                 >  6:201043/1‑149 (MQ=255)
                                                                                    cGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCggg                                                                <  4:81420/149‑1 (MQ=255)
                                                                                     ggggTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCgggg                                                               <  2:100858/149‑1 (MQ=255)
                                                                                       ggTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGc                                                              >  2:95444/1‑148 (MQ=255)
                                                                                            gcgGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCtt                                                       >  7:26598/1‑149 (MQ=255)
                                                                                              ggCAATTAAGTACGGGATCAGCTACAGCTCCTTGCCCGCTTCTTTAGAGAACGTCTCGAACAAGCGCTTGAAGCGGTCGTAGGTGACGATACCCGGCTTCATCATCTTGTCCA‑GGGCGGCCCGCTCGCTGTCGTGCGGGGCA‑ATCtt                                                       >  4:229105/1‑147 (MQ=255)
                                                                                              gccAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCa                                                     <  8:132802/147‑1 (MQ=255)
                                                                                               gCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAg                                                    <  1:168811/149‑1 (MQ=255)
                                                                                               gCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAg                                                    <  8:100518/149‑1 (MQ=255)
                                                                                                 aaTGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGTAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAgg                                                   >  2:32934/1‑148 (MQ=255)
                                                                                                   tGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGTAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAg                                                <  1:32934/149‑1 (MQ=255)
                                                                                                     aaGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAgcc                                              <  5:165275/149‑1 (MQ=255)
                                                                                                       gTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCtag                                                                            >  8:249248/1‑116 (MQ=255)
                                                                                                       gTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCg                                                                            <  7:249248/118‑1 (MQ=255)
                                                                                                            ggATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACg                                       >  3:206343/1‑149 (MQ=255)
                                                                                                             gATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACgt                                      >  3:12421/1‑149 (MQ=255)
                                                                                                                  ggTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGtggtgg                                 <  7:155271/149‑1 (MQ=255)
                                                                                                                       tGCTCCTTTCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGAc                            >  6:357738/1‑149 (MQ=255)
                                                                                                                           ccTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTc                                                      <  1:188245/119‑1 (MQ=255)
                                                                                                                           ccTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTc                                                      >  2:188245/1‑119 (MQ=255)
                                                                                                                            cTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCACATCTTCAGGTAGCCGCCAACGTGGTGGGTGAc                            <  1:44262/145‑1 (MQ=255)
                                                                                                                            cTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGAc                            >  2:44262/1‑144 (MQ=255)
                                                                                                                            cTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGtt                       >  7:181043/1‑149 (MQ=255)
                                                                                                                               gCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTc                      <  4:206343/147‑1 (MQ=255)
                                                                                                                                cccGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCt                    <  3:181765/148‑1 (MQ=255)
                                                                                                                                 ccGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCtt                   >  8:342143/1‑148 (MQ=255)
                                                                                                                                   gCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGtggtgg                                 <  5:21814/132‑1 (MQ=255)
                                                                                                                                   gCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGtggtgg                                 >  6:21814/1‑132 (MQ=255)
                                                                                                                                   gCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGAc                <  1:355045/149‑1 (MQ=255)
                                                                                                                                   gCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGAc                <  5:339007/149‑1 (MQ=255)
                                                                                                                                   cgcatgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGAc                >  5:111150/7‑149 (MQ=255)
                                                                                                                                    gcatgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCC‑GGGCA‑ATGGTAAGGTAGGCGCCAACGTTGTGGGTGGACAGTTGCGTGACg               >  2:328248/6‑148 (MQ=255)
                                                                                                                                      atgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTggtgggt                               >  6:58763/4‑131 (MQ=255)
                                                                                                                                      atgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTggtgggt                               <  5:58763/128‑1 (MQ=255)
                                                                                                                                         ttCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACACGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTc          >  5:199756/1‑149 (MQ=255)
                                                                                                                                           cGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCg        >  5:221933/1‑149 (MQ=255)
                                                                                                                                            gagaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACgtg                                     >  7:327391/1‑120 (MQ=255)
                                                                                                                                            gagaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACgtg                                     <  8:327391/120‑1 (MQ=255)
                                                                                                                                            gagaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCgg       >  4:343458/1‑149 (MQ=255)
                                                                                                                                            gagaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCgg       >  4:194180/1‑149 (MQ=255)
                                                                                                                                                aCTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGAt   >  8:70354/1‑149 (MQ=255)
                                                                                                                                                 cTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGAtt  >  1:327186/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
GGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCA‑GGGGGCCACGCTCGGTGTGCTCCGGGGCA‑ATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATT  >  NC_002947/5539134‑5539427

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