Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,600 A→C P141P (CCA→CCC creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,6000AC83.8% 84.4 / 7.3 37P141P (CCA→CCCcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base C (17/14);  total (20/17)
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.42e-01

ATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCT  >  NC_000913/4638380‑4638841
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
ataCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGcgcg                                                                                                                                                                                                                        <  1:308408/250‑1 (MQ=255)
  aCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGAGGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTg                                                                                                                                                                                                                      <  1:8672/250‑1 (MQ=255)
             agagGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATtggtgtgg                                                                                                                                                                                                          >  1:251194/1‑251 (MQ=255)
                    aaaaCCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTGTTGTGATTAGCGTCGGGGATGCGCGTGGTAGTGGTGTGGTGAAAg                                                                                                                                                                                                    <  1:135621/250‑1 (MQ=255)
                                 ggAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaac                                                                                                                                                                                      >  1:58784/1‑251 (MQ=255)
                                                 ccagttgtGGTTGCCTGAGTCATTGATTGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGGCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAcc                                                                                                                                                                       <  2:378130/243‑1 (MQ=255)
                                                          gTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCg                                                                                                                                                             <  1:100345/251‑1 (MQ=255)
                                                               cTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCAcc                                                                                                                                                        >  1:58432/1‑251 (MQ=255)
                                                                tGAGTCATTGATGGTTAATGGCCATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTAGCACAGGGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCAAAATATCATCTTAGGCAAGCCCTTTATTGTAATTAGCGTCGGGGATACGCGAGGTATTGGTGTGGTGAAAGCGCATGAAGTTAACAGAACGGGGCTGaaattttaagccggcgccg                                                                                                                                                       >  2:250878/1‑233 (MQ=255)
                                                                                gatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCgg                                                                                                                                       <  1:369332/251‑1 (MQ=255)
                                                                                 atggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTAAGGTCTGGCACAGTGATTTACCGCTAAACGCCGATTTAGATGTTTCGAGTCTTAGCGAGCTCAACGCGCCACATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGGGGAGGAAGCGCCTGAAATTAACGGAACGGCGCCGACGATTTAACCCGGCACCGGGGTAGAGCAAAGCggg                                                                                                                                      >  2:384838/1‑251 (MQ=255)
                                                                                      aaTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAggg                                                                                                                                 >  1:230212/1‑251 (MQ=255)
                                                                                             aCGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGAGCCAAATATCATCTTAGGCAAGACCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGAGATGACCATTGCACCCGACACCGGTTTAGAGCCAGGCGGGCAggggggggat                                                                                                                          >  2:259013/1‑244 (MQ=255)
                                                                                                          cAAGTTAGTGATTTATGCCGGTCAGGTCTGGTACAGTGATTTAACGCTAAAAACTGAGTTCGATGTTTCGCGTCTTAGCGAGCTCAACTCGCCAAATATCATCTTAGGCAAGCCCCTTATTGTGATTATCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCGGAAGTTAACGGAACTGCGCTGACCATTGAACCCGGCTCCGGGTTAGAGCAAGGCTGGCAGGGCGTGCATATCCCttttcccg                                                                                                             <  2:102508/251‑6 (MQ=255)
                                                                                                                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)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
ATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCT  >  NC_000913/4638380‑4638841

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