Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,320,543 G→A L41F (CTT→TTT)  phnI ← carbon‑phosphorus lyase core complex subunit PhnI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,320,5430GA96.8% 176.1 / ‑2.4 63L41F (CTT→TTT) phnIcarbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (32/28);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

AGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGG  >  NC_000913/4320300‑4320788
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
aGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAActgctg                                                                                                                                                                                                                                                 <  1:39385/250‑1 (MQ=255)
             cagcTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGCGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCg                                                                                                                                                                                                                                    >  1:184130/1‑250 (MQ=255)
             cagcTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCgg                                                                                                                                                                                                                                   >  1:166175/1‑250 (MQ=255)
                cTGGACGCAGGGAATGTCTTTATAAACGACGGCGATACGGCGTTCAAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCACCGTGATGCGGGAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGATTCAGCACCAGCGCCGCCCGTTAGATGGCGGCAATGCCTACTTCGGTCATCCGGCGATCTACCGCGAGGTTAAACAGCTGGTCAATCTGGGGGACGc                                                                                                                                                                                                                                >  2:200910/1‑251 (MQ=255)
                     cgccgGGGATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTTCTGTTCAATCTGCGCGACGCTCAGt                                                                                                                                                                                                                           >  1:347503/1‑251 (MQ=255)
                            atgggcTTTATAAACGCCGGAGATTCGGCGTTCGGGCCGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCa                                                                                                                                                                                                                    <  1:199623/246‑1 (MQ=255)
                                tCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATc                                                                                                                                                                                                                >  1:341021/1‑251 (MQ=255)
                                                  aTACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCg                                                                                                                                                                                              >  1:52666/1‑251 (MQ=255)
                                                        cGTTCGAGCCGCATCCCGGTGGTGTCGAGCGGCCCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGctctc                                                                                                                                                                                        <  1:4950/251‑1 (MQ=255)
                                                            cGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTAAATCTGGGCGACGCTCAGTTCGGGCAAATCGGAATCGCCTCGGCGTCGGCTCTCTTGc                                                                                                                                                                                    >  1:155186/1‑251 (MQ=255)
                                                                 cGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCCAATCGGTATCGCCTCGGCGTc                                                                                                                                                                                               >  1:31604/1‑235 (MQ=255)
                                                                  gCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCg                                                                                                                                                                              <  1:228159/251‑1 (MQ=255)
                                                                   cATCCCGGTGTTGTCAAGCGCCTTCCTTCCCGCCAGCTTCGCCACCGTGAGGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGt                                                                                                                                                                             <  1:140372/251‑1 (MQ=255)
                                                                      ccggggggggtggaggggcctgctcaCCGCCAGGTCCGCCAACTTGCTGCGGTAGGCGCGAGGTGGGAAAACCGCTTCAGTGGTATGGCCGCTGGCCTGTTTCGGCGCCGGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCAGCCCTCGATCTACCGCGAGGTTAAACTGCTGTTCAAGCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTggg                                                                                                                                                                          <  2:143917/226‑1 (MQ=255)
                                                                        cGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTggg                                                                                                                                                                          >  1:255174/1‑249 (MQ=255)
                                                                           gggggtggcgcgcaccGCTTACACCCAGCGTCGCCAGCGTGATGCGGTGGGCGCGCAGCAGGAAAATCGCTTAAACGTTACCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGCCGCTCAGTTCGGGCAAAGCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTTGgcggcg                                                                                                                                                                     <  1:206502/236‑1 (MQ=255)
                                                                                 gggtcGGCTCGCTCACCGCCAGTTTCGCCAACGTGAGGCGGTGGGCGCGAGGAAGGAAATTCCATGTATCGTAACCGCCGCTGGCCTGTTTCAGCGCCATCGCCGCCAGTTCGCGGTCGGCAATTCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGTTATCGCCTCGTCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATc                                                                                                                                                               <  2:371851/247‑1 (MQ=255)
                                                                                 cGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCGGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGAt                                                                                                                                                                <  1:234500/250‑1 (MQ=255)
                                                                                          cGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCg                                                                                                                                                      >  1:354313/1‑251 (MQ=255)
                                                                                            cTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGcc                                                                                                                                                    >  1:123253/1‑251 (MQ=255)
                                                                                              tACCGCCAGCTTCTCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGccc                                                                                                                                                   <  1:64164/250‑1 (MQ=255)
                                                                                              tACCGCCAGCTGCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCACCGTTATCGCCGCTGGCCGGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGccc                                                                                                                                                   <  1:64156/250‑1 (MQ=255)
                                                                                                ccGCCAGCTTCGCCAACGTGATACGGTAGGCGAGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGGTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTCTCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCt                                                                                                                                                >  2:252229/1‑251 (MQ=255)
                                                                                                   cccgattCGCCAACGTGATGCGGTAGACGAGCAGCCGGAAAATCGATTCAACGGTATAGCCGCTGGCCTGTTTCGGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACAGCGAGGTAAAACGGCTGTTAAATCTGGGCGACGATCAGTTCGGGAAAATCGGTATCGCCTCAGCGGCGGCTCTCTTGCAGGGCGTGGGCGGCGGCGATCGCCTTCTCGCCCCCtgtc                                                                                                                                             >  2:354971/6‑248 (MQ=255)
                                                                                                       cTTCGCCAACGAGATGCGGGAGGCGCACAGCAGGAAAATCGCTTCAACGATATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTAGGCACTGCCGCCTTCGGTCATCAAGCGATCTACCGCGAGGTAAAACTGCTGGTCAAGCTGGCCGACTCCAAGTTCGGGCAACTCGGAAACGCATCGGCGGCGGCTATCTTGCCGGGCGTGGGAGGCGGCGGTCGCATTCTCGCTCCCTTACACgg                                                                                                                                         >  2:304347/1‑251 (MQ=255)
                                                                                                                       cGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTGCCGCCTCTGGTCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTCCCGCGAGGTTAAACTGCTGGTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGcac                                                                                                                         <  2:341021/251‑1 (MQ=255)
                                                                                                                       cGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGcac                                                                                                                         <  1:347766/251‑1 (MQ=255)
                                                                                                                        ggTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGcaca                                                                                                                        <  1:4460/251‑1 (MQ=255)
                                                                                                                          tAGGCGCGCCGCCAGAAAATCGCTTCAACGATCTCGCCGCTGTCCTGTTTCAGCGCCAGCGCCGCCAGATAGCGGTCGGCAATGCCGCCTTCGGTCATAACGAGATCTACCGCGAGGGTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCCACTCGGGATCGCCTCGGCGTCGGCTCTCTAGCAAGACGTGGGCGGAGTCGATCGCCATCTCGCCCCCTTTCACGGCAACGTACATCAGCACAcc                                                                                                                      >  2:112424/1‑251 (MQ=255)
                                                                                                                          tAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACAcc                                                                                                                      >  1:108269/1‑251 (MQ=255)
                                                                                                                              cgcgcAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTcc                                                                                                                   >  1:213520/1‑250 (MQ=255)
                                                                                                                                     cagGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGAGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAgg                                                                                                           >  1:28176/1‑250 (MQ=255)
                                                                                                                                                 ttCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAAt                                                                                               >  1:255197/1‑251 (MQ=255)
                                                                                                                                                   cAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATa                                                                                              <  1:164122/250‑1 (MQ=255)
                                                                                                                                                      cGTTATCGGCGCCGGCCTGTTTCAGCGCCGGCGACGCCAGGTCTCGGTCGGCAATGCCGCATTCGGGCATCACGCGATCTACAGCGAGGTTAAAACGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCATCGACTCTCTGGCAGGGCGTGGGCGGCGGCGATCGCCTTGTCGCACCCTTTCACGGCAACGTACATCAGCACAGCTCCACAGGAGGGGTTCGCGGAAtcacc                                                                                           >  2:98465/1‑246 (MQ=255)
                                                                                                                                                       gTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGcc                                                                                           >  1:252802/1‑249 (MQ=255)
                                                                                                                                                        gttccgcgccgtgCCTGTTTGGGCGGCAGCGCCGACCGTCCGCAGTCGGCCATGCAGCCTTCGGTAATCGCTCGATCTACCGCGAGGTTAAACTGCTGTTCAAGCTGGGCGACGCACGGTGCTGGCAATTCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGTGCGGCGTCGGTCGCCTTCTCGGCCCCTTTCACCGCAACGTACATCAGCACACCTCCACAGGAGTGGTTCGCGGAAGAGCca                                                                                          <  2:97368/237‑1 (MQ=255)
                                                                                                                                                                 cTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTAATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGCGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTAACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGAGGAATCGCCAGCCGGCGCTcc                                                                               >  1:238417/1‑250 (MQ=255)
                                                                                                                                                                     ccTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCcac                                                                           >  1:101611/1‑251 (MQ=255)
                                                                                                                                                                         tttCAGAGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGCGGATAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAACTCGGTATCGACTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGCTCGCGGACTAGCCAGCAGGCGCTcccccccggg                                                                       >  2:165665/1‑242 (MQ=255)
                                                                                                                                                                            cagcgccagcgcCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGTGCGGCGTCTATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACGCCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCa                                                                     <  1:349817/250‑1 (MQ=255)
                                                                                                                                                                               cgacagcgccGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGAt                                                                 <  1:258495/248‑1 (MQ=255)
                                                                                                                                                                                gccagcgacGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTAAGTTCGGGAAACTCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTATCGCCCCCTTTCACGGCAACGTAAATCAGAACACCTCCACATGAGTGGTTCGCGGAATAGAAAGCAGGAGCTCGCCACAGGGCAGGaac                                                                >  2:368569/1‑249 (MQ=255)
                                                                                                                                                                                    gcgccGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGt                                                            <  1:307309/251‑1 (MQ=255)
                                                                                                                                                                                    gcgccGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGt                                                            <  1:177052/251‑1 (MQ=255)
                                                                                                                                                                                     cgccgcCAGTTCGCTGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCGCGGCAACGTACATCAGCACACCTTCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTc                                                           <  1:171239/251‑1 (MQ=255)
                                                                                                                                                                                       acgccAGTGCGGGGTCGGCTATGCCCCCATGGGTCCGCACGCGAGCTCCCGCGGGGGGAAACTGCTGTGCAATCTGGGCGACGCTCAGGTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGTGGCGTCGAGCGCCTTCTCGCCCCCTTTCACGTCAACGGACATCAGCACGCCTCCACAGGAGTGGTTCGCGGAATAGCCAGCAGGCTCTCGCTACAGGGCAGGATCAGGTc                                                           <  2:173611/248‑1 (MQ=255)
                                                                                                                                                                                                gcgggTCTGCAATGTCGGCTTCGGTCCTCCCGCGATCTTCCGCGTGTTTAAACTGCTGTTCAATCTGGGCGACGCTCTGTGCGGGCAAAGCGGTATCGCCTCGGCGTTGGCGCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCGTCGCGCCCCCTGTCACGGCAACGTACATCAGCACTCCTCCACATGAGTGGTTCGCAGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTAGATGCCGAg                                                  <  2:7836/246‑1 (MQ=255)
                                                                                                                                                                                                       ggCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACgg                                         >  1:327363/1‑251 (MQ=255)
                                                                                                                                                                                                                  cTTCGGGCATCACCCGATCTACCGCGAGGTTAAACTGCTGTTCACTCTGGGCGCCGCCCTGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCGCTTGAAGGGCGTGGGCGTGTTCGATCGCTTACGCGCCCCCTTTCACGGCAACGTACATCAGCACCCCTCCACATGCGTGGTTCTCGTAAGAGCCAGCCGGCGCTCCCCCCAGGTCAGGAACAGGTCGATGCCGAGCgggtgaggatgcggtcgcc                              >  2:18106/1‑235 (MQ=255)
                                                                                                                                                                                                                     cGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTGCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCgg                           <  1:284270/251‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCg                            <  1:42838/250‑1 (MQ=255)
                                                                                                                                                                                                                          atcaCGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGt                      <  1:101170/251‑1 (MQ=255)
                                                                                                                                                                                                                            caCGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGtt                     <  1:242136/250‑1 (MQ=255)
                                                                                                                                                                                                                               gcgATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCGCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGACCGGCTGCGGGCGCTCGGTGGGTTCGTg                 >  2:153954/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gATCTGCCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTAGCCACAGGTCAGGATCAGGTCGAGGCCGAGCGGGCACGGGTGCGGGCGCTCGGTGAGTTAGTGc                >  2:364307/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCGGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGc                >  1:291188/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTAAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACAACTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGCCAGGATCAGGTCGATGCCGAGCGGGGACGGATGCGGGCGCTCGGGGAGTTCgagc                >  2:327035/1‑247 (MQ=255)
                                                                                                                                                                                                                                   tCTACCGCGGGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGAATCGCCTCGGCGTCGTCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCCGCACACCTCCACATGAGTGGGTCGCGGAATAGCCACCAGGCGCTCGCCACAGGTACGGATCAGGTCGATGCCGAACGGGACAGGAAGCGTGCGCCCGGGGAGTTCgcgc                >  2:311452/1‑245 (MQ=255)
                                                                                                                                                                                                                                    cgaCCGCAAGGCTAACCTGCGGTTCAATCTGGCTGATCCGCAGTTCGGGCAAGCGGGTGTCGCCTCGGCTTCGGCTCCCTTGCAGTGCGTGGGTGGCTTCCATCGCCTTCTCGCCCCCGTTCACGGCAACGTACATCAGCACACCTCCTCATGAGTGGTTCGCGGAAGAGCCAGCAGGCGATCGCCACAGGTCAGGATCAGGTCGATGCCGAGAGGGAACGGATGCGGGCGCCCGGTGAGTTCGTGCAGaa            <  2:213520/249‑1 (MQ=255)
                                                                                                                                                                                                                                      gccGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAg              <  1:35193/246‑1 (MQ=255)
                                                                                                                                                                                                                                           gAGGTTAAACTGCTGTTCAAGCTGGGCGACGCTCAGTTCGGGCAAATCGGGATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTGGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTc     <  1:263226/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCgg  <  1:280614/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
AGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGG  >  NC_000913/4320300‑4320788

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