Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,173 A→T 6.1% intergenic (+35/+32) 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,1730AT 6.1% 225.0 / 4.7 82intergenic (+35/+32)PP_5376/PP_5377hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base T (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

CCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATA  >  NC_002947/6129027‑6129319
                                                                                                                                                     |                                                                                                                                                   
ccGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTaa                                                                                                                                                    >  8:323099/1‑148 (MQ=255)
 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                                                                                                             <  2:267644/147‑1 (MQ=255)
                                       tGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGc                                                                                                             <  7:103276/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                                                                                     >  8:292380/1‑149 (MQ=255)
                                                             aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCAt                                                                                     >  7:95741/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                                                            <  6:360709/149‑1 (MQ=255)
                                                                                      gCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGAt                                                            <  4:61798/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                                              >  1:177339/1‑147 (MQ=255)
                                                                                                      ccA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCAt                                              >  3:145140/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                                            <  2:55709/145‑1 (MQ=255)
                                                                                                            aaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATga                                            >  1:55709/1‑145 (MQ=255)
                                                                                                             aTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCGGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGac                                       <  6:67467/149‑1 (MQ=255)
                                                                                                                 aGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCa                                   <  6:193036/149‑1 (MQ=255)
                                                                                                                         aTGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGGCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACtt                            <  8:172837/148‑1 (MQ=255)
                                                                                                                            cTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATcc                                >  3:218989/1‑141 (MQ=255)
                                                                                                                            cTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATcc                                <  4:218989/141‑1 (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)
                                                                                                                                                     |                                                                                                                                                   
CCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATA  >  NC_002947/6129027‑6129319

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: