Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,893,219 C→T S188F (TCC→TTC)  ygcN → putative oxidoreductase with FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,893,2190CT98.2% 172.2 / NA 54S188F (TCC→TTC) ygcNputative oxidoreductase with FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/25);  minor base G (0/1);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GGCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATCCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTA  >  NC_000913/2892976‑2893465
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
ggCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTAt                                                                                                                                                                                                                                                  <  1:140577/250‑1 (MQ=255)
    ggtggccgaagccgaaAAAGAAGGTGTCGAATGCAGCCCCGGTGCTGCGGGGGATGCACTGTATGAAGAAAACGGCAGAGTCTGGGGCGTTATTTGTGGTAACGATATTCTCCGCGCCCGTGATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGTAACATGCGCTATTgaag                                                                                                                                                                                                                                             <  2:63925/247‑1 (MQ=255)
    ggTTGCCGAAGCCGAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTga                                                                                                                                                                                                                                               >  2:409283/1‑249 (MQ=255)
            aagctgaaAAAGAAGGTGTCGAATGCATCCCCGGTGCGACTGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGtttt                                                                                                                                                                                                                                     <  1:197342/251‑1 (MQ=255)
             agccgaaAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTc                                                                                                                                                                                                                                    >  1:129752/1‑251 (MQ=255)
               ccgaaAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCAt                                                                                                                                                                                                                                  >  1:110938/1‑251 (MQ=255)
                   aaaaGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGGATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTgg                                                                                                                                                                                                                              <  2:180183/251‑1 (MQ=255)
                         agGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACGGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATa                                                                                                                                                                                                                        <  2:18409/251‑1 (MQ=255)
                                  aTGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCACCAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGgcg                                                                                                                                                                                                               >  2:3988/1‑251 (MQ=255)
                                        ccccGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGGCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCg                                                                                                                                                                                                          <  1:270394/250‑1 (MQ=255)
                                              aGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGt                                                                                                                                                                                                   >  1:39050/1‑251 (MQ=255)
                                                 gACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGCGACGATATTCTCCGCGCCCGCTATGTGGTGCTGGCAGAAGGTGCCCACCGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGACGCGCTGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGATGCTGTTct                                                                                                                                                                                                >  1:434499/1‑250 (MQ=255)
                                                     gTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCgg                                                                                                                                                                                            >  1:473302/1‑251 (MQ=255)
                                                           gCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGa                                                                                                                                                                                       <  1:2511/250‑1 (MQ=255)
                                                            cACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATc                                                                                                                                                                                     >  1:374583/1‑251 (MQ=255)
                                                             aCTGTAGGAAGAAAACGGCAGTGTCTGGTGCGTTAGTTGTGGTGACGATATTCTCCGCGCCCGGTATGTGGGTGTGGCAGAAGGTGCCAACAGCGTTCTGGCTGAACGTCACGTGTTAGGGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTTGAGAATAACGAATGCGCAGCGTTGCTGTCCAGCGGCAGGATCt                                                                                                                                                                                    <  2:204653/251‑1 (MQ=255)
                                                               tGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCg                                                                                                                                                                                             <  2:206137/240‑1 (MQ=255)
                                                               tGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCg                                                                                                                                                                                             >  1:206137/1‑240 (MQ=255)
                                                                     aagaaAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGAc                                                                                                                                                                            >  1:325642/1‑251 (MQ=255)
                                                                        aaaaCGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTATTCAGCGGCAGGGTCTGGGATGACtt                                                                                                                                                                          >  2:329026/1‑250 (MQ=255)
                                                                         aaaCGGCAGAGTGTGTGGCGTTATGTGTGGGGACGATATTCACCGCGCTCGTTATGTGGTGCTGGCAGAAGGTGCCAACCGCGGCCTGGCGGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGGTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTAGTGAAGATCGTTTTCATCTGGATAATAACGAAGGCGCAGCTTTGCTGTTCAGCGGCAGGATCTGTGATGACTTAc                                                                                                                                                                        <  2:431749/251‑1 (MQ=255)
                                                                                  agTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTTCTGCCGGCGAAGCGATGGCGGTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTACAGCGGCAGGATCTGTGATGACTTACCCggcggc                                                                                                                                                                <  2:182/250‑1 (MQ=255)
                                                                                          gCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCtt                                                                                                                                                       >  1:12618/1‑251 (MQ=255)
                                                                                           cGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCttt                                                                                                                                                      >  1:217515/1‑251 (MQ=255)
                                                                                                     ggtgACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGTATCTGTGATGACTTACTCGGCGGCGCCTTTCTTTTTACTAATc                                                                                                                                             >  2:323107/1‑250 (MQ=255)
                                                                                                          cGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTACTcaacaa                                                                                                                                        >  1:478404/1‑250 (MQ=255)
                                                                                                           gATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGACGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTAGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCACAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCGGCGTTGCTGTTCAGCGGCAGGGTCTGGGATGACTTACCCGGCGGCGCATTTCTTTATACTCATCAACAAAc                                                                                                                                      >  2:53639/1‑251 (MQ=255)
                                                                                                                           gTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCGCTCGTTAGGGAtt                                                                                                                      <  1:429307/251‑1 (MQ=255)
                                                                                                                            ttATGTGGTGCTGGCAGAAGTTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACGTACCCGGCGGCGCATTTCTTTTTACTAATCAACAAACGCTCTCGTTAGGGATTg                                                                                                                     >  2:264745/1‑251 (MQ=255)
                                                                                                                             tATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTg                                                                                                                     <  1:132493/250‑1 (MQ=255)
                                                                                                                                                    ccAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACg                                                                                              <  1:3988/251‑1 (MQ=255)
                                                                                                                                                                 ggCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTc                                                                                 >  1:476182/1‑251 (MQ=255)
                                                                                                                                                                    tGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCgg                                                                              >  2:422167/1‑251 (MQ=255)
                                                                                                                                                                             cGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTAGACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAgc                                                                     >  1:159381/1‑251 (MQ=255)
                                                                                                                                                                                   aGTGACTCGTCCTTCTGGCGAAGTGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGa                                                               <  2:342098/251‑1 (MQ=255)
                                                                                                                                                                                       aCTCGTCCTGCTGGCGAAGCGATTGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCg                                                           <  2:484623/251‑1 (MQ=255)
                                                                                                                                                                                        cTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGc                                                          >  2:71006/1‑251 (MQ=255)
                                                                                                                                                                                                gctgGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTaaa                                                    <  1:483885/249‑1 (MQ=255)
                                                                                                                                                                                                gctgGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAg                                                   <  1:409283/250‑1 (MQ=255)
                                                                                                                                                                                                 ctgGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAgcg                                                 >  1:414495/1‑251 (MQ=255)
                                                                                                                                                                                                   ggCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCCGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCAGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGATGACTCGCTTTAAAgcgc                                                >  2:238637/1‑250 (MQ=255)
                                                                                                                                                                                                   ggCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCAGGAAACATTCGCTATTGAAGAACGTTGTCATCTGTAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGTCAAGCGAGCTGCTGACTCGCTTTAAAgcg                                                 <  2:284828/249‑1 (MQ=255)
                                                                                                                                                                                                      gAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGCCTCGCTTTAAAGCGCATcc                                            >  2:168281/1‑251 (MQ=255)
                                                                                                                                                                                                              ggCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTgc                                    >  2:460630/1‑251 (MQ=255)
                                                                                                                                                                                                                 gTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGcc                                 <  1:422270/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggggATCAAAGAAGTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCACATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCg                               <  1:70540/250‑1 (MQ=255)
                                                                                                                                                                                                                        gtCAAAGAAGTTCTGCCGCGGGAAACATTCGCTATTGAAGACCGTTTTCAGCTGGAGAAAAACGAAGGCGCAGCGTTGCTGTTCGGCGGCAGGATCTGTGATGTCTTACCCGGCGGCGCATTTCTTTATACTAATCAACAACCGCTCTCGTTAGGGATTGTTGGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTAt                          <  2:80413/250‑1 (MQ=255)
                                                                                                                                                                                                                              gaagTGCTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCaaaa                     >  2:11364/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ggCTGGCGCTGGAAACATTCGCTAGTGAAGAACGTTTGCAGCTGGAGAATAACGAAGGCGCAGGGTTGCGGTTCAGCGGCAGGATCTGTGATGACGTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCCCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGGGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTAGCAAAAACACg                <  2:434499/250‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTg                                                                                                                     <  2:377101/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTg                                                                                                                     >  1:377101/1‑148 (MQ=255)
                                                                                                                                                                                                                                      gTCGCTGGAAACATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGa              <  1:460630/249‑1 (MQ=255)
                                                                                                                                                                                                                                              aaaCATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTgg      >  1:292954/1‑249 (MQ=255)
                                                                                                                                                                                                                                                aCATTCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGTTCGGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAATCGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTTCGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTa  <  2:330689/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
GGCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGAAACATCCGCTATTGAAGAACGTTTTCATCTGGAGAATAACGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTA  >  NC_000913/2892976‑2893465

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