Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,701,625 A→G 14.9% G223G (GGA→GGG malY → PLP‑dependent beta‑cystathionase and maltose regulon regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,701,6250AG14.9% 153.2 / 2.3 70G223G (GGA→GGGmalYPLP‑dependent beta‑cystathionase and maltose regulon regulator
Reads supporting (aligned to +/- strand):  ref base A (20/37);  new base G (8/4);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 5.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCATCTATCTGAAAGATAACCTGACGTATA  >  NC_000913/1701390‑1701860
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
aTGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGgaga                                                                                                                                                                                                                                               <  2:664598/238‑1 (MQ=255)
     tggttTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTggg                                                                                                                                                                                                                                         <  1:680740/238‑1 (MQ=255)
      ggttTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCTAGCACCCGCATATTCCCTGGAGTAAGGGGGCTCGCGGGGAcccggtg                                                                                                                                                                                                                                        >  2:229295/1‑233 (MQ=255)
      ggttTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAAATCATGGCTGACCTGTGCGAGCGTCATGGGGTGCGGGTATTTTCCGTGGAAATCCATATGGGAATGGTTTGGGGCGAGCAGCCCCCTATTCCCTGGGGCAATGGGGGCCCCGGGGAATggggg                                                                                                                                                                                                                                       >  2:265234/1‑237 (MQ=255)
                tatGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATGATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                                                                                                                                                                                                             <  1:277036/239‑1 (MQ=255)
                  tGGGCAAGTTGGAAGCCGGTTTGGGGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCg                                                                                                                                                                                                                           <  2:480511/239‑1 (MQ=255)
                          ttGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGGGGAGACTGGGCGTTGCTAACGTCGGGCTCGaa                                                                                                                                                                                                                   >  1:590841/1‑239 (MQ=255)
                          ttGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGGGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGGGACTGGGCGTTGTTAAAGTCGGGGTCGaa                                                                                                                                                                                                                   >  1:1045053/1‑239 (MQ=255)
                          ttGAAAGGCGTTTTGGGGAACCCAGATGTTAAAATTATGCCCCTTGGGAGCCCACAGAATCCTGCCGGGAAAGGTGGGACGTGCGATGAGCTGGAGATCATGGCTGGCCTGTGCGAGCGTCATGGTGTGCGGGTTATTTTCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGaa                                                                                                                                                                                                                   <  2:193186/239‑1 (MQ=255)
                           tGGAAGCCGTGTTGGCAAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCTGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGaa                                                                                                                                                                                                                   <  1:820095/238‑1 (MQ=255)
                                   gtgtTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGGTGCTAACGTCGGGCTCGAAAAGTTTCaa                                                                                                                                                                                                        >  1:633310/1‑239 (MQ=255)
                                    tgtTGGCGAAACCAGAGTGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCCTATTCCCTGGAGTAATGTGGTTCGCGGGACATGGGGGTTGCTAACGCCGGGCCCGAAAAGTTTCaaa                                                                                                                                                                                                       >  1:798454/1‑238 (MQ=255)
                                                    agggaaaaTTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg                                                                                                                                                                                       <  1:195524/235‑1 (MQ=255)
                                                    aTGTAAAATTATGCCCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGGGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATAGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCgg                                                                                                                                                                                      <  1:640340/239‑1 (MQ=255)
                                                      gttaaaTTTTTCTCCTGTGGTGCCCACAGAAGCCGACCGGGAAGGTGTGGACGTTCGATGAGCTGGAGATCATGGCTGACCTGTGGGAGCGTCATGGTGTGGGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGt                                                                                                                                                                                     <  2:590841/236‑1 (MQ=255)
                                                           ttttttcTCCTGGTTGGCCCCAAGAAGCCTACGGGGAAAGGTTGGAGGGTGGATGGGCTGGAGTTCATGGCTGACCTGGGCGGGGGTTATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTa                                                                                                                                                                                <  2:1072670/233‑1 (MQ=255)
                                                           aTTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGGGATGAGCTGGAGATCATGGCTGACCTGTGCGGGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTa                                                                                                                                                                                <  1:630805/239‑1 (MQ=255)
                                                           aTTATGCGCCTGTGTAGCCCACAGAATCCTAACGGGAAAGTGGGGACGTTCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTa                                                                                                                                                                                <  1:1018810/239‑1 (MQ=255)
                                                           aTTAGGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCCTGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGTCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTa                                                                                                                                                                                <  1:709727/239‑1 (MQ=255)
                                                             tATGCTCCTGTGTAGCCCACAGAATCCTTCCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGGGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGGGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCaa                                                                                                                                                                                                        <  1:837618/213‑1 (MQ=255)
                                                               tGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATTTGGCTCGCGGGGACTGGGGGTTGCGAACGTCGGGCTCGAAAAGTTTCAAATATTCCCCCCTGACCGGGGCTTACggg                                                                                                                                                                            >  1:41135/1‑239 (MQ=255)
                                                                gCTCCTGTGTAGCCCAAAGAATCCTACCGGGAAAGGGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGGGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGa                                                                                                                                                                           <  1:1011401/239‑1 (MQ=255)
                                                                    cTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGGGCTTACGGGAtttt                                                                                                                                                                       >  1:206230/1‑237 (MQ=255)
                                                                            ccccaaGACTCCTCCCGGGAAGGTGTGTGCTGGGGATGGGCGGGAAGTCATGGTTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAAt                                                                                                                                                                <  2:1045053/233‑1 (MQ=255)
                                                                              cacaGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGGCGGGCTCGAAAAGGTTTAATATTCCCGCCCTGACCGGGGCTTACGGGGTTTTTGAAAATaac                                                                                                                                                             >  1:280612/1‑237 (MQ=255)
                                                                                caGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGGGCTTACGGGATTATAGAAAAAAGCat                                                                                                                                                           >  1:893545/1‑238 (MQ=255)
                                                                                 aGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGCGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGAGTATAGAAAATAGCAGt                                                                                                                                                          >  1:886076/1‑239 (MQ=255)
                                                                                  cattCCGACCGGGACAGTGGGGGAGTGGGATGGTGTGGAGAGCATGGCTCACCTGGGCGAGCGGCATGGTGTGCGGGTTATGTCGGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATGCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGTGTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTa                                                                                                                                                         <  2:1005739/236‑1 (MQ=255)
                                                                                   aaTCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCTTATTCCCTGGAGTAATGTTGCTCTCGTGGACTCGGCTGTGCTAACCCCGGGCTCGCAAAGTTTCCAATATCCCCCCCTCCCCGGTGCTTCCGGGGTTTAAGaaaaaaaccaca                                                                                                                                                         >  2:329745/1‑230 (MQ=255)
                                                                                   aaTCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACAGGGCGTTGCTACCGTCGGGCTCGCAAAGTTTCAATATTTCCGCCCTTCCCGGGGTTTCCGGGAGTATAGAAAATAGcgctt                                                                                                                                                         >  1:687017/1‑234 (MQ=255)
                                                                                      ccTTCCGGGAAAGTGGGGAGGTGGGATGGGCTGGAGATCATGGCTGACCTGTGCGGGGGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTGGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCcg                                                                                                                                                     <  2:280612/239‑1 (MQ=255)
                                                                                        tACCGGGAAAGGGGGGACGTGCGATGAGCGGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTGTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCcgcg                                                                                                                                                   <  1:759065/239‑1 (MQ=255)
                                                                                         tccGGGAAATTTTGGACGTGTGTTGGGCTGGAGATTGTGGCTGCCCTGTGCGAGTGGCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGa                                                                                                                                                  <  2:348628/238‑1 (MQ=255)
                                                                                          ccGGGAAAGTGTTGGGGGGCGATGAGCTGGAGATCATGGGTTACCTGTGCGAGGGTCATGGTGTGGGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGAt                                                                                                                                                 <  2:893545/239‑1 (MQ=255)
                                                                                                             cGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGGGACTGGGCGTTGCTCACGTCGGGCCCGCAAAGGTTCACAATTTCCGCCCCGCCCGGGGCTTCCGGGGTTATAGAAAAAAGCAGAAGCCGCCGCGCCGACTTATCCGCAcctg                                                                                                                              >  2:60505/1‑236 (MQ=255)
                                                                                                             cGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGAc                                                                                                                                                                                         <  1:988684/180‑1 (MQ=255)
                                                                                                             cGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGCCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGAc                                                                                                                                                                                         >  2:988684/1‑180 (MQ=255)
                                                                                                                tGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTGGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATtata                                                                                                                                                                      <  1:340599/196‑1 (MQ=255)
                                                                                                                      ggggggtgtGGCTGTCCTGGGGGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGAGATGGTTTGGGGCGAGCAGCCGAATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTg                                                                                                                     <  2:1053715/231‑1 (MQ=255)
                                                                                                                            ccgtgggtcccTGGGCGGGTGTCATGGTTGGGGGGTTCTTTCCAATGAAATCACTGTGGGTAGGGGTTGGGTGGAGCAGCCGCATATTCCCTGGAGTAAGGTGGCTCGGGGGGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATAGTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATggg                                                                                                                <  2:614768/229‑1 (MQ=255)
                                                                                                                              tGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTTTAGAAAATACCAGTAGCCCCGGTGCCTATTTTTCGGCCCTGAAAGGCCGTGATTGGCtt                                                                                                            >  1:956913/1‑239 (MQ=255)
                                                                                                                                    aCCTGTGCGATCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGGTTTTAGAAAATACACATAACCACCAGGCCCATTTTTAGGGACTGAAAAGCAGAGATGGGCGtttttc                                                                                                       >  2:937489/1‑238 (MQ=255)
                                                                                                                                     ccTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGt                                                                                                                                                                                     <  1:161829/160‑1 (MQ=255)
                                                                                                                                         tgCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGGGACTGGGCGATGCCAACCTCGGGGTCGGAAAGTTTCAATATTCCCGCCCTGACAGGTACGTACGGGGTTCTAGAAAAAAAAAGTAACCCCGATGCCCTTTTTTCGGCACGTAAAAGCCGTGGTGGGAGTTTTTTCCCCtc                                                                                                  >  2:1053893/1‑238 (MQ=255)
                                                                                                                                         tgCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAACCCCGATGCCTATTTTTCGGCACTTAAAAGCCGTGATGGGCTTTTTTCCCCTTc                                                                                                 >  1:374486/1‑239 (MQ=255)
                                                                                                                                                        gtgGGGGTTTTTTCCATGAAAATCCTTATGGTTTTGGGTTGGGGGGAGCATGCGAATATTCCCTTGAGTAATGTGGCGCGCGTGGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTaa                                                                                   <  2:798454/238‑1 (MQ=255)
                                                                                                                                                        gtgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCTAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTaa                                                                                   <  1:158211/238‑1 (MQ=255)
                                                                                                                                                        gggggggTTAATTCCGTTGAAACCCATATGGATATGGTTTGGGGCGAGCAGCCGGATATTCCCGGGAGTAATGTGGCTTGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTaa                                                                                   <  1:955181/234‑1 (MQ=255)
                                                                                                                                                           cGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACCGGGCGTTGCTAACGTCCGGCTCGGAAAGGTTCCATTTTCCCCCCCTGCCCCGCGCTTACGGGGTTATAGAAAATAGCCGTCTCCGCGGCTCCCATTTTTCGGCCCTGAAAGGCCGTGGTGGGCTTTCTTCCCCTTCGGGACCGGGGGTAAATg                                                                                >  2:667627/1‑238 (MQ=255)
                                                                                                                                                                                 cATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGaa                                                                                                                             <  1:895444/173‑1 (MQ=255)
                                                                                                                                                                                         tatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGGTTATAGAAAATAGCAGTAGCCGCGATGCCTATTTTTCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGGTAACTGCCCATAACGCCGCCCATCAGCAAggggggc                                                  >  2:1069245/1‑234 (MQ=255)
                                                                                                                                                                                                  gggCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGGGACTGGGCGTTGCTTACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGGGCTTACTGGGTTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGGACTGAAAGGGCGTGATGGGCTTTCTTCCCCTTCGGGACTGGCGTTAACTGCCCACACTGCCCCCTGCCAGCAAGGAGGCCCGGGGCTgg                                         >  2:175416/1‑239 (MQ=255)
                                                                                                                                                                                                   ggCGAGCAGGCGAATATTCACTGGAGTAATGTGGTTCGCGGAGACTGGGCGTTGCTAACGGCGGGCGCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTgg                                         <  1:883585/238‑1 (MQ=255)
                                                                                                                                                                                                      gAGCAGGCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACGGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATGGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATg                                      <  1:1083506/238‑1 (MQ=255)
                                                                                                                                                                                                       agcagcCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGGTAACCTCGGGCTCGAAAAGGTTCAATATTACCGCCCCGCCCGGGGCTTACGGGATTATAGAAAATAACAGTAGCCCCGCTGCCTATTTCTCGGCCCCGCAAATCACAGATGGGCTTTCTTACCCCTCCGTACCGGCGTTTACTGACCCAATTCCCGCCCATCACCAAAGGGCGCCGCGGGTGGCTg                                      >  2:478993/1‑237 (MQ=255)
                                                                                                                                                                                                       agcagcCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCgg                                                                                                                                    <  1:1049868/144‑1 (MQ=255)
                                                                                                                                                                                                       agcagcCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCgg                                                                                                                                    >  2:1049868/1‑144 (MQ=255)
                                                                                                                                                                                                       agcagcCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTTGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATg                                      >  1:892638/1‑237 (MQ=255)
                                                                                                                                                                                                               aatatTCCCTGGAGTAAGGTGGCTCGCGGAGACTGGGCTTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCat                            <  1:1069245/238‑1 (MQ=255)
                                                                                                                                                                                                                atatTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTTGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCat                            <  1:754170/238‑1 (MQ=255)
                                                                                                                                                                                                                 tatTCCCTGGCGTAATGTGGATCGCGGAGACTGGGCGTTGCTAACGTAGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCatct                          <  1:373356/239‑1 (MQ=255)
                                                                                                                                                                                                                 tatCCCCTGGATTAAGGTGGCTGGGGGAGACTGGGCGTTGCTAACGTCGGTGTCGAAAAGTTTCAATATTCCCGCCCTTACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCatct                          <  2:956913/239‑1 (MQ=255)
                                                                                                                                                                                                                      ccTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGGTGCCTTACGCATCTATCTg                     >  1:296455/1‑239 (MQ=255)
                                                                                                                                                                                                                             ttaTGTGGCTCGCGGGGACTGGGCGTTGCTAACGTCGGGCCCGAAAATTTTTAATATCCCCGCCCTGACCGGTGCTTACGGGATTAAAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCATCTATCTGAAAGAt               <  2:363171/236‑1 (MQ=255)
                                                                                                                                                                                                                             tAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGGGCTTACGGGATTATAAAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCGTTGTTCCCCTTCGGTACTGTCGTTAACTGCCCATAATCCCCCCTTTCAGCAAAGAGGCCCGTGGCTGGATGCCTTACCCACCTATCTCAAAGaaa              >  2:21242/1‑237 (MQ=255)
                                                                                                                                                                                                                               aggtgGCTCCAGGAGACTGGGCGTTGCAAACGTCGGGCTCGAAAAGTTTCAATGTTCCCGCCCTGACCGGGGCTTACGTGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCATCTATCTGAAAGATAAc            <  1:841506/237‑1 (MQ=255)
                                                                                                                                                                                                                               aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCATCTATCTGAAAGATAAc            <  2:892638/239‑1 (MQ=255)
                                                                                                                                                                                                                               aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAAGGCCGCCGTGGCTGGATGCCTTACGCCTTTATCTGaaaagtaa             >  1:278131/1‑233 (MQ=255)
                                                                                                                                                                                                                               aTGTGGCGCGCGGAGACGGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCATCTATCTGAAAGATAAc            <  1:667627/239‑1 (MQ=255)
                                                                                                                                                                                                                                  gggCTCGCGGAGCCTGGGCGTTGCTAACGTTGGGCTCGAAAAGTTTCAATATTCCCGGCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCATCTATCTGAAAGATAACCTg         <  1:329745/238‑1 (MQ=255)
                                                                                                                                                                                                                                         cgGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCCCCCTGACCGGGGCTTACGGGGTTTTAGAAAAAAGCAGTAGCCGCGATGCCTATTTATCGGCCCTGAAAGGCCGTGAGTGGCTTTCTTCCCCCTCGGTACTGGCGTTAAATGCCCCTATCGCCGCCCTTCCACCAGGCGCGCCCGGGCTGGGTGCCTTCCGCCTCCATCTGAAAGaaaaaccaccgaaaa  >  2:1074363/1‑228 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
ATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTGCTTACGGGATTATAGAAAATAGCAGTAGCCGCGATGCCTATTTATCGGCACTGAAAGGCCGTGATGGGCTTTCTTCCCCTTCGGTACTGGCGTTAACTGCCCATATCGCCGCCTATCAGCAAGGCGCGCCGTGGCTGGATGCCTTACGCATCTATCTGAAAGATAACCTGACGTATA  >  NC_000913/1701390‑1701860

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