Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,391 A→G E160E (GAA→GAG)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,3910AG87.2% 85.3 / 3.7 39E160E (GAA→GAG) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (15/19);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

CCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACAC‑TGTATCGACTGTTCGGTGC  >  NC_000913/4626147‑4626635
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
ccGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATc                                                                                                                                                                                                                                                  <  1:10587/251‑1 (MQ=255)
 cGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGCAACCGCTGTATGTCACCGGCGGAGAGTCGCTGCAACAGGGGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAgctat                                                                                                                                                                                                                                                 >  2:135832/1‑247 (MQ=255)
            agTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGcc                                                                                                                                                                                                                                      <  1:314593/251‑1 (MQ=255)
                cGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCga                                                                                                                                                                                                                                    >  1:148797/1‑249 (MQ=255)
                            aGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAAcc                                                                                                                                                                                                                        <  1:51014/249‑1 (MQ=255)
                             ggcggcggcGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGaa                                                                                                                                                                                                                     <  2:148797/251‑1 (MQ=255)
                                ggcggcgTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCt                                                                                                                                                                                                                  >  1:166659/1‑251 (MQ=255)
                                       gggtgCCTGGAAGTGCCCTTCTGATTTGCGGTAACCCTGGGGCGGGGAAATCCACTCTGGTGCTTCTAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGTGAAGAGTCGCTTCAACGGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGGTGTCGGAAACCAGCAGCGAGCATATCTGCCTGATTGTCGAAGATGAGCAACAGAAGCTGATGGTa                                                                                                                                                                                                           <  2:254462/250‑1 (MQ=255)
                                           gCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATAGCCTTGGCCTGCCGCCTGACAATCTCAATATGATGTCGGAAACCATCATCGAGCAGATCTGCCTGATTGCCGAAGAAAAGCAACCGAAGCTGCTGGTACTTg                                                                                                                                                                                                       >  1:356177/1‑251 (MQ=255)
                                               ggAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGGCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTc                                                                                                                                                                                                   <  1:286953/251‑1 (MQ=255)
                                                   gggccattgggatgggaggTAACCGAGGTTCGGGGAAATCCACTCTGCTGCTGCAAACGCTGTGCAAACTGGCCGAGCAGATGAAAACTCTGTATGTCACCTGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGAta                                                                                                                                                                                               <  2:206053/234‑2 (MQ=255)
                                                        aTTCTGATTGGCGGTAACCCTGGTGCGGGGAAAGCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCAGCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGt                                                                                                                                                                                          <  2:159612/251‑1 (MQ=255)
                                                          tCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTG‑CCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGa                                                                                                                                                                                        <  1:78805/250‑1 (MQ=255)
                                                                       aaCCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCgg                                                                                                                                                                            <  1:66711/250‑1 (MQ=255)
                                                                               ggCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGt                                                                                                                                                                   <  1:6628/250‑1 (MQ=255)
                                                                                                 gttgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCGGTATGTCTCCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTg                                                                                                                                                  <  2:75724/246‑1 (MQ=255)
                                                                                                  ttgatGCAAACGCTGTGCCAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGAGACGAGTCGCTGCAACAGGTGTCAATGCGCGCTCATCGCCTTGGCCTGCCGACGGAAAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTACCGAAGAAGAACAACCGAAGCTGATGGTAATTGACTCGCGCCAGGGGATGCAGATGGCGGTTGTATAGttcactaccagcatatgggg                                                                                                                                                >  2:61018/5‑232 (MQ=255)
                                                                                                            cGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCGTCGCCTGGCAGCGTGGCGCAGGTGCgg                                                                                                                                      >  1:283575/1‑250 (MQ=255)
                                                                                                                 tgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAAc                                                                                                                                 >  1:230308/1‑251 (MQ=255)
                                                                                                                           cccAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTAt                                                                                                                      >  1:305858/1‑251 (MQ=255)
                                                                                                                                cagATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGac                                                                                                                 <  1:117966/251‑1 (MQ=255)
                                                                                                                                       aaacgcgttATGTCGCCGGCGAAGAGTCGCTTCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATAGGTTGTCGTAAACCAGCATCGAGCAGATCAGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGTGGC‑TTATCTGACGCGCTTc                                                                                                          <  2:220261/243‑1 (MQ=255)
                                                                                                                                                  gTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATACGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGACTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATGGCCTGGCAGCGTGGAGCAGGTGCGTGAAACGGCGGC‑TTATATGACGCGATTCGCCAAAAagcg                                                                                               >  1:186084/1‑251 (MQ=255)
                                                                                                                                                              gggtttcttgtacatggggtgaTGCGCGCTCATGGCCTGTGCATGCCGACTTATAATCTCAAGAGGTGGTAGGAAACCATCATCTAACTGATCGTCCTGATGTCCGAAGAAGAGCAACCGATGGTGATGGTAATTGACGCGATCCAGGTGATGCAGAGTGCGGATGTACATGCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGTGGC‑TTATTTTGCACTCGTCGTCAATCCGCGCGTTGTGGCGaa                                                                                   <  2:197336/230‑2 (MQ=255)
                                                                                                                                                                      gCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAt                                                                             <  1:61018/249‑1 (MQ=255)
                                                                                                                                                                       acacatgTGGCAATGCGCGCACATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTAGGTGACGCGCATCGCCAAAACGCGCGGTTTGGCGATTGTCAtggt                                                                          <  2:61531/245‑1 (MQ=255)
                                                                                                                                                                              tGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATATCACTATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGCTGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAAAGTCATCGCCTGGCAGCGCGGCGACGGTGCGTGAAACGGAGGCGTTACT‑GACGCGCTTCGCCAAAAAGCGCGGTGTGGCgtctgtcctagttgtgtac                                                                   >  2:164634/1‑233 (MQ=255)
                                                                                                                                                                                         gcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCaa                                                           >  1:198132/1‑248 (MQ=255)
                                                                                                                                                                                                          gCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCtggcgggc                                        <  1:311356/250‑1 (MQ=255)
                                                                                                                                                                                                          gCCGACAGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGTACACGTAACCAAAGATGGCTCGCTGGCGGGcc                                       <  1:324760/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGc                              >  1:115983/1‑251 (MQ=255)
                                                                                                                                                                                                                       ctcAATATGTTATCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGTAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCCGGTGATGCATATGGCGGATGTTCAGGCATCGCCTGGCAGCGTGGCGCAGGGGCGGGAAACGGCGGC‑TTATCTGACGCGACTGGCCACAACGCGCGTTGTGGTGATTGTCAGGGTTGGACACGTAACCACATAAGGCTCgctggcggggccaaacgttgttgg                          >  2:318914/1‑233 (MQ=255)
                                                                                                                                                                                                                         cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTCAAGTGCTGGAac                        >  1:275427/1‑251 (MQ=255)
                                                                                                                                                                                                                         cAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGGAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGAGTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTg                            <  2:166659/247‑1 (MQ=255)
                                                                                                                                                                                                                                  gtCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACAC‑TGTATCg              >  1:263840/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGGGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGGGGCGATTGTCATGGTTGGACACGTAACCAAAGATTGGTCGCTGGCGGGCCCTAAAATGC‑GGGACACCTGTATCGAc            >  1:317503/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaCCAGCAGCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGGCTCGGTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGGCACGTAACCAAAGATGGCTCGCTGGCGGGCCCTGAAGTGCTGGAACAC‑TGTATCGACTGtt        <  1:189813/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aCCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACAC‑TGTATCGACTGTTc       <  1:135832/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATCTGACGCGCTTCGCCAAAACGCGCGGTGAGGCGGTTGTAATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAAGGCTGGAACAC‑TTTATCAACTGTTCGGTgg  >  1:334434/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGC‑TTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACAC‑TGTATCGACTGTTCGGTGC  >  NC_000913/4626147‑4626635

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: ‑