Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,238,232 C→T P200L (CCT→CTT)  yjbH → putative lipoprotein YjbH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,238,2320CT90.9% 127.4 / 5.5 56P200L (CCT→CTT) yjbHputative lipoprotein YjbH
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (30/20);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-02

CAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCGTC  >  NC_000913/4237997‑4238475
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
cAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGCTCCACGGTCTTGCCTCACTGTTTgg                                                                                                                                                                                                                                          <  1:343243/251‑1 (MQ=255)
                          cggcggTACGGTGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACTGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGc                                                                                                                                                                                                                <  2:530601/251‑1 (MQ=255)
                          cggcggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGc                                                                                                                                                                                                                <  1:142154/251‑1 (MQ=255)
                             cggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCTCCAGCGGTAATGTGAAAAATCAGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAACCAGGCGGGCACAATCGAAGGTTACCAGAGGTTCCACGGTCTTGCATCACAGTTTGGCGGCGTggtattacagacgtcctagaaac                                                                                                                                                                                                             >  2:547588/1‑230 (MQ=255)
                              ggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGGTCTATCGACGGTAGCCAGGTGTTCCACGGTCTTGCCTCCCTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAAcc                                                                                                                                                                                                            >  2:266421/1‑251 (MQ=255)
                              ggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAGAACCAGACGCCCTGGCAAcc                                                                                                                                                                                                            >  1:205018/1‑251 (MQ=255)
                              ggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTAGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGGGTGGAATACCAGACGCCCAGGCAAcc                                                                                                                                                                                                            >  1:558941/1‑251 (MQ=255)
                                   ggggCTGTTTGATGCGGAATATCGTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAGCCGCTgc                                                                                                                                                                                                       <  2:139082/251‑1 (MQ=255)
                                    gggCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGTCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTgcg                                                                                                                                                                                                      <  1:632788/251‑1 (MQ=255)
                                     ggCTGATTGATGCGGAATATCTTGTTGACAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGCTAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGCCGGGCGCCAGATGTTCCACGGGATTGCCTCACTGTTTTGCGGCGGGGAATCACAGCCGCCCTGGCAACCGCTgcgc                                                                                                                                                                                                     >  2:613120/1‑251 (MQ=255)
                                       cTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGGTGAATACCAGACGCCCTGGCCACCGCTGCCCCt                                                                                                                                                                                                   >  2:592489/1‑251 (MQ=255)
                                        tGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGttt                                                                                                                                                                                                                                            >  1:392443/1‑209 (MQ=255)
                                        ggTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGttt                                                                                                                                                                                                                                            <  2:392443/208‑1 (MQ=255)
                                         gTTTGATGCGGAATATCTTGCTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACACTCGGCCTGGGCTGGGGGGATTTGGGCACCAGCGGTAATGTGAGAAATCCGCTCTGTT‑CCGCCAGTGATAAATATTACTATCGAGATAACAGATACAAACATGCGAGAGTTATCGACGGAAGACAGCTGCTGCCCGGTCATGCCTCCCTGTTTGCCGGCGGGGAATACCAGACGCCCTGGCAAcccccaccccgt                                                                                                                                                                                                  >  2:411713/1‑240 (MQ=255)
                                                                 cAGCAAAGCCTGGGGGCCGTTCGAGTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAAATCCAGACGCCCTGGCAACCGCTGCGCCTTAAACTGGatgctgaaggaaataaat                                                                                                                                                                         >  2:209856/1‑234 (MQ=255)
                                                                         ccTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATcagcag                                                                                                                                                                 <  1:592489/251‑1 (MQ=255)
                                                                           tGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAgg                                                                                                                                                                >  1:31521/1‑250 (MQ=255)
                                                                            gggggCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGTTATTTGGGCACCAGCGGTAATGTGAAAAAGCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAgg                                                                                                                                                                <  2:643047/249‑1 (MQ=255)
                                                                              gggCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGAttt                                                                                                                                                            >  1:542617/1‑251 (MQ=255)
                                                                                        aTTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGa                                                                                                                                                    >  1:316395/1‑249 (MQ=255)
                                                                                           ttACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTg                                                                                                                                               >  1:261251/1‑251 (MQ=255)
                                                                                             aCGCTCGGCTTGGGCTGGGTGGATTTGGTCACCAGCGGTAATGTGAAAAATCCGCT‑TGTTTCAGCCTGTGATAAAGATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑Ga                                                                                                                                            <  2:27022/251‑1 (MQ=255)
                                                                                                            tgggGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTTTCGAAGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGCATAACAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGTGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGCAGCTG‑GGGAAAAAAAGT‑AAGtt                                                                                                                            >  2:299296/1‑251 (MQ=255)
                                                                                                                ggTATTTGGGCAACAGCGGTAAGGTTAAAAAACCGCTCTGTT‑CAGCCATTGTTAAATATGGCTATCTCGATAACAGCTAACAACAGGCGGGACCTATCGACGTGAGCCAGATGTGCCACGGTCTTGCCTCATTGTTTGTCGGTGTGGAATACCAGACGCCCTAGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAAc                                                                                                                        <  2:558941/251‑1 (MQ=255)
                                                                                                                        ggCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGCGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGA‑AAATTTAACGTCggagc                                                                                                                >  1:408181/1‑248 (MQ=255)
                                                                                                                              aGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAATTTAACGTCGG‑GGCGGttt                                                                                                          >  1:368497/1‑250 (MQ=255)
                                                                                                                               gCGGGAATGTGAAAAATCCGCGCTGTT‑CAGCCAGTGATAAATATTGCTAACGCGTTAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGGTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTAt                                                                                                        <  2:26410/251‑1 (MQ=255)
                                                                                                                               gCGGAAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGGTTTTGCCGGGAAGCTG‑GGGCAAAAAAAT‑AAGTTTAACGTCGG‑GGCGCTTTAt                                                                                                        >  2:510003/1‑251 (MQ=255)
                                                                                                                                 ggTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGATATCGCGATAACAGCTACAAACAGGCGGGATCTATAGAAGGAAGCCAGATGTTCACCGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACAGCTGCGCCTGAAACTGGAGGATGAAGGCAATAATTATAAGCAGGATTTTGCCGGGAAGCTG‑GAGAAAAAAAGTCAAGTT‑AAAGTCGG‑TGCGGTTTAt                                                                                                        >  2:386152/1‑249 (MQ=255)
                                                                                                                                  gTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGGTAGATGTTCCACGGTCTTGCCTCTCTGTTTGGCGGCGTGGAATACCAGACGCCCGGGCAACCGCGGCGCCTGAAACTGTAGTAGGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATcg‑c                                                                                                    <  1:95619/251‑1 (MQ=255)
                                                                                                                                     aTGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAACCAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGAAGCCCTGGCAACCGCTGCGCCTGATACTGGAGTATGAAGGCAATCATTATCAGCAGGATTTTGACGGGAAGCTG‑GAGCAAAAAAAT‑AAGTTTAACGTCGG‑GGCGaagatcgggagg                                                                                                  >  2:601342/1‑240 (MQ=255)
                                                                                                                                     aTGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGa                                                                                                             <  1:601342/240‑1 (MQ=255)
                                                                                                                                       gtgAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTAc                                                                                               >  1:309818/1‑251 (MQ=255)
                                                                                                                                       gtgAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACAGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGATG‑GAGCAAAAAAGT‑AAGTTTAACGTCGGGTG‑GATTTATcggggta                                                                                                 >  1:356242/1‑245 (MQ=255)
                                                                                                                                                acGCTCTGTT‑CAGCCAGGTATAAATATTGCTATCGCGATAACAGCTACAAACAGTCAGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGTTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCTGGGAAGATG‑GATCAAAAAAGT‑AAGTTTGACGTCGG‑TGCGATTTATCG‑CGTTCCCGATTggg                                                                                       <  1:229843/249‑1 (MQ=255)
                                                                                                                                                 cGCTCTGTC‑CAGCCAGGGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGcc                                                                                     <  1:427002/251‑1 (MQ=255)
                                                                                                                                                    tctGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGTCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAACACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGAc                                                                                  >  1:37381/1‑251 (MQ=255)
                                                                                                                                                     ttatt‑CAGCCAGTGATAACTATTGCTATCGCGATAACAGCTACAAACAGGCGGGATTTATCGACGGTCGCCAAATGTTCCACGGTCTTGACTCACTGTTTGGTGGCGGGTAATACCAGACGCCCGGGCAACCGCTGCGCCTGAAACTTGAGAATAAAAGCCATACTTATCCGCAGGACTTTGCCGGGAAACTGCGA‑CAAACACGA‑AAGTTTAACATCGG‑TTCTATTTATCG‑CGTTACCGagagggtccaac                                                                                 >  2:225469/4‑240 (MQ=255)
                                                                                                                                                            aGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTAACGACGGTAGCCAGATGATCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCCGACGCCCTGGCAACCGCTGCGCCTGAAAATGGAGTATGAAGGCAATAATTAGCAGCAGGATTTTGCCGTGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑TGTTACCGATTGGGCCGACGTTAAcc                                                                           >  1:132729/1‑251 (MQ=255)
                                                                                                                                                                       aaTATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATAACAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCT‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGa                                                                 >  1:196662/1‑250 (MQ=255)
                                                                                                                                                                          atTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCAAGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGGAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACg                                                              <  1:324079/250‑1 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTACACACAGGCGGGATCTATCGCCGGTTACCATATGTTCCACGGGCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATAAAGGAAATAATTATAAGTAGGGTTTTGCCGGGAAGCTG‑GAGAAAAAAAGT‑TAGTTTAACGTCGG‑GGGGATTTATCGGCGT‑ACCGGTTGGGCCggagtgtaactttgatatggatgtgg                                                           >  2:445175/1‑225 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGGTTGGGCCGACGTTAACCTTAGATATGAACGTggg                                                          >  2:240267/1‑250 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGc                                                          <  1:240267/251‑1 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTAAAAACAGGTGGGATTTATGCACGGTAGCTAGATGTTCCAATGTCTTGCCTCACTGTTTGGCGGCGATGAATTCCAGACGACTTGGCAACCGCTGCGCCTGAAACGGGAGTATGAAGGCAATAATTATCAGAAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGAGATTTATCG‑CGTTACGGATTGTGTCGACGTTAACCTTAGTTATGTACGTGGc                                                          <  2:132729/251‑1 (MQ=255)
                                                                                                                                                                                   gcgATAACAGCTACAAACAGGCGGGAGCTATCGACGGTAGCCAGAGGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGGGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑ACTTACCCATTGGGCCGACGTTCAccattgctaataaaggggcatcacc                                                    >  2:248492/1‑228 (MQ=255)
                                                                                                                                                                                       tAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGGTTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCtttt                                                >  1:142209/1‑250 (MQ=255)
                                                                                                                                                                                                     aGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGA‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTCACg                                  >  1:186024/1‑251 (MQ=255)
                                                                                                                                                                                                      ggCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGaa                                                                <  1:69660/220‑1 (MQ=255)
                                                                                                                                                                                                              cTATCGAAGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGCAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACATTAGTTATGAACGTTGCAACACCTTTATGTTTGGCGTTACGTTgcgcccc                         >  2:604517/1‑248 (MQ=255)
                                                                                                                                                                                                                   gACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACtt                    >  1:203930/1‑251 (MQ=255)
                                                                                                                                                                                                                   gACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACtt                    >  1:555380/1‑251 (MQ=255)
                                                                                                                                                                                                                     cGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTAGGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTa                  <  1:386152/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGAt             <  1:179502/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGAt             <  1:299296/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCGTc  >  1:581307/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
CAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTT‑CAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTG‑GAGCAAAAAAGT‑AAGTTTAACGTCGG‑TGCGATTTATCG‑CGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCGTC  >  NC_000913/4237997‑4238475

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: ‑