Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 488,403 C→T T623I (ACC→ATC)  mscK → potassium dependent, small conductance mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913488,4030CT95.3% 122.9 / NA 43T623I (ACC→ATC) mscKpotassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/18);  minor base G (0/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAACGACAGCCAGCTCAACACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGCGACTGGTCACCATTACCGTGCTGTCGAT  >  NC_000913/488160‑488637
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
gTAACGACAGCCAGCTCAACACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGCGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGCCCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCAc                                                                                                                                                                                                                                     >  2:85804/1‑251 (MQ=255)
  aaCGACGGCCAGCTCAACACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCCCGCGGTTGGCCGGATTCTGTTGACCATGCAGCGCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCGGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACt                                                                                                                                                                                                                                    <  2:201617/250‑1 (MQ=255)
    cGACAGCCAGCTCAAAACACCAACAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGACCTGCTATGGTCGTTCCGCAAAAAACTGGCGCTATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATAAGCCACTGGc                                                                                                                                                                                                                                 >  2:230281/1‑251 (MQ=255)
               tCAACACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAAtt                                                                                                                                                                                                                      >  1:143104/1‑251 (MQ=255)
                cAACCCACCAAAAGCGATCGTTAACGACCGGACCCGGGCGCTGCCGGTGTGCCTGATTATTCGCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGTTATGGTCATTCAGCAAAAAACTGCCGATATTCTGGCTGGTGGTTCGCGTGTGCTGGGAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCTCTTCGGCAGGCCGTAACAGCAGATCAGCCACTGGCGTCGGCAAAtt                                                                                                                                                                                                                      <  2:153853/250‑1 (MQ=255)
                 aaCACCCCAAAAGCGGTCCTTATCGACCTGATACGGGCGCTGCCGGAGTGCCTGATTATTCTCGCGGTTGGCCTGATTATGTTGACAATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATCTTCTGGCTGGGGTATGGCCTGGGCTGGCAGGTACTGGAGAAACACGGCGTTGCCGTCAGTCACTTCGGCATGCCGGAACAGCAGTTCAGCAACTGGCGTCGGCACATTGt                                                                                                                                                                                                                    >  2:164349/1‑251 (MQ=255)
                   cacacCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGCGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTcc                                                                                                                                                                                                                  >  2:246803/1‑251 (MQ=255)
                        cAAAAGCGATCCTTATCGACCTGAACCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGGTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTCTGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCCTCGGCAAATTGTCCGCATc                                                                                                                                                                                                             >  2:48245/1‑251 (MQ=255)
                                  ccTTAGCGACTTGATCAGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGTCGAGAGTCTGGCTGGTGTTTGGCCTGGGCTGGAAGGTACTGTATAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCTCATCAGTCTCGCa                                                                                                                                                                                                    <  2:129945/250‑1 (MQ=255)
                                   cTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCCGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCAtt                                                                                                                                                                                                  >  1:252111/1‑251 (MQ=255)
                                     tATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATtgc                                                                                                                                                                                                >  1:232196/1‑251 (MQ=255)
                                             tGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCACAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGACGGTACTGGAGACAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCCTCAGTCTCGCATTGCTGCCTATc                                                                                                                                                                                        >  1:256305/1‑251 (MQ=255)
                                              gATCCGTGCTCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGTACAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAGCTGGCGATAGTCTGGCTGGTGTTTGGCCGGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATTCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCACTCTCGCATTGCTGCCTATc                                                                                                                                                                                        <  2:37316/250‑1 (MQ=255)
                                               aTCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATcc                                                                                                                                                                                       >  1:180332/1‑250 (MQ=255)
                                                 ccGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGGGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAGCAGCAGATCAGCCACTGGCTTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCAt                                                                                                                                                                                     <  2:256305/250‑1 (MQ=255)
                                                            cGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGAACATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCCAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGAAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATACATTTCTGGTCtgtg                                                                                                                                                                         >  2:226483/1‑251 (MQ=255)
                                                                                           gATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTgatgg                                                                                                                                          <  1:7575/251‑1 (MQ=255)
                                                                                                                         cGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTttct                                                                                                            >  1:42173/1‑251 (MQ=255)
                                                                                                                               gCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGTTGATTTTCTTCAAcc                                                                                                      >  2:249181/1‑251 (MQ=255)
                                                                                                                                  aTGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCACCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACctgc                                                                                                   >  1:159759/1‑251 (MQ=255)
                                                                                                                                            aGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACATGCTGCTGATTGc                                                                                         <  1:53109/251‑1 (MQ=255)
                                                                                                                                             gCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGcc                                                                                        <  1:201892/251‑1 (MQ=255)
                                                                                                                                              cAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGAGCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCt                                                                                       <  2:57191/251‑1 (MQ=255)
                                                                                                                                                 aaaaCTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTcc                                                                                    >  1:227866/1‑251 (MQ=255)
                                                                                                                                                        gCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATg                                                                             >  1:34963/1‑251 (MQ=255)
                                                                                                                                                          gATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCGGTCTCGCATGGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGc                                                                           <  1:156791/251‑1 (MQ=255)
                                                                                                                                                             atTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCg                                                                         >  1:123607/1‑250 (MQ=255)
                                                                                                                                                                  ggctggTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGAtgt                                                                     <  2:78910/249‑1 (MQ=255)
                                                                                                                                                                       gtgtTTGGCCTGTGCTGGAAGGTACGGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGa                                                              <  2:143104/251‑1 (MQ=255)
                                                                                                                                                                                  gtgCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCg                                                     <  1:66101/249‑1 (MQ=255)
                                                                                                                                                                                   tgCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTg                                                   >  1:23438/1‑250 (MQ=255)
                                                                                                                                                                                      tGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAACTTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATa                                                >  1:245447/1‑250 (MQ=255)
                                                                                                                                                                                                gggAGAAAAACGGCGTTGCCGTACGACACTTCGGCGGGCCGGAACAGCAGAGCAGCCACTGGCGTCGTCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCGGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTGCTTCAACCTGCTGCTGAATGCCTTCCTGGGATGGCCGATGTGCCGCGAATGCTGGCGTGATAAAGAATCGcac                                     <  2:109729/250‑1 (MQ=255)
                                                                                                                                                                                                     aaaaaCGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCAt                                <  2:123607/251‑1 (MQ=255)
                                                                                                                                                                                                       aaaCGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGc                              >  1:174655/1‑251 (MQ=255)
                                                                                                                                                                                                         aCGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGAtgtg                                                                    <  1:77651/211‑1 (MQ=255)
                                                                                                                                                                                                         aCGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGAtgtg                                                                    >  2:77651/1‑211 (MQ=255)
                                                                                                                                                                                                           ggCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGCGACt                          <  1:249181/251‑1 (MQ=255)
                                                                                                                                                                                                             cGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATGTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGCGACTgg                        <  1:85804/251‑1 (MQ=255)
                                                                                                                                                                                                                   cgtccctCACTCCGGCTTGCCTGAACAGCAGTTCACCTACTGGGGTCGGCAGATTGTCCGCAGCAGTCTCGCATTGCTGCCTATCCATTTATGGTCTGTGGTGGCGGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTGTCTCCAACCTGGTGCTGATTGCCTGCCGGGTATGGCCGATTTGCCGCGAAAGCTGGCGTGATAAATCGGCGCACACCGTGCGACTggactcca                  <  2:129706/245‑7 (MQ=255)
                                                                                                                                                                                                                      aCGTCACTTCGGCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGCGACTGGTCACCATTa               <  1:111364/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggCATGCCGGAACAGCAGATCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGCGACTGGTCACCATTACCGTGCTGTc     >  2:190837/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tGCCGGAACAGCAGATAAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGCGACTGGTCACCATTACCGTGCTGTCGAt  >  1:91988/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
GTAACGACAGCCAGCTCAACACACCAAAAGCGATCCTTATCGACCTGATCCGTGCGCTGCCGGTGTGCCTGATTATTCTCGCGGTTGGCCTGATTCTGTTGACCATGCAGCTCAACATCAGCGAACTGCTATGGTCGTTCAGCAAAAAACTGGCGATATTCTGGCTGGTGTTTGGCCTGTGCTGGAAGGTACTGGAGAAAAACGGCGTTGCCGTACGTCACTTCGGCATGCCGGAACAGCAGACCAGCCACTGGCGTCGGCAAATTGTCCGCATCAGTCTCGCATTGCTGCCTATCCATTTCTGGTCTGTGGTGGCAGAACTTTCCCCGCTGCATCTGATGGATGATGTGCTGGGGCAAGCGATGATTTTCTTCAACCTGCTGCTGATTGCCTTCCTGGTATGGCCGATGTGCCGCGAAAGCTGGCGTGATAAAGAGTCGCACACCATGCGACTGGTCACCATTACCGTGCTGTCGAT  >  NC_000913/488160‑488637

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