Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,131 G→A 20.3% A337T (GCC→ACC)  icd → isocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,1310GA20.3% 86.9 / 21.5 51A337T (GCC→ACC) icdisocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
Reads supporting (aligned to +/- strand):  ref base G (25/15);  new base A (4/6);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACG  >  NC_000913/1195903‑1196353
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
ggCCCGTGGCTGAAAGTTAAAAACCCGAACCCTGGCAAAGGGATCGGCATTAAAGCCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAg                                                                                                                                                                                                                                <  1:325064/229‑1 (MQ=255)
ggCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGGTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAg                                                                                                                                                                                                                                <  1:478260/229‑1 (MQ=255)
  cccGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAg                                                                                                                                                                                                                                >  1:171300/1‑227 (MQ=255)
  cccGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAg                                                                                                                                                                                                                                <  2:171300/227‑1 (MQ=255)
   ccGTGGCTGAAAGTTAAAACCCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAAcc                                                                                                                                                                                                                              <  2:70402/228‑1 (MQ=255)
             aaGTTAAAACCCCGAACAGTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTa                                                                                                                                                                                                                    <  1:626002/228‑1 (MQ=255)
                   aaaaCCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGCATCGGTATCCCCCCTGGTGCAAACATCGGTGACGAAGGCGCCCTGTTTGAAGCACCCCCCGGTACACCGcc                                                                                                                                                                                                             >  2:338850/1‑229 (MQ=255)
                                gggcaaagggtccttgaTTAAAGAGGTGATTGGCGTTGCATTCCTGGAACAGAGCCGGCTGCGTCCGGCTGAATATGATGTTATCGCCCGTATGCACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTTGTGCAAACATCGGGTACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGt                                                                                                                                                                                                <  2:109251/213‑1 (MQ=255)
                                  gCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATCTGATGTTATCGCCTGTATGAACCTGAACGGTGACTCCATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGGGACCAATAGGCCCTGTTTGAGGCCACCCACGGAACTgcgcagcaccaggccggcag                                                                                                                                                                                              >  1:595630/1‑213 (MQ=255)
                                          aTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGGGACGAATGCGCCCGTTGTGAAGCCACCCACGGTACTGCGCCGCAACATGCCGGTCGGGCAAAg                                                                                                                                                                                      >  1:599777/1‑228 (MQ=255)
                                                  tAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCt                                                                                                                                                                              >  1:695097/1‑229 (MQ=255)
                                                      ggcGTGAATGCTGATGCATTCCTGCAACAGATCCTGCTGCGGCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGct                                                                                                                                                                          <  1:228646/227‑1 (MQ=255)
                                                         gTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTTACGAATGCGCCCTGTTTGAAGCCACCCCCGGTACTTCCCCGAAATATGCCGGTCCGGCAAAAGGAAAACCAGGcccg                                                                                                                                                                        >  2:302450/1‑225 (MQ=255)
                                                         gTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGGCAAAGTAAAACCTGGctct                                                                                                                                                                        >  1:673498/1‑228 (MQ=255)
                                                         gTGATTACTCCTGCATTCCTGCAACAGATCCTGCTGCCTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCATCGCCGGTTGGCGGTATCGGTATTGACCCTGGTGCAAACCACGGGGACGAATGCGGCCTGTTTGAAGCACCCCACGGTACTGGCCCGACAAATGCCGGTCAGGAACAAGAAAAACCGGGCGCta                                                                                                                                                                       >  2:234804/1‑229 (MQ=255)
                                                          tGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCCACCATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGCCAACGTAAATCCTGCCTCta                                                                                                                                                                       >  1:306348/1‑228 (MQ=255)
                                                                  ggtGCTTTCCTGCAACAGATCCTGCCGCGTCCGGCTGAAGATGATGTTGTCGCCTATATGAACCTGAACGGTGACTTCATTTCTGACGACCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAATATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtct                                                                                                                                                                <  2:467248/225‑1 (MQ=255)
                                                                  gaggCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtctct                                                                                                                                                              <  1:486997/226‑1 (MQ=255)
                                                                   aTGCATTCCTGCCAAAGATCCTGCTGCGTCCGGTTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtctctc                                                                                                                                                             <  2:306348/229‑1 (MQ=255)
                                                                      ccttCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATTGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                          <  1:20147/227‑1 (MQ=255)
                                                                      cATTCCTGCAACAGTTCCTGCTGCGTCCGGCTGAATATGATGGTATCGCCTGCATGAACCTGAACGGTGGCTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                          <  1:325596/229‑1 (MQ=255)
                                                                                               tCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCCCCGCTGAGCTGATGCTGCGCCACATggtt                                                                                                                                 >  1:491627/1‑227 (MQ=255)
                                                                                               tCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCACCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAACCCACCCACGGTACTGCGCCGAAATATGCCGGGCCGGCAAAAGTACTACCTGGCTTTCTTCTTCTCACCGCTGAAGTGATGGTGGGCCAAATGGGt                                                                                                                                 >  2:693355/1‑229 (MQ=255)
                                                                                                  ggCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTCCTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTTTCCGCTTAGATGGTGCTGCGCCACCTGGGTTgt                                                                                                                              >  2:40540/1‑228 (MQ=255)
                                                                                                  ggCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTg                                                                                                                               >  1:173239/1‑228 (MQ=255)
                                                                                                              tgTTATCGCCTGTATGAACCTGAACGGAGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGGTTGAAGCCACCCACGGTACTTCGCCGAAATATGCCGGTCAGGAGAACGTAAATCCTGGCTCAATTATTCTCTCCGCGGGGATGATGCTGCGACACAGGGGTTGGACCGAAGCgccc                                                                                                                  >  1:747388/1‑226 (MQ=255)
                                                                                                                  aTCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                                                                              >  1:477542/1‑229 (MQ=255)
                                                                                                                      cttttttGAACCTGACAGGTTACTACTTTTCGTGCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAACCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTaa                                                                                                           <  2:327795/223‑1 (MQ=255)
                                                                                                                          tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTTACGAATGCGCCCTGTTTGAAACCACCCACGGCACTGCGCCGAAATAAGCCTGTCCGGACAAAGTAAATCCTGGCTCTATTATTCACTCCGCTTCGATGAGTCTGCGCCACCTGGGTTGGACCggaccagccgcacttattttt                                                                                                      >  2:351727/1‑209 (MQ=255)
                                                                                                                          tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAAAATGCCGGTCAGGACAACGAAAATCCTGGCCCTATTATTCTCTCCGCTGGGATGATGTTCAGCCAAATGGGGTGGACCGAACGGGCTGACTTAATTGtt                                                                                                      >  2:250968/1‑229 (MQ=255)
                                                                                                                          tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCATGTTTGAAGCCCCCCGCGGGACTGCGCCGAATTATGCCGGTCAGGACAACGTAAATCCTGGCTCTGTTCTTCTCTCCGCTGAGATGATGCTGCACCACATGGGTTGGGCGCAATAGCCTGCCTTAATtttt                                                                                                      >  2:606439/1‑229 (MQ=255)
                                                                                                                          tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGGATCGCCCCTGGTGCAAACATCGGGGAAGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATTATGCGCCACATGGGTTGTACCGAAGCGGCTGACTTAATTGtt                                                                                                      >  2:673414/1‑229 (MQ=255)
                                                                                                                           aTGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                     <  2:477542/229‑1 (MQ=255)
                                                                                                                                ccTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCCATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGCAGCGGCTGCCTTAATTGTTTAAGGt                                                                                                >  1:698624/1‑229 (MQ=255)
                                                                                                                                    aaCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCAGCTGGATTACTTGTTACAggactgg                                                                                            >  1:299598/1‑224 (MQ=255)
                                                                                                                                    aaCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGGGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGGGATGATGCTGCGCCACATGGGTTGGACGGAAGCGGCTGACTTAATTGTTAAAggaaggg                                                                                            >  1:250237/1‑224 (MQ=255)
                                                                                                                                               cATTTCTGACGCCCTGGCAGCTCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc                                                                                 <  1:302450/229‑1 (MQ=255)
                                                                                                                                                   tCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAAGGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCCCTATTATTCGCTCCCCTGAAATGATGCTGCGCCACATGGGTTGGACCGAGGCGGCTGCCTTAATAGTTAAGGGGATGGAGGGCGCAATAAAcg                                                                             >  1:344937/1‑229 (MQ=255)
                                                                                                                                                      gACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGa                                                                          <  1:184958/229‑1 (MQ=255)
                                                                                                                                                       gcGCCCTGGCCGCGCAGGTTGGCGGTTTCGGTATCGCCCCTGGTGCAAACATCGGGGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACCGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCGCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGa                                                                          <  2:698623/227‑1 (MQ=255)
                                                                                                                                                       aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAAcgcg                                                                           >  1:509026/1‑227 (MQ=255)
                                                                                                                                                                  gcgcGGGTTGGGGGTTTCGGTATCGCCCCTGGTGCAAACCTCGGTGACGAATGCGCCCTGGTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                                                              <  2:305732/229‑1 (MQ=255)
                                                                                                                                                                   cgcAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGGCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATCCTGCGCCACAAGGGTTGGACCGAAGCGGCTGACCGAATTGTTAAAGGGATGAAAGGGGCAACCAAAGGCAAAACCTTAAccat                                                             >  2:352040/1‑227 (MQ=255)
                                                                                                                                                                    gcAGTTTGGCTGTATCGGTATCGCCCCTGGTGCAAACAACGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCGCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                            <  2:684346/229‑1 (MQ=255)
                                                                                                                                                                              ggtaacggtatcGTCCCTGGTGCAAACATCGGTGACGAATGCGACCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTACCCTATGACTTCGa                                                    <  1:606440/227‑1 (MQ=255)
                                                                                                                                                                                         tccccccTGGTGCAAACAGCGGTGAAGAATGCGCCCTGTTTGAAGCCACCCACGGTAATGCCCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                        <  2:173239/226‑1 (MQ=255)
                                                                                                                                                                                         cgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                        <  1:673415/228‑1 (MQ=255)
                                                                                                                                                                                          gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTAGAGTGTCTgctgg                                         >  1:220528/1‑226 (MQ=255)
                                                                                                                                                                                            ccccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTGTTAAAGGTATGGAAGGGGCAACCACCCCGGAAACCGTAACCCATAACTTagaggtgcggagggtggg                                     >  2:351459/1‑215 (MQ=255)
                                                                                                                                                                                                   gtcAAACCTCGGTGACGAATGCGCCCTGGTGTAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTaaa                              <  1:338850/227‑1 (MQ=255)
                                                                                                                                                                                                       aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATCCTGCGCAACATGTGTTCGACCGAAGCGTCTGCCTTAATTGTTAAAGGTATGAAAGGCGCAATCAAAGCGAAAACCATAAACAATGACTTAGCGAGGCTgatggcaacggccaaaaggc                          >  2:423007/1‑214 (MQ=255)
                                                                                                                                                                                                            cGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCCATCAACGCGAAAACCGTAAACTATGACTTCGAGTGTCGGATGGTGGGCGCTAAAACGCGGaaa                     >  2:636194/1‑229 (MQ=255)
                                                                                                                                                                                                             ggTGACGAATGCGCCCTGTTTGAAACCACCCGCGGTACTGCGCCGAAATATGCCGGTCCGGACAACGTAACACCTGGCTCTATTCTTCTCTCCGCTGAGCTGATGCTGCGCCACATGGGTTGGACCGCACCGGCTGACTTAGTTGTTAAAGTGATGGCAGGGGCAATCAACGCGAAAACCACAACCTATGACTCCGAGCGTGTTAGGGATGGTGCTAAACAGCTGaaaa                    >  2:497512/1‑228 (MQ=255)
                                                                                                                                                                                                                    aaTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGCCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGACATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAAACCCTGAAATGTTCaga             >  1:401861/1‑229 (MQ=255)
                                                                                                                                                                                                                               tttGAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg  <  2:446582/229‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
GGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACG  >  NC_000913/1195903‑1196353

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: