Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,927 C→A H250Q (CAC→CAA 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,9270CA95.5% 129.8 / ‑2.7 45H250Q (CAC→CAAcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (18/24);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGG  >  NC_000913/4638700‑4639162
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
gagCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCtaa                                                                                                                                                                                                                       <  1:132353/251‑1 (MQ=255)
     aGGCGGGCAGGGAGGGCATATCCCTTTACCTGAAGGGGAATGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAACGGTACCGGCGATCTTTCTGGGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCAGCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGTAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAAtc                                                                                                                                                                                                                  <  2:283632/251‑1 (MQ=255)
       gCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAAtctc                                                                                                                                                                                                                <  1:386817/251‑1 (MQ=255)
               ggCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCg                                                                                                                                                                                                        <  1:112640/251‑1 (MQ=255)
                    gCATATCCCTTTACCCGAAGGGGACTGGCGTAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTg                                                                                                                                                                                                   <  1:658643/251‑1 (MQ=255)
                              ttACCCGAAGGGGACTGGCTGAAGCATAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCTATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCa                                                                                                                                                                                          <  1:555033/250‑1 (MQ=255)
                                aCCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGCTCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATg                                                                                                                                                                                       >  1:649176/1‑251 (MQ=255)
                                  ccgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGAt                                                                                                                                                                                     <  1:503504/249‑1 (MQ=255)
                                      aGGGGACTTGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGGGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCg                                                                                                                                                                                 <  2:234905/251‑1 (MQ=255)
                                            cTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTggg                                                                                                                                                                          >  1:528808/1‑251 (MQ=255)
                                               gCGGAAGCAGAACCTTAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGaaa                                                                                                                                                                       <  1:352476/251‑1 (MQ=255)
                                               gCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGaa                                                                                                                                                                        >  1:509598/1‑250 (MQ=255)
                                                   aaGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTg                                                                                                                                                                            <  1:365952/242‑1 (MQ=255)
                                                   aaGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGATTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTg                                                                                                                                                                            >  2:365952/1‑242 (MQ=255)
                                                         aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTCCCGGCGATCTTTCTGTGGTGACTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGCGATATTTTGCCAGCCAAACGGGCAGTTAGCGAGTCAGGTTTTCCGGAGCAATGGCAAAGCAGCTGGTTTGCTAATACTCTCGGTTAACTTTTTGCTTAAGGCTATGATACCGGCTAGGGAAACTTTCCCgg                                                                                                                                                              >  2:236609/1‑251 (MQ=255)
                                                            cTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCcgcggt                                                                                                                                                          >  1:631748/1‑249 (MQ=255)
                                                                           gCCCTGAATTTAAGCGGTACCGGCCATCTTTCTGGGGTGTCTGCGGTGCGAATTGGCAAAATTACCTTAACCAGCAACTGGCCGCATCCCATTTTTTTAGGTGATTTTTTGCAAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGAAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGTAGtgac                                                                                                                                           <  2:404862/251‑4 (MQ=255)
                                                                                       aGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGGCGCccgccg                                                                                                                                >  1:397812/1‑250 (MQ=255)
                                                                                             cccGGCGATCTTTCTGTGGTGCCTCCAGGGTGTAATAGCTAAATGCCCTTAACCAGCAACTGGCCGTATCCCAGTTTTTTAGGTGATTTTTTGCCCGCCAAACGTGAAGTTAGCGAGTCAGTTTTTCAGGCTCAATGTCAAAGCCGCTGGTTTGTTAATAATCTCGGTGAGCTTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTTGCGTCGCAGTGACGACGCCAGCCGATCAAt                                                                                                                          <  2:524681/249‑1 (MQ=255)
                                                                                                  cGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGGAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCCGCCCGCCCGTCAATACCact                                                                                                                    >  2:351637/1‑249 (MQ=255)
                                                                                                        ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTa                                                                                                                  <  1:20705/247‑1 (MQ=255)
                                                                                                          gtccgatctctgccgggcgTAATAGCGAAATGACCTTAACCAGCAACCGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAAcgga                                                                                                             <  2:17474/237‑4 (MQ=255)
                                                                                                            gtggtgCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGAc                                                                                                            <  1:420908/249‑1 (MQ=255)
                                                                                                                   ctgccgggcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAAcgga                                                                                                             >  1:17474/5‑238 (MQ=255)
                                                                                                                     gccgggcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTa                                                                                                  <  1:342657/248‑1 (MQ=255)
                                                                                                                                acgtgttctcttccggtctcAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                                                       <  1:443572/232‑1 (MQ=255)
                                                                                                                                  aGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAAt                                                                                                                          >  1:503138/1‑214 (MQ=255)
                                                                                                                                      aaaTGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGAAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATAACGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGAGACGACGCCAGCCGATCAATACCCATTAACGGACCGGGCGAATAAGTACGCCATTctgctg                                                                                 >  1:622077/1‑251 (MQ=255)
                                                                                                                                           aCCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTaataat                                                                                                                                                                                                                   <  1:296238/117‑1 (MQ=255)
                                                                                                                                           aCCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTaataat                                                                                                                                                                                                                   >  2:296238/1‑117 (MQ=255)
                                                                                                                                               tagccAGAAACTGGCCTCATCCCATTTTTTTAGGTGATTTTCTACCATCCAAACGGGAAGGTATCTAGTCAGGGTGTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAAAAATCGCGGTGGGCGTTTGGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCTTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACt                                                                         <  2:649176/247‑1 (MQ=255)
                                                                                                                                                   cAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTCCCACCCAACCGGGAGGTTAGCGGGTCAGGTTTTCAGGCGCATTGCCAAAGCAGCTGGTTTGATAATAATCTCGGTGAGCGTTTTGCTTCAGGCAAGGATACCGGCTGGGAAAACTTCCCCGCCTTTAGCGCCCCAGGGCCAACGCCACCCGATCAATACCAAATAAGGGAACGGGCGAATAAGTACGACAtataacaggaagagcgccggg                                                                    >  2:443572/1‑231 (MQ=255)
                                                                                                                                                    aGCATCTGGCCGCAGCCCAGTTTTTGAGGTGATTTGTTGCGATCCAAACGGGAAGTGAGCGTGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGGTCGTTTTGCTTCTGGCAATGATACCGGCTGGGAAAACGTCCCCGGGTTTAGCGGCGCAGTGTGGACGTCAGCCGATCAACACCAATTAACGGTCCGGGCGCCTAAGGACGCCATTCTGCTGATTGCACTGACttt                                                                   <  1:415290/251‑1 (MQ=255)
                                                                                                                                                                   cccAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTg                                                     <  1:477872/250‑1 (MQ=255)
                                                                                                                                                                      aGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTTCTGATTGCCCTGACTTTTTTGGCGTTCTTTGtttt                                                 >  2:220162/1‑251 (MQ=255)
                                                                                                                                                                                  gATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCa                                      <  1:219591/250‑1 (MQ=255)
                                                                                                                                                                                           ccAGCCAAACGGGACGTTAGCGAGTCAGGTTTTCCGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGCCGCCAGCCGATCAATCCCCATTCACGGCCAGGGCGCCTAAGAACGCCATTCTGCTGATTGCACTTAATTTTATGGCGTTTTTTGttttttacaagctcacagcgacaac                            >  2:614152/1‑231 (MQ=255)
                                                                                                                                                                                             aGCCAGACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                          >  1:86822/1‑251 (MQ=255)
                                                                                                                                                                                             aGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                          >  1:203731/1‑251 (MQ=255)
                                                                                                                                                                                              gCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGttt                         <  1:351637/251‑1 (MQ=255)
                                                                                                                                                                                                     gggAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGCTCAATACCAATTAACGGACCCGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGACAAGCTCACCGCGCAACGTTTACAccc                   >  2:471254/1‑250 (MQ=255)
                                                                                                                                                                                                              gCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCCatat         >  1:400170/1‑251 (MQ=255)
                                                                                                                                                                                                                     aGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCtgg  <  1:220162/251‑1 (MQ=255)
                                                                                                                                                                                                                       gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCCATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGCCtt                                                                    >  2:108677/1‑183 (MQ=255)
                                                                                                                                                                                                                       gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACtt                                                                    <  1:108677/183‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
GAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGG  >  NC_000913/4638700‑4639162

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: