Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,649,989 C→A 5.9% G309C (GGC→TGC)  PP_4113 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,649,9890CA 5.9% 270.4 / 6.0 85G309C (GGC→TGC) PP_4113hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base A (2/3);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTC  >  NC_002947/4649843‑4650136
                                                                                                                                                   |                                                                                                                                                     
gcgcgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGccc                                                                                                                                                      <  8:307864/148‑1 (MQ=255)
gcgcgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGcc                                                                                                                                                       >  3:410343/1‑147 (MQ=255)
     cgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCGGGGCTGGGCAAGCACAGcc‑ct                                                                                                                                               <  2:9454/149‑2 (MQ=255)
       ccGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCgg                                                                                                                                             <  4:380932/149‑1 (MQ=255)
       ccGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGGGCAAGCACAGCC‑CTgg                                                                                                                                             <  7:20345/149‑1 (MQ=255)
          aCCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCa                                                                                                                                           >  7:367449/1‑148 (MQ=255)
           cccGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCaaa                                                                                                                                         <  1:323466/149‑1 (MQ=255)
             ccgTGAATGCCCCAGCTCCCGTTCCGCACCCGCACGCCATCTTCGACCAGGTGCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCaaaaa                                                                                                                                       <  5:327478/147‑1 (MQ=255)
             cGGGGAATGCCCCAGCTCCAGTTCAGCAGCGGCACGCCATCTTCGACGAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACG‑ATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCC‑CTGGcttcga                                                                                                                                       <  6:340467/149‑6 (MQ=255)
                  aaTGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCAATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑cc                                                                                                                                   <  1:95420/148‑1 (MQ=255)
                    tGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCa                                                                                                                               >  4:355624/1‑149 (MQ=255)
                       cccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCAccc                                                                                                                            >  2:76907/1‑149 (MQ=255)
                       cccAGATACAGTTCAGCACCCGCAAGCCATCTTCGACCAGGTTCACCGAGGCAGAGATGTTGGAGGTGATACCGGCATCCGAGTTGAGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTTCTTGCCCAGCC‑ACGGCAAAAAG‑CCCTCAcc                                                                                                                             >  2:231031/1‑148 (MQ=255)
                        ccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCAACGAGGAAGCGATGTTGGCGGTGATACCGGCATCCGATTTGCGCTCTACC‑ATCACCTCGAAGCCCGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑ACCTCACACt                                                                                                                           >  1:117543/1‑149 (MQ=255)
                              ccAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAg                                                                                                                     >  7:191943/1‑149 (MQ=255)
                               cAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGc                                                                                                                    >  2:416155/1‑149 (MQ=255)
                                aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCg                                                                                                                   <  7:376606/149‑1 (MQ=255)
                                aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCg                                                                                                                   <  8:191943/149‑1 (MQ=255)
                                    cAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGaa                                                                                                               >  3:364908/1‑149 (MQ=255)
                                     aGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAAt                                                                                                              <  3:152778/149‑1 (MQ=255)
                                     aGCACCCGCACGCCATCTTCGACCAGCTTCACCGAGTCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAC‑CCCTCACCCTGCGCAGCgagact                                                                                                              >  5:62293/1‑147 (MQ=255)
                                       cACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACA‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGc                                                                                                            <  5:367779/149‑1 (MQ=255)
                                          ccGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                         >  2:306037/1‑149 (MQ=255)
                                          ccGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                         >  5:291365/1‑149 (MQ=255)
                                                 ccATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACAt                                                                                                   <  3:175014/148‑1 (MQ=255)
                                                 ccATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACAt                                                                                                   >  4:175014/1‑148 (MQ=255)
                                                  cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATgg                                                                                                 >  2:326656/1‑149 (MQ=255)
                                                  cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATgg                                                                                                 >  4:409273/1‑149 (MQ=255)
                                                   aTCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGc                                                                                                >  7:280068/1‑149 (MQ=255)
                                                      ttCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGc                                                                                             >  5:156240/1‑149 (MQ=255)
                                                          aCCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTg                                                                                         <  4:4913/149‑1 (MQ=255)
                                                           ccAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGATGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                         >  7:32795/1‑132 (MQ=255)
                                                           ccAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGATGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                         <  8:32795/132‑1 (MQ=255)
                                                            cAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGc                                                                                       >  5:116838/1‑149 (MQ=255)
                                                               gTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTg                                                                                    >  3:292721/1‑149 (MQ=255)
                                                                ttCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGTGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGc                                                                                   <  8:123857/149‑1 (MQ=255)
                                                                    ccGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAGCCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCgg                                                                                <  1:326656/149‑1 (MQ=255)
                                                                      gAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                             <  5:89724/149‑1 (MQ=255)
                                                                      gAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                             <  1:76907/149‑1 (MQ=255)
                                                                      gAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                             <  3:355624/149‑1 (MQ=255)
                                                                          cAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTa                                                                         <  4:410343/149‑1 (MQ=255)
                                                                           aGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAc                                                                        >  8:395881/1‑149 (MQ=255)
                                                                           aGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAc                                                                        <  6:158692/149‑1 (MQ=255)
                                                                             cGATGTTGGCGGTGATACCCGCATCCGAGTTGCGCACTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGAGACATGGCTGCCGTGGCTTGCGGGGCGGTTACCg                                                                      >  4:258610/1‑149 (MQ=255)
                                                                                tGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGACGTGGCTTGCCGGCCGGTTACCGTcg                                                                   >  6:425761/1‑149 (MQ=255)
                                                                                     gCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGaa                                                                                                               <  5:192069/100‑1 (MQ=255)
                                                                                     gCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGaa                                                                                                               >  6:192069/1‑100 (MQ=255)
                                                                                     gCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGc                                                              <  6:291365/149‑1 (MQ=255)
                                                                                      cGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTTCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCg                                                             <  7:259409/149‑1 (MQ=255)
                                                                                           aTACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTaga                                                         <  1:416155/148‑1 (MQ=255)
                                                                                            tACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTagag                                                        <  6:286739/148‑1 (MQ=255)
                                                                                               cGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                                    >  3:48131/1‑149 (MQ=255)
                                                                                                 gCATCCGAGTTGGGCTCTACC‑ATCACATCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAT‑CCAGCACCCTGGGGAGCGAGCTTGCCTGCGACATGGCTGCCGTGGCTTGCGGGGCGGTGAGCGTCGCGGGCGGAGAGGCATGt                                                  >  4:237798/1‑149 (MQ=255)
                                                                                                  cATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTc                                                 >  5:393559/1‑149 (MQ=255)
                                                                                                    tCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTg                                               <  6:174502/149‑1 (MQ=255)
                                                                                                     ccGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTgg                                              >  1:134918/1‑149 (MQ=255)
                                                                                                      cGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCTGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                             >  1:332212/1‑149 (MQ=255)
                                                                                                      cGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                             <  5:160788/149‑1 (MQ=255)
                                                                                                      cGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                             <  4:292721/149‑1 (MQ=255)
                                                                                                             cgcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGtta                                      >  7:107449/1‑148 (MQ=255)
                                                                                                             cgcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTc                                      <  4:48131/149‑1 (MQ=255)
                                                                                                              gcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTCCCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCg                                     <  2:316234/149‑1 (MQ=255)
                                                                                                              gcTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCg                                     <  6:156240/149‑1 (MQ=255)
                                                                                                                   aCC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTg                                         >  8:412607/1‑140 (MQ=255)
                                                                                                                   aCC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTg                                         <  7:412607/140‑1 (MQ=255)
                                                                                                                    cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGccc                               >  7:318892/1‑149 (MQ=255)
                                                                                                                    cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGccc                               >  1:363172/1‑149 (MQ=255)
                                                                                                                    cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcc                                <  6:393559/148‑1 (MQ=255)
                                                                                                                    cc‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGc                                 <  7:298781/147‑1 (MQ=255)
                                                                                                                     c‑ATCACCTCGAAGCGAGGGCTGTTCTTGCCGAGCG‑CCGTCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGACTGCTTCATGGCTGCCGTGTCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcgcc                              <  8:367449/149‑1 (MQ=255)
                                                                                                                        tCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCTGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCt                            <  1:231031/149‑1 (MQ=255)
                                                                                                                        tCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCt                            <  8:280068/149‑1 (MQ=255)
                                                                                                                          aCCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                           <  3:409273/148‑1 (MQ=255)
                                                                                                                           gctttttgCCGGGGCTGGGCAAGCACAGCCTCTGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                           >  1:259491/8‑148 (MQ=255)
                                                                                                                               gAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAAt                     <  5:425761/149‑1 (MQ=255)
                                                                                                                               gAAGCCAGGGATTTGCTTGCCCAGCC‑CCGGCAAATAG‑CCCTCAGCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTTGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTGGATaa                      <  6:62293/148‑1 (MQ=255)
                                                                                                                                   ccggggCTGGGCAAGCACAGCC‑CTGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt                 >  2:121712/4‑149 (MQ=255)
                                                                                                                                   ccAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt                 >  2:117593/1‑149 (MQ=255)
                                                                                                                                   ccAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGAAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCTCAGCCCCCTCGATAATGGCt                 >  4:372999/1‑149 (MQ=255)
                                                                                                                                     aGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTg                >  7:53254/1‑148 (MQ=255)
                                                                                                                                     aGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTg                >  1:190698/1‑148 (MQ=255)
                                                                                                                                       ggCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAAc             >  8:252758/1‑149 (MQ=255)
                                                                                                                                            tgCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCgg        >  5:231727/1‑149 (MQ=255)
                                                                                                                                            tgCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCgg        >  1:46362/1‑149 (MQ=255)
                                                                                                                                              cTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCggg       <  5:211283/148‑1 (MQ=255)
                                                                                                                                              cTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCggg       <  2:332212/148‑1 (MQ=255)
                                                                                                                                                 gCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGt   >  7:383097/1‑149 (MQ=255)
                                                                                                                                                   ccAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTc  >  6:442959/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACC‑ATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCC‑CCGGCAAAAAG‑CCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTC  >  NC_002947/4649843‑4650136

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