Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,903,069 G→T 13.6% A392S (GCC→TCC)  ygcE → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,0690GT13.6% 199.5 / 4.4 82A392S (GCC→TCC) ygcEputative kinase
Reads supporting (aligned to +/- strand):  ref base G (22/47);  new base T (10/2);  total (33/49)
Fisher's exact test for biased strand distribution p-value = 1.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATGACGGTGCTGGACTGGCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTG  >  NC_000913/2902835‑2903302
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
tttatgaCGGGGTTGGAGTGGCTTACCAATCCGGGGGAACGCTAAAAACGGGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCAGGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCg                                                                                                                                                                                                                                          <  2:938077/235‑1 (MQ=255)
 tgatgaCGGTGCTGGACTGGCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAAAAGCGGACTGTTTATTCAGATTTTTGCCgcc                                                                                                                                                                                                                                        >  1:356776/1‑237 (MQ=255)
   atgaCGGTGCTGGACTGGCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGCCAATATGTGTAATGAATTGAATCCCTTTGCGCAGCATGTGATCATTACTGGCGGCGGTTCCAAAAACGATCTGTTTTATCAAATTTTTTCCGACgg                                                                                                                                                                                                                                     >  2:515741/1‑238 (MQ=255)
   atgaCGGTGCTGGACTGGCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACCGCGATCTGTTTATGCAAATTTTTGCCGACgg                                                                                                                                                                                                                                     >  1:1027364/1‑238 (MQ=255)
        ggTGCTGGACTGGCTGACCAATCCGGGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTc                                                                                                                                                                                                                                <  1:702406/238‑1 (MQ=255)
               gACTGGCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTTGAACAACGAACTATTTTTGGAAATTTTTTCCGCAGTGTTTAAACTTcc                                                                                                                                                                                                                        >  1:1011622/1‑239 (MQ=255)
                  tGGCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACACCGATCTGGTTATGCCGATTTTTTCCCACGGGTTCAACCTTCCgc                                                                                                                                                                                                                      >  1:921123/1‑237 (MQ=255)
                    gCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTTTTCCGAATTTTTCCGACGGGTTCAACCTTCCggacc                                                                                                                                                                                                                   >  2:93309/1‑236 (MQ=255)
                    gCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTTCCGACGTGTTCAACCTTCCGGCac                                                                                                                                                                                                                   >  1:988839/1‑239 (MQ=255)
                    gCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCCGTTCGAAAAGCGAGCTTGTTTAGCAAGTTTTTTGCGACGGTGTTAACCTTCCGGCac                                                                                                                                                                                                                   >  2:211137/1‑239 (MQ=255)
                    gCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGAGGTTCGACCAGCGATCTGTTTATGCGGATTTTTTCCGACGTGTTTAACCTTCCGGCac                                                                                                                                                                                                                   >  1:602592/1‑239 (MQ=255)
                    gCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATAGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACACTTACGACAATATGTGTAATGCAATGAATCACTTTGCGAAGCATGTCATCATTATGTGAGGCGGATTGAAAAGCACTCCGTTTTTTCCGATTTTTACCCCCCGGTTTAACCTTctgtgac                                                                                                                                                                                                                   >  2:512948/1‑233 (MQ=255)
                    gCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCACGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGCAGAACAATTCCGCCAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGTTCATTACTGCCGGCGGTTTGAACAGCGCCTTGTTTTTGAAGGTTTTTTCCGGCGGGTTTAAACCTCCCggac                                                                                                                                                                                                                    >  2:763517/1‑237 (MQ=255)
                     ctatccATTCCGTGGGAACCGTACAAACGGGGGATTATGATCGGTTTTGATTCCAGAATGGATTACCCATGGATATATCGTTCGAGACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCac                                                                                                                                                                                                                   <  2:837174/234‑1 (MQ=255)
                             tCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACg                                                                                                                                                                                                          <  1:263859/239‑1 (MQ=255)
                                  gggAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCAtt                                                                                                                                                                                                     <  1:952352/239‑1 (MQ=255)
                                     aaCCGTGAAAACGCGGGATTATGATCGGCTTTGATTCCAGCATTGATTAAGCATGGATATAGCGTTCGATACTGGAAAGCGTGGCGCTTACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTaa                                                                                                                                                                                                   <  1:31768/238‑1 (MQ=255)
                                       ccgttaaaaCGCGGGTTTATGGTCGGGTTTGATTCCAGCATGGATTCCGCATGGATATATCGGTCGATACTGGAAAGCGTTGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACgg                                                                                                                                                                                                <  1:619548/233‑1 (MQ=255)
                                       ccGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTCACGCCCTTAAAgg                                                                                                                                                                                                >  1:383872/1‑239 (MQ=255)
                                         gTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTTCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCCTTAACGGtt                                                                                                                                                                                              >  1:308262/1‑239 (MQ=255)
                                            aaaaCGCGGGATTATGGTCGGGTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTtgtg                                                                                                                                                                                           <  1:331177/238‑1 (MQ=255)
                                                        tATGATCGGCTTTGATTCCAGCATGGATTAAGCATGGATATATCGGTCGATACTGGAAAACGTGGCGCTGACGCTGAAGAACAATTACGAAAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTggggg                                                                                                                                                                               <  1:15357/239‑1 (MQ=255)
                                                          tggtCGGCTTGGATTCCCGCTGGGGTTACGGTTGGATATATCGTGCGAAACTGGAAAGCGGGGCGCTGACGCCGAAGAACAATTTCGGCAATATGTGTAAGGAAATGAATCCCTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCa                                                                                                                                                                             <  2:988839/236‑1 (MQ=255)
                                                          tGATCGGGTTTGATTCCAGCATGGATTACGGATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCa                                                                                                                                                                             <  1:577834/239‑1 (MQ=255)
                                                          tGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATAACGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCa                                                                                                                                                                             <  1:1079397/239‑1 (MQ=255)
                                                               gggtttGATGCCAGCTAGGATTACGCATGGATATATCGTGCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCTGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGAt                                                                                                                                                                        <  1:512948/236‑1 (MQ=255)
                                                               ggCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGGTTTTTCAGATTTTTTCCGCAGGGTTTAACCTTCCGGCACCACGTAACGCCATTAAAGGTTGTGCAAACCCgggggggggggg                                                                                                                                                                        >  2:1080201/1‑232 (MQ=255)
                                                                 gtttGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGAtt                                                                                                                                                                       <  1:217913/237‑1 (MQ=255)
                                                                           aGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGGTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGCGAGCGGTTAATAAAGCggg                                                                                                                                                            >  1:788065/1‑238 (MQ=255)
                                                                              cttgagtACGGATGGATTAGTCTTTCGATACTGGAAAGCGGTGGGCTGACGCTGAAGAACAATTACGCCAATATTGGTAATGAAAGGAATCACTTTGCGAAGCATGTGATCATTGCTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGt                                                                                                                                                            <  2:790335/230‑1 (MQ=255)
                                                                                  aTTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCt                                                                                                                                                      <  1:711761/238‑1 (MQ=255)
                                                                                           ggataatgtgTTCGATACGGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATAc                                                                                                                                            <  1:763517/230‑1 (MQ=255)
                                                                                           ggATATAGCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAAGACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATAc                                                                                                                                            <  1:361862/239‑1 (MQ=255)
                                                                                             atatatCGTGCGATACTGGAAAGCGTGGCGCTTACACTGAAGAACAATTACTACAATATGTGTAATGAAATGAATCCCTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCa                                                                                                                                                                                                                    <  2:820525/165‑1 (MQ=255)
                                                                                                     ttCGATACTGGAAAGCGTGGCGCTGACGGTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGc                                                                                                                                                                           <  2:683641/198‑1 (MQ=255)
                                                                                                     ttCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGc                                                                                                                                                                           >  1:683641/1‑198 (MQ=255)
                                                                                                       cGATACTGGAAAGCGTGGCGCTGCCGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGAATTTTGCCGGCGTGTTCAACCTTCCGGCACCACGGAAGGCCATTAACGGTTGTGCAAGCCTGGGGGCGGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGCTTACGc                                                                                                                                 >  1:901258/1‑238 (MQ=255)
                                                                                                       cGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACg                                                                                                                                                                                                          >  1:820525/1‑165 (MQ=255)
                                                                                                       cGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATAAAGCGGTAGGTCTGGGGCTATACCCGGa                                                                                                                                       >  1:713627/1‑232 (MQ=255)
                                                                                                       cGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGAAGCGCTTAATACACCGGTAGGCCTGGGGCATTAACCGGATTAccaa                                                                                                                                >  1:436550/1‑236 (MQ=255)
                                                                                                       cGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGGGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGGCTGGGGCTATACCCGGCTTACGc                                                                                                                                 >  1:782549/1‑238 (MQ=255)
                                                                                                       cGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGAACATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGa                                                                                                                                       <  2:713627/232‑1 (MQ=255)
                                                                                                           tcTGGAAAGGGTGGCGTTGACGGTGAATAACAATTACGACAATATGTGTAATTAAATGAATCCCTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTGGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACgg                                                                                                                            <  2:383872/238‑1 (MQ=255)
                                                                                                           gatGGAAAGCGTTGCGTTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTTTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACgg                                                                                                                            <  1:226833/237‑1 (MQ=255)
                                                                                                           aCTGGAAAGCGTTGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACgg                                                                                                                            <  1:1037295/239‑1 (MQ=255)
                                                                                                           aCTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACgg                                                                                                                            <  1:1016599/239‑1 (MQ=255)
                                                                                                                    cGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCCGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATa                                                                                                                   >  1:326901/1‑239 (MQ=255)
                                                                                                                    cGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTAAACCCGGATTAAGCAACGGccgcgcgta                                                                                                                   >  1:403817/1‑231 (MQ=255)
                                                                                                                    cGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGGCGGGGGCTATACCCGGGTTACGCAACGGCTGtgcgta                                                                                                                   >  1:350524/1‑234 (MQ=255)
                                                                                                                       ggCGCGGAGGCGTAAGAAAAATTGCCAAAATATGGGTAATTAAAGTAAACGCTTTGCGAAGCATGTGATCAGTACGGGCGGCGGGTCGAACAGCGATCTGTGTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACa                                                                                                                <  2:1067203/239‑1 (MQ=255)
                                                                                                                             tacgatgaAGAACAATTACGGCAATATGTGTAATGAAATGAATCACTTTGCGAAGCCTGGTAGCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGtt                                                                                                           <  1:814495/237‑1 (MQ=255)
                                                                                                                                gttgaAGAAAAATTTCGAAAATATGTGGAATGAAAGGAACTACTTTGCGAAGCGTGTGATGATTACCGGGGGCGGTTCGACCAGCGATTCGTTTAGGAAGATGTTTGCCGACGTGTTCAACCTGCCGGCACGACGTAACGCCATTAACGGGTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTcgcg                                                                                                       <  2:782549/237‑1 (MQ=255)
                                                                                                                                       aatcaTTTCGATATTATGTGTAATTAAAGGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTTTTTAGGCAGATTTTTGCCGGCGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                                                                                                <  2:18129/235‑1 (MQ=255)
                                                                                                                                       aaaaaTTAAGAAAATAAGTGTAATGAAAGAGATCACTTTGCGAAACAGTTGATCTTTACGGAGGGCGGTTCAAACAGGGATCGTTTTATGCGGATTTTTGGCGAAGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                                                                                                <  2:1027364/236‑1 (MQ=255)
                                                                                                                                       aaCAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATTTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                                                                                                <  1:1080201/239‑1 (MQ=255)
                                                                                                                                       aaCAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACCGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                                                                                                <  1:494656/239‑1 (MQ=255)
                                                                                                                                       aaCAATGCCGACAATATGTGTAATGAAATGAATCCCTTTGGGAAGGATGTGATCATTAATGGCGGCGGTTCGAACAGGGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                                                                                                <  1:124565/239‑1 (MQ=255)
                                                                                                                                        aacaTTACTACAATGTGGGTAATGAAATTAATACCTTTCGGAGCGCTGTGAGAATTTCTGGGGGGGGGTCGTACAGCGAATGTTTTTGTCAGTATTTTTCCGACGTGCTCAACCTCCCGGCACGAGCTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGTCTATAACCGGATTACGCAACGGCTGTCGATAACGTGGTTCGCGTGAAAGa                                                                                                <  2:389116/235‑1 (MQ=255)
                                                                                                                                               cGACAATAGGTGTAATGAAATGAATCCCTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTta                                                                                        <  1:383877/239‑1 (MQ=255)
                                                                                                                                                    atatGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGGGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGTTTACGCAACGGCTGTCGGTAACa                                                                                                                >  1:781480/1‑210 (MQ=255)
                                                                                                                                                    atatGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGGTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTTACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAAAATGGGTCGCGTGAAAAAAAACTTTAAAcag                                                                                   >  2:538476/1‑237 (MQ=255)
                                                                                                                                                          gtAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGGAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTtag                                                                             >  1:805065/1‑239 (MQ=255)
                                                                                                                                                           tAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATTCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAAACCCATTAACGGTTGTGCAAGCCTGGGGGCCGCGCTTAATACAGCGGTAGGGCTGGGGCTATACCCGGATTACGCAACGGCTGTCGGAAAAATGGGTGGCGGGAAAAAAAAATTTATAACGCTTaaga                                                                            >  2:340966/1‑239 (MQ=255)
                                                                                                                                                            tatgaaattaaTTCTTTTGGGAAGCATGTGATCATTACTGGCGGGGGTTCGAACGGGGTTCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGagag                                                                           <  2:326901/238‑1 (MQ=255)
                                                                                                                                                             atgaaatgaaTCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGGGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCGCGCTTAATACAGCGGTTGGTCTTGGGGTATACCCGGGTTACGGAAAGGGTGTCGCAAACATGGTTCGGGGGAAAAATGTATTTTTAccatttaggaga                                                                          >  2:672028/1‑229 (MQ=255)
                                                                                                                                                             atgaaatgaaTCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGGGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGGCTGGGGCTATACCCGGATTACGCCACGGCTGTCGTAAACATGGGTCGCGTGCAAAAAATCTTTATACCCATTggggaa                                                                          >  2:692762/1‑234 (MQ=255)
                                                                                                                                                              tgaaatgaaTTCCTTTTGGAATGATGTGATCATTTCTTGCGGCGGTTCGAACCGGGGACTGTTTAGGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGGAACGCCATTAACGGGTGGTCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCa                                                                         <  2:350524/239‑1 (MQ=255)
                                                                                                                                                                   tttaccccTTTGGGGAGAATGTGATCAGTTCGGGGGGGGTTTCGAACAGCGATTTGTTTAGGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGcc                                                                    <  2:921123/232‑1 (MQ=255)
                                                                                                                                                                   tgaaTCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCa                                                                                                                                                                             <  2:727186/134‑1 (MQ=255)
                                                                                                                                                                   tgaaTCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCa                                                                                                                                                                             >  1:727186/1‑134 (MQ=255)
                                                                                                                                                                             ttcGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGGTTAGGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGGCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTacg                                                          <  1:892562/237‑1 (MQ=255)
                                                                                                                                                                             tGCGAAGCAGGGGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTacg                                                          <  1:211137/239‑1 (MQ=255)
                                                                                                                                                                                  aGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                                                     <  1:93309/239‑1 (MQ=255)
                                                                                                                                                                                    cATGTGATCATTCCTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                                                                                                <  1:1051384/194‑1 (MQ=255)
                                                                                                                                                                                     catggtgtCATTTCGGGGGGCGGTTCGAACAGCGATCTTTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGa                                                    <  2:436550/231‑1 (MQ=255)
                                                                                                                                                                                      atggtaacATTACTGGCGGGCGGTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGGGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGa                                                    <  1:494711/230‑1 (MQ=255)
                                                                                                                                                                                                          ggTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCaa                                                                                                                               <  2:318913/142‑1 (MQ=255)
                                                                                                                                                                                                          ggTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCaa                                                                                                                               >  1:318913/1‑142 (MQ=255)
                                                                                                                                                                                                           ggtCGAACAGCGATCTGTTTATGCAGATTTTTGTCGACGTGTTCAACCTTCCGGCCCGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGa                              <  1:900439/236‑1 (MQ=255)
                                                                                                                                                                                                            tccgatccgCGATCTGTTTATTAAGATTTTTTGTGAGGGTGTCAACCTTCAGGAACGACGAAACGCCAGTAACGGTTGGCCAAGCCGGGGGGCAGGTATTAATACCGCGGTAGGGCTGGGGCTATAACGTGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGAc                             <  2:730757/230‑1 (MQ=255)
                                                                                                                                                                                                              cGAACAGGCATCTGGTTAAGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCTa                          <  1:340966/239‑1 (MQ=255)
                                                                                                                                                                                                                                        tttttGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGTCGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAAAATACTCATGTGATCCTg  >  1:1067696/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
CTGATGACGGTGCTGGACTGGCTGACCAATCCGTGGGAACCGTACAAACGCGGGATTATGATCGGCTTTGATTCCAGCATGGATTACGCATGGATATATCGTTCGATACTGGAAAGCGTGGCGCTGACGCTGAAGAACAATTACGACAATATGTGTAATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTG  >  NC_000913/2902835‑2903302

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