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

CAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCGTC  >  NC_000913/4237997‑4238475
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
cAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGCTCCACGGTCTTGCCTCACTGTTTgg                                                                                                                                                                                                                                          <  1:343243/251‑1 (MQ=255)
                          cggcggTACGGTGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACTGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGc                                                                                                                                                                                                                <  2:530601/251‑1 (MQ=255)
                          cggcggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGc                                                                                                                                                                                                                <  1:142154/251‑1 (MQ=255)
                             cggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCTCCAGCGGTAATGTGAAAAATCAGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAACCAGGCGGGCACAATCGAAGGTTACCAGAGGTTCCACGGTCTTGCATCACAGTTTGGCGGCGTggtattacagacgtcctagaaac                                                                                                                                                                                                             >  2:547588/1‑230 (MQ=255)
                              ggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGGTCTATCGACGGTAGCCAGGTGTTCCACGGTCTTGCCTCCCTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAAcc                                                                                                                                                                                                            >  2:266421/1‑251 (MQ=255)
                              ggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAGAACCAGACGCCCTGGCAAcc                                                                                                                                                                                                            >  1:205018/1‑251 (MQ=255)
                              ggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTAGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGGGTGGAATACCAGACGCCCAGGCAAcc                                                                                                                                                                                                            >  1:558941/1‑251 (MQ=255)
                                   ggggCTGTTTGATGCGGAATATCGTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAGCCGCTgc                                                                                                                                                                                                       <  2:139082/251‑1 (MQ=255)
                                    gggCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGTCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTgcg                                                                                                                                                                                                      <  1:632788/251‑1 (MQ=255)
                                     ggCTGATTGATGCGGAATATCTTGTTGACAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGCTAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGCCGGGCGCCAGATGTTCCACGGGATTGCCTCACTGTTTTGCGGCGGGGAATCACAGCCGCCCTGGCAACCGCTgcgc                                                                                                                                                                                                     >  2:613120/1‑251 (MQ=255)
                                       cTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGGTGAATACCAGACGCCCTGGCCACCGCTGCCCCt                                                                                                                                                                                                   >  2:592489/1‑251 (MQ=255)
                                        tGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGttt                                                                                                                                                                                                                                            >  1:392443/1‑209 (MQ=255)
                                        ggTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGttt                                                                                                                                                                                                                                            <  2:392443/208‑1 (MQ=255)
                                         gTTTGATGCGGAATATCTTGCTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACACTCGGCCTGGGCTGGGGGGATTTGGGCACCAGCGGTAATGTGAGAAATCCGCTCTGTTCCGCCAGTGATAAATATTACTATCGAGATAACAGATACAAACATGCGAGAGTTATCGACGGAAGACAGCTGCTGCCCGGTCATGCCTCCCTGTTTGCCGGCGGGGAATACCAGACGCCCTGGCAAcccccaccccgt                                                                                                                                                                                                  >  2:411713/1‑240 (MQ=255)
                                                                 cAGCAAAGCCTGGGGGCCGTTCGAGTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAAATCCAGACGCCCTGGCAACCGCTGCGCCTTAAACTGGatgctgaaggaaataaat                                                                                                                                                                         >  2:209856/1‑234 (MQ=255)
                                                                         ccTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATcagcag                                                                                                                                                                 <  1:592489/251‑1 (MQ=255)
                                                                           tGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAgg                                                                                                                                                                >  1:31521/1‑250 (MQ=255)
                                                                            gggggCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGTTATTTGGGCACCAGCGGTAATGTGAAAAAGCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAgg                                                                                                                                                                <  2:643047/249‑1 (MQ=255)
                                                                              gggCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGAttt                                                                                                                                                            >  1:542617/1‑251 (MQ=255)
                                                                                        aTTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGa                                                                                                                                                    >  1:316395/1‑249 (MQ=255)
                                                                                           ttACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTg                                                                                                                                               >  1:261251/1‑251 (MQ=255)
                                                                                             aCGCTCGGCTTGGGCTGGGTGGATTTGGTCACCAGCGGTAATGTGAAAAATCCGCTTGTTTCAGCCTGTGATAAAGATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGa                                                                                                                                            <  2:27022/251‑1 (MQ=255)
                                                                                                            tgggGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTTTCGAAGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGCATAACAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGTGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGCAGCTGGGGAAAAAAAGTAAGtt                                                                                                                            >  2:299296/1‑251 (MQ=255)
                                                                                                                ggTATTTGGGCAACAGCGGTAAGGTTAAAAAACCGCTCTGTTCAGCCATTGTTAAATATGGCTATCTCGATAACAGCTAACAACAGGCGGGACCTATCGACGTGAGCCAGATGTGCCACGGTCTTGCCTCATTGTTTGTCGGTGTGGAATACCAGACGCCCTAGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAAc                                                                                                                        <  2:558941/251‑1 (MQ=255)
                                                                                                                        ggCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGCGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGAAAATTTAACGTCggagc                                                                                                                >  1:408181/1‑248 (MQ=255)
                                                                                                                              aGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAATTTAACGTCGGGGCGGttt                                                                                                          >  1:368497/1‑250 (MQ=255)
                                                                                                                               gCGGGAATGTGAAAAATCCGCGCTGTTCAGCCAGTGATAAATATTGCTAACGCGTTAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGGTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTAt                                                                                                        <  2:26410/251‑1 (MQ=255)
                                                                                                                               gCGGAAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGGTTTTGCCGGGAAGCTGGGGCAAAAAAATAAGTTTAACGTCGGGGCGCTTTAt                                                                                                        >  2:510003/1‑251 (MQ=255)
                                                                                                                                 ggTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGATATCGCGATAACAGCTACAAACAGGCGGGATCTATAGAAGGAAGCCAGATGTTCACCGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACAGCTGCGCCTGAAACTGGAGGATGAAGGCAATAATTATAAGCAGGATTTTGCCGGGAAGCTGGAGAAAAAAAGTCAAGTTAAAGTCGGTGCGGTTTAt                                                                                                        >  2:386152/1‑249 (MQ=255)
                                                                                                                                  gTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGGTAGATGTTCCACGGTCTTGCCTCTCTGTTTGGCGGCGTGGAATACCAGACGCCCGGGCAACCGCGGCGCCTGAAACTGTAGTAGGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATcgc                                                                                                    <  1:95619/251‑1 (MQ=255)
                                                                                                                                     aTGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAACCAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGAAGCCCTGGCAACCGCTGCGCCTGATACTGGAGTATGAAGGCAATCATTATCAGCAGGATTTTGACGGGAAGCTGGAGCAAAAAAATAAGTTTAACGTCGGGGCGaagatcgggagg                                                                                                  >  2:601342/1‑240 (MQ=255)
                                                                                                                                     aTGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGa                                                                                                             <  1:601342/240‑1 (MQ=255)
                                                                                                                                       gtgAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTAc                                                                                               >  1:309818/1‑251 (MQ=255)
                                                                                                                                       gtgAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACAGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGATGGAGCAAAAAAGTAAGTTTAACGTCGGGTGGATTTATcggggta                                                                                                 >  1:356242/1‑245 (MQ=255)
                                                                                                                                                acGCTCTGTTCAGCCAGGTATAAATATTGCTATCGCGATAACAGCTACAAACAGTCAGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGTTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCTGGGAAGATGGATCAAAAAAGTAAGTTTGACGTCGGTGCGATTTATCGCGTTCCCGATTggg                                                                                       <  1:229843/249‑1 (MQ=255)
                                                                                                                                                 cGCTCTGTCCAGCCAGGGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGcc                                                                                     <  1:427002/251‑1 (MQ=255)
                                                                                                                                                    tctGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGTCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAACACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGAc                                                                                  >  1:37381/1‑251 (MQ=255)
                                                                                                                                                     ttattCAGCCAGTGATAACTATTGCTATCGCGATAACAGCTACAAACAGGCGGGATTTATCGACGGTCGCCAAATGTTCCACGGTCTTGACTCACTGTTTGGTGGCGGGTAATACCAGACGCCCGGGCAACCGCTGCGCCTGAAACTTGAGAATAAAAGCCATACTTATCCGCAGGACTTTGCCGGGAAACTGCGACAAACACGAAAGTTTAACATCGGTTCTATTTATCGCGTTACCGagagggtccaac                                                                                 >  2:225469/4‑240 (MQ=255)
                                                                                                                                                            aGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTAACGACGGTAGCCAGATGATCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCCGACGCCCTGGCAACCGCTGCGCCTGAAAATGGAGTATGAAGGCAATAATTAGCAGCAGGATTTTGCCGTGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGTGTTACCGATTGGGCCGACGTTAAcc                                                                           >  1:132729/1‑251 (MQ=255)
                                                                                                                                                                       aaTATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATAACAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCTCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGa                                                                 >  1:196662/1‑250 (MQ=255)
                                                                                                                                                                          atTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCAAGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGGAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACg                                                              <  1:324079/250‑1 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTACACACAGGCGGGATCTATCGCCGGTTACCATATGTTCCACGGGCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATAAAGGAAATAATTATAAGTAGGGTTTTGCCGGGAAGCTGGAGAAAAAAAGTTAGTTTAACGTCGGGGGGATTTATCGGCGTACCGGTTGGGCCggagtgtaactttgatatggatgtgg                                                           >  2:445175/1‑225 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGGTTGGGCCGACGTTAACCTTAGATATGAACGTggg                                                          >  2:240267/1‑250 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGc                                                          <  1:240267/251‑1 (MQ=255)
                                                                                                                                                                             gCTATCGCGATAACAGCTAAAAACAGGTGGGATTTATGCACGGTAGCTAGATGTTCCAATGTCTTGCCTCACTGTTTGGCGGCGATGAATTCCAGACGACTTGGCAACCGCTGCGCCTGAAACGGGAGTATGAAGGCAATAATTATCAGAAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGAGATTTATCGCGTTACGGATTGTGTCGACGTTAACCTTAGTTATGTACGTGGc                                                          <  2:132729/251‑1 (MQ=255)
                                                                                                                                                                                   gcgATAACAGCTACAAACAGGCGGGAGCTATCGACGGTAGCCAGAGGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGGGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGACTTACCCATTGGGCCGACGTTCAccattgctaataaaggggcatcacc                                                    >  2:248492/1‑228 (MQ=255)
                                                                                                                                                                                       tAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGGTTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCtttt                                                >  1:142209/1‑250 (MQ=255)
                                                                                                                                                                                                     aGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGATGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTCACg                                  >  1:186024/1‑251 (MQ=255)
                                                                                                                                                                                                      ggCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGaa                                                                <  1:69660/220‑1 (MQ=255)
                                                                                                                                                                                                              cTATCGAAGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGCAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACATTAGTTATGAACGTTGCAACACCTTTATGTTTGGCGTTACGTTgcgcccc                         >  2:604517/1‑248 (MQ=255)
                                                                                                                                                                                                                   gACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACtt                    >  1:203930/1‑251 (MQ=255)
                                                                                                                                                                                                                   gACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACtt                    >  1:555380/1‑251 (MQ=255)
                                                                                                                                                                                                                     cGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTAGGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTa                  <  1:386152/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGAt             <  1:179502/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGAt             <  1:299296/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCGTc  >  1:581307/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
CAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCGTC  >  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: