Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,600 A→C P141P (CCA→CCC 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,6000AC92.5% 160.9 / ‑0.7 54P141P (CCA→CCCcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base C (20/29);  total (22/32)
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.00e+00

ATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACC  >  NC_000913/4638361‑4638846
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
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)
                                                                                        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)
                                                                                                                                                                                                         agggagcTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTAGTGTGATTAGCGACGGGGAAGCGCGTGGTATTGGTGTGGTGAAAGCGCATGAAGTTAACGGAACGGCGCTGACCATTGAGCCCGGCACCGGGTTAGAGCAAGGCGGGAAGGGCGTGCATATCCCTTTGCCCGAAGGGCACTGGCGGAAGCAGAACCTTAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTTTgtgg                                     <  2:250719/244‑1 (MQ=255)
                                                                                                                                                                                                         aGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGGCTGGCGGAAGGAGAACCTGAAGCAGAAAATGGCCCTGAAATTAAGCGGGAACGGCGaattttttggggg                                    >  2:319578/1‑239 (MQ=255)
                                                                                                                                                                                                            gagcTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTTTGTGGTGcc                                 >  1:211676/4‑251 (MQ=255)
                                                                                                                                                                                                              gcTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTg                               <  1:406591/250‑1 (MQ=255)
                                                                                                                                                                                                                                     tAGGCACGCCCTTTATTGTGATTTGTGACGGGGATGCGCGTTGGATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGTGCTGGCCATTGAACCCGGCACCGGGTTAGAGCTAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGTATAACCTGAAGTTGAATATGGCCCTGAATTTAATCGGTACCGTCGATCTTTCTATGGTGCCTGCCGGGCGTAATAGCGAAAt            <  2:434649/247‑1 (MQ=255)
                                                                                                                                                                                                                                      aGGCAAGCCCTTTATTGGGATTAGCGTCGGGGAGGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCTTTGAACCCGGCACCGGGTTAGAGAAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCATGCCGGGCGTAATAGCGAAATGCCCt       <  2:357375/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCtt      <  2:177072/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGAGGTATTCGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGGTAGCGCAAGGCGGGCCGGGCGTGCATATCACTTTACCCGAAGGGGACTGGCGGAGCAGAACCTGAAGCTGAATATTGCCCTGTCTTTTAGCGGTACCTGGGATCATTATGTGGTGACTGCCGGGCGTAATAGCGaaaagccctt      >  1:138318/1‑243 (MQ=255)
                                                                                                                                                                                                                                           aGCCCTTTATTGGGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAAcc  <  2:82346/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
ATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACC  >  NC_000913/4638361‑4638846

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