Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,735 T→C P186P (CCT→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,7350TC84.1% 101.4 / 6.2 44P186P (CCT→CCCcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (18/19);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

AGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGG  >  NC_000913/4638492‑4638977
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAgggg                                                                                                                                                                                                                                               <  1:195437/251‑1 (MQ=255)
                     cTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTg                                                                                                                                                                                                                          >  1:258279/1‑250 (MQ=255)
                                aTTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCCGAACCTGAAGCTGAATATg                                                                                                                                                                                                              >  1:210066/1‑251 (MQ=255)
                                    aaCGCTAAAAGCCGATTTCTATGTGTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCATGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAatattgcac                                                                                                                                                                                                          <  2:58784/251‑6 (MQ=255)
                                      cactAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGTAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGACAGGGGGTGCATATCCCTTTACCCGAAGGAGAACGGCGGAAACAGaaacagacaattgaaatgctccc                                                                                                                                                                                                          >  1:219066/4‑229 (MQ=255)
                                        ctAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGCAGTTAACGGAACGGCGCTGACCCTTGACCCCGGCACCGGGTTAGAGCAAGGCGGGAAGGGAGAGCATATCCCTTTACCCGAAGGGGACCGGCGGGAGCAGGATCTGGAGCTggaaatgggcatgac                                                                                                                                                                                                      >  2:22136/2‑237 (MQ=255)
                                        ctAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGGACCCAAAGCTGAATATGGCCCTGa                                                                                                                                                                                                       >  1:74445/2‑250 (MQ=255)
                                                  aTTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAg                                                                                                                                                                                                                       <  1:166537/224‑1 (MQ=255)
                                                    ttctaTGTTTCGCGTCTTAGCTATCCCATCGCGTCAAATATCAGCTTAGGCTTTTCCTTTTGTGTGATTAGCGTCGGGGATGCGCGTGGTATTGTTGTGGTGAAAGCGCCCGAAGTTAACGGAACGGCGCTTACCATTGAACTCGGCGCCGGGTTAGAGCAAGGCTGGCTGGGCGTTCATATCTCTTTATCCTAAGGGGACTGGCGTAAGCAGATCTTGAAGCGGAATATGGCTCTGAATTTaggtgata                                                                                                                                                                                           <  2:346976/246‑8 (MQ=255)
                                                    ttCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGGAGCTGGATAAGGCCCTGAATTTAAgggggac                                                                                                                                                                                          >  1:292915/1‑245 (MQ=255)
                                                     tCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAGCgggaac                                                                                                                                                                                         >  1:318478/1‑247 (MQ=255)
                                                          gTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCg                                                                                                                                                                                     <  1:259013/250‑1 (MQ=255)
                                                                cttCTTAGCGAGCACAACGCGCCAAATATCATCTTAGGCAAGCTCTTTATTGTGATTAGCTTCGGGGATGCTCGAGGTATTGGCGTGTTGAAAGCGCCGGAAGTTAACGGATCGGCGCTGTCCATTGAACCCGGCTCCTGTTTAGAGCAAGGCGGGCAGGGCGTGCATAGCCCTTTCCCCGGAGTGGACTGGCGGAAGCAGTACCTTAAGCTGAATATGGCCCTGAATTTAAGCGGTACCTGTGATCtt                                                                                                                                                                                <  2:228434/247‑1 (MQ=255)
                                                                 gTCTCAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCTGAAGGGGACTGGCGGGAGCAGGACCATGAGCTGGATATGGCCCTGGATTTTAGCGGGACCAGCGGTCttttt                                                                                                                                                                             >  2:189292/1‑249 (MQ=255)
                                                                   cTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCCGAATTTTAGCGGTACCGGCGATCTTTCtgt                                                                                                                                                                           >  1:179235/1‑251 (MQ=255)
                                                                        cgagcTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGc                                                                                                                                                                      <  1:384838/247‑1 (MQ=255)
                                                                             acAACGCGCCAAAGATCACCTTAGGCAAGCCACTTAGTGTGATTAGCGTCGGGGATGCGCGTGGGAGTGGTGTGGTGAAAGCGCCGGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTggc                                                                                                                                                                 <  2:210066/250‑1 (MQ=255)
                                                                                           atCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGGGTAATAGCgg                                                                                                                                                   >  1:142481/1‑250 (MQ=255)
                                                                                             catcTTAGGCAAGCCCTTTATTGAGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAGGGGGACTGGCGGGAGCAGAACCTGAAGCTGGATATGGCCCTGAATTTAAGCGGAACCGGCGATCATTCTTGGGTGCCTGCCGGGCGGAATAGctac                                                                                                                                                  >  1:214374/4‑247 (MQ=255)
                                                                                               tcTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGGAATg                                                                                                                                               >  1:333462/2‑251 (MQ=255)
                                                                                                  tAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAAt                                                                                                                                                <  1:285870/247‑1 (MQ=255)
                                                                                                  tAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGa                                                                                                                                              <  1:189292/249‑1 (MQ=255)
                                                                                                   aGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAACGGGACCGGCGATATTTCTGTGGTGCCTTCCGGGCGTAATAGCGAAATTACCt                                                                                                                                           >  1:283933/1‑251 (MQ=255)
                                                                                                         gccctttATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGGAGGGGCCTGGCGGAAGCGGGACCTGGAGATGAATTAGGCCCCGGAGTTATACGGGACCGGCGGTCTTTTTGTGGTGCCTGCCGTGGATAATAGCGagagacacttaacc                                                                                                                                      >  2:59510/5‑237 (MQ=255)
                                                                                                                   gtgATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGACTATGGCCCTGAATTTAAGCGGTACCGGCGATATTTTTGTGGTGCCTGCCGGGCGTAATAGCGGAATGACCTTAACCAGCAACTGGcc                                                                                                                           >  1:241807/1‑251 (MQ=255)
                                                                                                                           gggcGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTTAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCa                                                                                                                   <  2:333462/248‑1 (MQ=255)
                                                                                                                           cGTCGGGGTTGAGCGTGGTATTGGTGTTGTGAAAGCGCCGTAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCCGAAGGGGACTGGCGGAATCAGCACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTTCCGGCGATCTTTCTGTGGTGCCTGCCGGTCGTAATAGCGAAATGACCTTAACCAGCAACAGGCCGCATCCCa                                                                                                                   <  2:291218/251‑1 (MQ=255)
                                                                                                                            gTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGGAATGACCTTAACCAGCAACTGGCCGCCTccc                                                                                                                    >  1:118128/1‑249 (MQ=255)
                                                                                                                             gcGGGGATGCGCGGGGTATTGGTGTGGTGAAAGCGCCTGAATGTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAATCATAACCTGAAGCTGAATATGGCCCTGAATTGAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGt                                                                                                                 <  2:230212/250‑1 (MQ=255)
                                                                                                                                   aTGCGAGTGGTATGGGTGTTGTGAAAGCGCATGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCTGGCAGGGCGTGCATATCCCTTTCCCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTTAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCAGCCCAGttttgta                                                                                                           <  2:193414/251‑4 (MQ=255)
                                                                                                                                              aTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAtttttt                                                                                                >  1:235578/1‑249 (MQ=255)
                                                                                                                                              aTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttt                                                                                                 <  1:22136/250‑2 (MQ=255)
                                                                                                                                                              cgcCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGaa                                                                               <  1:351465/251‑1 (MQ=255)
                                                                                                                                                                   gAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAAACAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACggggaggt                                                                            >  1:243596/1‑244 (MQ=255)
                                                                                                                                                                        tAACGGAACGGCGCTGCCCATTGACCCCTTCCCCGGGTTAGAGCACGGCGGGCAGAGCGTGCATATCCCTTTACCCGAAGGGGACTGGCTGAAGCAGAACCTGAATGTGAATATGGACCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGTCGTAATAGCGCAATGACCTTAACCAGCAACTGGCCGAATCCCAGTTTTTTAGGTGATTTTGTGCCAGCCAACCGGTAAGTTAGCgcg                                                                      <  2:284779/250‑3 (MQ=255)
                                                                                                                                                                            ggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGGGGTACCGGCGATCTTTCTGTGGATCCTGCCGGGCGTAATAGCGAAATTGCCTTAACCAGCAACTGGCCGCATCCCAGGTTTTTAGGTGATTTTTTGCCCGCCAAACGGGAAATTaacggggaggg                                                                 >  2:309966/1‑241 (MQ=255)
                                                                                                                                                                                  gcgcTGACCATTGAACCCGGCTCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGCAGAACCGGAATCTGAATATTGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTTATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCa                                                           <  2:243596/251‑1 (MQ=255)
                                                                                                                                                                                       gACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGgcgc                                                      <  1:93614/251‑1 (MQ=255)
                                                                                                                                                                                        cccATTGAACCCGGCACCGGGTTAGAGCAAGACGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAATCTGAATATGGCCCTTAAGTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCa                                                     <  2:328199/250‑1 (MQ=255)
                                                                                                                                                                                               aaCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCa                                               >  1:109061/1‑250 (MQ=255)
                                                                                                                                                                                                         cGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTGTCCCGAAGGTGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCATGCGCAATGGCAAAGCAGCTGGt                                    <  2:258279/251‑1 (MQ=255)
                                                                                                                                                                                                               agagCAAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCta                              <  1:299413/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCtaataa                          <  1:59510/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGAACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTTGGTGATTTTTTGCACACCAAACGGGAAGGTAGCGCGTCAGGTTTTCCGGCGCAATGGCAAAGCAGCTGGTTTTATAATAATCTCGGTGAGCGTTTTGCTTCAgg  >  2:126504/1‑247 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
AGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGG  >  NC_000913/4638492‑4638977

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