Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,129,160 A→C 5.1% intergenic (+22/+45) 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,1600AC 5.1% 234.0 / 3.5 79intergenic (+22/+45)PP_5376/PP_5377hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base C (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

TGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAG  >  NC_002947/6129015‑6129307
                                                                                                                                                    |                                                                                                                                                    
tgtgACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGaaaa                                                                                                                                                   >  3:37313/1‑149 (MQ=255)
     ccTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGgcg                                                                                                                                              >  4:364384/1‑149 (MQ=255)
       tCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGCAAACGgcgct                                                                                                                                            >  6:169161/1‑148 (MQ=255)
       tCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGcc                                                                                                                                            >  5:193036/1‑149 (MQ=255)
            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                                                                         >  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                                  >  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)
                                                                                                                                                    |                                                                                                                                                    
TGTGACCTCGACCCGGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCA‑‑AAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAG  >  NC_002947/6129015‑6129307

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: