Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,567 A→G E130E (GAA→GAG creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,5670AG96.4% 159.8 / ‑2.1 56E130E (GAA→GAGcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (27/27);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-01

GCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGT  >  NC_000913/4638333‑4638812
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
gCAAAAACTAGTTAGGCCGCTTATCGCTAGTCCTGTGAACGCGCTTTATCCGGTGCAGTAAGAGGATAACACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGCTGATGGCAATCAGAACGTGGAGGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTTACGCTCAAAGCCGTGTACGATGTTTCGCGTCTTAGCGAGATCAACGCGCCacat                                                                                                                                                                                                                                        >  2:242792/1‑247 (MQ=255)
 cAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAatat                                                                                                                                                                                                                                      >  1:56709/1‑251 (MQ=255)
  aaaaaCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATc                                                                                                                                                                                                                                     >  1:152849/1‑251 (MQ=255)
   aaaaCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATc                                                                                                                                                                                                                                     >  1:19010/1‑250 (MQ=255)
     aaCTAGTTGGGCCGCTTATCGCTATTCCTGTGACAGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCACTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCAACGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCGCCGCGCCAAATATCatc                                                                                                                                                                                                                                  >  2:266253/1‑249 (MQ=255)
       cTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCtt                                                                                                                                                                                                                                >  1:233656/1‑251 (MQ=255)
         cGTTGGGCAGCTTATCGCTATTCCTGTGTCCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAg                                                                                                                                                                                                                              <  1:307783/251‑1 (MQ=255)
             gggCCGCTTATCGCTCTTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGCAGTTTTATCCATTTTTGGTTACCTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCa                                                                                                                                                                                                                           >  1:284805/1‑250 (MQ=255)
                            aTTCCTGTGACCTAGCCTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATtgtg                                                                                                                                                                                                             <  1:26692/249‑1 (MQ=255)
                            aTTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATtgtg                                                                                                                                                                                                             <  2:152849/249‑1 (MQ=255)
                                cTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTCAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGc                                                                                                                                                                                                       >  1:16417/1‑251 (MQ=255)
                                  gtgACCGAGCTTTATTCGGTGCAGGAAGAGGATAAAACCTTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTGGATGGCAATCAGAACGTGGAAGAACGCAAGATATGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCTCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCg                                                                                                                                                                                                      <  2:486487/250‑1 (MQ=255)
                                          gCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCAAAATATCATCTTAGGCAAGCCCTTTATTGTGCTTAGCGTCGGGGATg                                                                                                                                                                                             >  2:418483/1‑251 (MQ=255)
                                               ataCGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGACGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGcgcgg                                                                                                                                                                                        >  1:149178/1‑250 (MQ=255)
                                                        aGGAAGAGGATAAAACCGTGGAGCGGATACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATtggt                                                                                                                                                                               <  2:348797/251‑1 (MQ=255)
                                                                gATAAAACCGGGGAGCGGAAACGAAGTTTTATCCATCTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACATTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCTGTTATTGTGATGAGCGTCGGGGATGCGCGTTGTATTGGTGTGGTGaa                                                                                                                                                                       <  2:312091/251‑1 (MQ=255)
                                                                                ggAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaac                                                                                                                                                       >  1:396787/1‑251 (MQ=255)
                                                                                  aaaCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAAGCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaacgg                                                                                                                                                     <  2:19010/251‑1 (MQ=255)
                                                                                      gAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaacggaacg                                                                                                                                                 >  1:177072/1‑251 (MQ=255)
                                                                                        aGTTTTATCCAGTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGGGATTAGCGTCGGGGATGCGCTTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGgc                                                                                                                                               <  2:173303/251‑1 (MQ=255)
                                                                                              aTCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGTTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGGGTGTTGAAAGCGCCTGAAGTTAACGGAACGTCGCTGAc                                                                                                                                         <  1:331233/251‑1 (MQ=255)
                                                                                              aTCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAc                                                                                                                                         >  1:222853/1‑251 (MQ=255)
                                                                                                 catgtttGGTTACCTGAGTCACTGATGGTTGATGGCAATCAGAACGTGGAAGAAAGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAt                                                                                                                                      <  2:222853/247‑1 (MQ=255)
                                                                                                        ggTTACCTGAGTCATTGATGGTTGATGGCAAGCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTCAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAccc                                                                                                                               <  2:16417/251‑1 (MQ=255)
                                                                                                           tACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGcc                                                                                                                                                                                                                       <  1:193986/160‑1 (MQ=255)
                                                                                                           tACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGcc                                                                                                                                                                                                                       >  2:193986/1‑160 (MQ=255)
                                                                                                                    cATTGATGGTTGAGGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTGTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGtt                                                                                                                    <  2:56709/250‑1 (MQ=255)
                                                                                                                    cATTGATGGGTGCTGGCAAACAGAACGTGGAAGAAAGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTCACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCACCTATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGACGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGGAGTTAACGGAACGGCGCTGCCCATTGAACCCGGATCCGGGtt                                                                                                                    >  2:251357/1‑250 (MQ=255)
                                                                                                                      ttgatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGtt                                                                                                                    >  1:105091/1‑248 (MQ=255)
                                                                                                                           ggttgatgtCAATAAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGAGTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAg                                                                                                            <  2:316273/251‑1 (MQ=255)
                                                                                                                              tgatggCAATCAGAACGTGGAAGAACGCAAGATATGGGTTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGCGGGGAAAGCGCCTGAAGTTAACGGACCGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCCAGGCg                                                                                                         >  2:24101/1‑251 (MQ=255)
                                                                                                                                atggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCggg                                                                                                       >  1:45351/1‑251 (MQ=255)
                                                                                                                                        taGAACGTGGCCGAACGTAAGATAGGGATTTAACCCGGTCTGGTCTGTCACAGTGATTTAACCATAAAAGCCGATTTCGATGTTTCGCGTCTTAGCTAGTTCAACGCGCCAAATATCCTCTTAGGCAAGCCATTTATTGTGAATAGCGCCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGTTGCTGACCATTGAACCCGGCCCCGGGTTAGAGCACGGCGGGccggtct                                                                                                <  2:130830/249‑7 (MQ=255)
                                                                                                                                                ggAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATcc                                                                                       >  1:445126/1‑251 (MQ=255)
                                                                                                                                                  aagaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGGGGGCATATcccc                                                                                     >  2:324496/1‑250 (MQ=255)
                                                                                                                                                     aaCGCAAGATAGGGATTTATACCGGGCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTGCGCGTCGTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAGTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTa                                                                                  <  2:488736/251‑1 (MQ=255)
                                                                                                                                                      aCGCAAGATAGGGATTTATTCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCTCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTGTAc                                                                                 <  2:300284/251‑1 (MQ=255)
                                                                                                                                                                      tataCCGGTTGGGTCTGGCACAGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGTGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGc                                                                  <  2:345690/250‑1 (MQ=255)
                                                                                                                                                                                 ggTCTGTCACAGTGTTTTAACCCTAAAAGCCGATTTCGTTGTTTCTCGTCTTATCGAGCTCAACGCGCCAAATATCATCTTTGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCACTTTACCCGAAGGAGACAGGCGGACGCAGaccc                                                      <  2:283673/251‑4 (MQ=255)
                                                                                                                                                                                 ggTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAcc                                                      <  2:165483/251‑1 (MQ=255)
                                                                                                                                                                                   tCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGAAGGGGGAGCATATCCCTTTACCCGAAAGGGACTTGCGGagagagaacc                                                      >  2:83155/1‑240 (MQ=255)
                                                                                                                                                                                        cacaGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCt                                               <  1:187730/251‑1 (MQ=255)
                                                                                                                                                                                                         aaaaGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGGACCTGAAGGTGAATATGGCCCTGAAtt                              >  1:82723/1‑251 (MQ=255)
                                                                                                                                                                                                           aaGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATAGCCCTTTACCCGAAGGGGACTGGCGGGAGCAGGACCTGAAGCttaaaattgaccttcatttt                            >  2:4354/1‑232 (MQ=255)
                                                                                                                                                                                                            aGCCGATTTCTATGTTTCGCGTCTTAGTGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTaa                           <  2:342239/251‑1 (MQ=255)
                                                                                                                                                                                                            aGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTaa                           <  2:233656/251‑1 (MQ=255)
                                                                                                                                                                                                               cGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACTCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCGTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCg                        <  1:390290/251‑1 (MQ=255)
                                                                                                                                                                                                               cGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGGGACTGGCGGAGGGAGAACCTGGAGATGAATATGGCCCTTGATTTAAggg                        >  1:441729/1‑249 (MQ=255)
                                                                                                                                                                                                                     cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAc                                                                                                                                         <  1:435075/132‑1 (MQ=255)
                                                                                                                                                                                                                     cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAc                                                                                                                                         >  2:435075/1‑132 (MQ=255)
                                                                                                                                                                                                                     cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAcc                   <  2:396787/250‑1 (MQ=255)
                                                                                                                                                                                                                     cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCg                  <  1:9801/251‑1 (MQ=255)
                                                                                                                                                                                                                          tttagcgTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAATGGTGTGGTGAAAGCGCCTGAAGTTAAAGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCATGCATATCCCTTTACCCGCAGGGGACTGGCGGGAGCAGAACCTGAAGCTGAAAATGGCACTGACttttaacggaacagtcgat             >  1:382547/5‑235 (MQ=255)
                                                                                                                                                                                                                          tttCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGCAGCTGAATATGGCCCTGAGTTTAAGCGGTACCGGCGAt             >  1:382535/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCt       <  2:82723/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGGCTGGCGGAAGGAGAACCTGAAGCAGAAAATGGCCCTGAAATTAAGCGGGAACGGCGaattttttggggg  >  2:319578/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
GCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGT  >  NC_000913/4638333‑4638812

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: