Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,186 G→T 5.8% intergenic (+48/+19) PP_5376 → / ← PP_5377 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,129,1860GT 5.8% 270.8 / 5.9 86intergenic (+48/+19)PP_5376/PP_5377hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base T (2/3);  total (42/44)
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.22e-01

TGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGGGACGGCGA  >  NC_002947/6129039‑6129332
                                                                                                                                                      |                                                                                                                                                   
tgACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTg                                                                                                                                                     <  2:303915/148‑1 (MQ=255)
   cGTGATTCCTGCGCAAATGACCTATTAAGACTCGACAGGCAAGATCAACATTCTCGAATACCTCGCAGCTGGAATCAATTGCCGAGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtt                                                                                                                                                 <  4:224631/149‑1 (MQ=38)
    gggctGCCTGCGCAAATGACCTATGAAGACTCGACAGCCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTg                                                                                                                                                <  8:128340/145‑1 (MQ=255)
    gTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAACCGCAAACGGCGCTTTTAAGGCGCCGTTTTCGttt                                                                                                                                                <  5:169161/149‑2 (MQ=255)
      gATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTgg                                                                                                                                               >  3:61798/1‑148 (MQ=255)
      gATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTg                                                                                                                                                <  1:333776/147‑1 (MQ=255)
       aTGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCAGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTggg                                                                                                                                              >  7:172837/1‑148 (MQ=255)
         gCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGaa                                                                                                                                            >  1:267644/1‑148 (MQ=255)
          ccTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGaaaa                                                                                                                                         >  8:330225/1‑149 (MQ=255)
               gcAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTg                                                                                                                                    <  6:235365/149‑1 (MQ=255)
                cAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTgg                                                                                                                                   <  4:168604/149‑1 (MQ=255)
                 aaaTGACCTATTAATACTCGAAACGCAAGCTCAACATTCTAGTATACCGCGCACCTTGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTTAATTCTGTCCAAATCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTgg                                                                                                                                   <  1:280816/149‑1 (MQ=255)
                 aaaTGACCTATGAAGACTCGCCAGGCAAGCTCAACATTCTCGAATACCGCGCCGCTGGCATGAAGTGCCGCGGGCA‑‑AAAGTGCGCGCTGAATGCTCTCCACACGAAAACGGCGCAGTACAAGCGCCGTTTCAGTTGGGACGACGGAggg                                                                                                                                  >  2:327072/1‑149 (MQ=255)
                 aaaTGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTggg                                                                                                                                  >  2:18614/1‑149 (MQ=255)
                           tGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGc                                                                                                                          <  7:103276/147‑1 (MQ=255)
                           tGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGc                                                                                                                          <  2:267644/147‑1 (MQ=255)
                              aGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGa                                                                                                                     <  6:382796/149‑1 (MQ=255)
                                 cTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATAc                                                                                                                  >  5:137401/1‑149 (MQ=255)
                                    gACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAGGAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATAct                                                                                                                 >  1:300644/1‑149 (MQ=255)
                                       aGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGTCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTg                                                                                                                            >  6:47473/1‑133 (MQ=255)
                                       aGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGTCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTg                                                                                                                            <  5:47473/133‑1 (MQ=255)
                                            aGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCAc                                                                                                       <  3:364384/149‑1 (MQ=255)
                                               tCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACtt                                                                                                     >  7:197966/1‑148 (MQ=255)
                                               tCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACtt                                                                                                     >  1:93861/1‑148 (MQ=255)
                                                 aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCAt                                                                                                  >  7:95741/1‑149 (MQ=255)
                                                 aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCAt                                                                                                  >  8:292380/1‑149 (MQ=255)
                                                 aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCAt                                                                                                  >  5:360709/1‑149 (MQ=255)
                                                   cATTCACGTATACCGCGCACCTGTCATTAATTGCGTCGGCCA‑‑TAATTGAGCTCTGATAGGTGTCCACACGAAAACCGCGCCTTAAAAGCGCCGTTTGCGTTGGTATACCGGTGGGGCAGTGGCCTGGATACTCTTCAACGTCTTCa                                                                                                   <  8:185656/146‑1 (MQ=255)
                                                     ttCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCg                                                                                              <  7:292380/149‑1 (MQ=255)
                                                      tctcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAACACGGTGGGTCAGTGGCCTGGATActc                                                                                                                <  1:2975/131‑1 (MQ=255)
                                                      tctcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATActc                                                                                                                >  2:2975/1‑130 (MQ=255)
                                                        tcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGt                                                                                                                             >  6:362036/1‑115 (MQ=255)
                                                        tcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGt                                                                                                                             <  5:362036/115‑1 (MQ=255)
                                                        tcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTcc                                                                                           >  4:169658/1‑149 (MQ=255)
                                                          gAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGc                                                                                         >  4:352413/1‑149 (MQ=255)
                                                                 gcgcAGCTGGCATCAATTGCCGCTGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                                                  <  3:352413/149‑1 (MQ=255)
                                                                     aGCTTGCATCAATTGCCGCTGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTga                                                                                <  5:207563/147‑1 (MQ=255)
                                                                      gCTGGCATCAATTGCCGCGGACA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCTGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGAc                                                                             <  4:37313/149‑1 (MQ=255)
                                                                          gCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGa                                                                          <  6:137401/148‑1 (MQ=255)
                                                                          gCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGAt                                                                         <  4:61798/149‑1 (MQ=255)
                                                                          gCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGAt                                                                         <  6:360709/149‑1 (MQ=255)
                                                                               aaTTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGAcc                                                                     <  1:18614/148‑1 (MQ=255)
                                                                                aTTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAACCGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTg                                                                   <  2:300644/149‑1 (MQ=255)
                                                                                aTTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTg                                                                   <  8:95741/149‑1 (MQ=255)
                                                                                 ttGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGagtcgg                                                                  >  1:333897/1‑144 (MQ=255)
                                                                                     cgcgGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATg                                                              >  2:267406/1‑149 (MQ=255)
                                                                                     cgcgGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATg                                                              >  3:295861/1‑149 (MQ=255)
                                                                                          ccA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCAt                                                           >  3:145140/1‑147 (MQ=255)
                                                                                          ccA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCAt                                                           >  1:177339/1‑147 (MQ=255)
                                                                                          ccA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATga                                                         <  4:295861/149‑1 (MQ=255)
                                                                                           cA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCAt                                                           >  3:307398/1‑146 (MQ=255)
                                                                                           cA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGAACGAGGACGTGGTATGCAt                                                           >  8:124701/1‑146 (MQ=255)
                                                                                           cA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACAAATTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTAGGCaa                                                           >  3:283266/1‑145 (MQ=255)
                                                                                            a‑‑aaaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgaga                                                       <  1:267406/149‑1 (MQ=255)
                                                                                            a‑‑aaaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACCGTCGGTCAGAGGCCTGGATACTCAGCAACCACTTCATCCCGTCCCCTCTTGTTGAGACCGAAGACCTGCTATGCATga                                                         >  1:388222/1‑147 (MQ=255)
                                                                                               aaaTTGATCTCTGAATGCTGTCCAAACGTATTCGGCGCCTTAAATGCGCGGTTTGCGTTGGGAAAACGGTGGGTGAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTATTGTTGAGACCGATGAGCTGGTATGCATgagt                                                       <  2:388222/148‑2 (MQ=255)
                                                                                               aaaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCGGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgagag                                                      <  3:239401/149‑1 (MQ=255)
                                                                                                aaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATga                                                         >  1:55709/1‑145 (MQ=255)
                                                                                                aaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATga                                                         <  2:55709/145‑1 (MQ=255)
                                                                                                 aTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCGGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGac                                                    <  6:67467/149‑1 (MQ=255)
                                                                                                     aGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCa                                                <  6:193036/149‑1 (MQ=255)
                                                                                                             aTGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGGCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACtt                                         <  8:172837/148‑1 (MQ=255)
                                                                                                                cTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATcc                                             <  4:218989/141‑1 (MQ=255)
                                                                                                                cTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATcc                                             >  3:218989/1‑141 (MQ=255)
                                                                                                                     ccaaCGCAAACGGCGCTTTTAAGGCGCCGTTTTCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTg                                >  5:172236/3‑149 (MQ=39)
                                                                                                                     cAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCt                                 <  2:177339/148‑1 (MQ=255)
                                                                                                                     cAAACGAAAACGGCGCCTTAAAAGAGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTATTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTg                                >  5:145014/1‑149 (MQ=255)
                                                                                                                      aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCTACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTg                                >  3:221941/1‑148 (MQ=255)
                                                                                                                      aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCACCCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTg                                >  5:113908/1‑148 (MQ=255)
                                                                                                                      aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACATg                                >  3:108994/1‑148 (MQ=255)
                                                                                                                          cGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAAGCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACAATCCACTTCCATACCTGga                              >  7:306814/1‑147 (MQ=255)
                                                                                                                           gAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGagag                           >  8:142171/1‑149 (MQ=255)
                                                                                                                                 ggCGCTTTTAAGGCGCCGTTTTCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATc                                              >  4:371238/1‑124 (MQ=255)
                                                                                                                                 ggCGCTTTTAAGGCGCCGTTTTCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATc                                              <  3:371238/124‑1 (MQ=255)
                                                                                                                                  gcgcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCg                    <  4:307398/149‑1 (MQ=255)
                                                                                                                                  gcgcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCg                    <  4:145140/149‑1 (MQ=255)
                                                                                                                                    gcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCCCTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCgg                   >  1:262369/1‑148 (MQ=255)
                                                                                                                                     ccTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCa                 <  8:197966/149‑1 (MQ=255)
                                                                                                                                      cTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCa                 >  8:83542/1‑148 (MQ=255)
                                                                                                                                       ttAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATa               >  1:29229/1‑149 (MQ=255)
                                                                                                                                            aGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTgg          <  2:93861/149‑1 (MQ=255)
                                                                                                                                            aGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTgg          <  3:169658/149‑1 (MQ=255)
                                                                                                                                             gcgcCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTgg          >  3:352626/1‑148 (MQ=255)
                                                                                                                                              cgcCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGGGa        <  7:330225/149‑1 (MQ=255)
                                                                                                                                              cgcCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGGGa        <  7:124701/149‑1 (MQ=255)
                                                                                                                                                  gTTTGCGTTGGGAAAACGGGGGGTCAGTGGCCTGGATACTCTGCAGCCACTTCATCCAGTCCGCTCTTGTTGAGCCCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGGGACGGc    >  8:147656/1‑149 (MQ=255)
                                                                                                                                                     tGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGGGACGGCGa  <  6:145014/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGGGACGGCGA  >  NC_002947/6129039‑6129332

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: