Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,176 A→G 6.2% intergenic (+38/+29) 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,1760AG 6.2% 242.4 / 4.4 82intergenic (+38/+29)PP_5376/PP_5377hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base G (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGG  >  NC_002947/6129028‑6129324
                                                                                                                                                       |                                                                                                                                                     
cGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTaaaa                                                                                                                                                       <  2:293422/149‑1 (MQ=255)
        cTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGtt                                                                                                                                               >  5:67467/1‑149 (MQ=255)
           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                                                                                                                    <  5:47473/133‑1 (MQ=255)
                                                  aGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGTCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTg                                                                                                                    >  6:47473/1‑133 (MQ=255)
                                                  aGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGttt                                                                                                                                              >  2:361238/1‑108 (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                                                                                          >  5:360709/1‑149 (MQ=255)
                                                            aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCAt                                                                                          >  7:95741/1‑149 (MQ=255)
                                                            aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCAt                                                                                          >  8:292380/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                                      <  3:371238/124‑1 (MQ=255)
                                                                                                                                            ggCGCTTTTAAGGCGCCGTTTTCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATc                                      >  4:371238/1‑124 (MQ=255)
                                                                                                                                             gcgcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCg            <  4:145140/149‑1 (MQ=255)
                                                                                                                                             gcgcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCg            <  4:307398/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  <  3:169658/149‑1 (MQ=255)
                                                                                                                                                       aGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTgg  <  2:93861/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATACCTGG  >  NC_002947/6129028‑6129324

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: