Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,298 C→T 13.5% L58F (CTT→TTT)  yfeN → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,2980CT13.5% 141.0 / 2.4 64L58F (CTT→TTT) yfeNputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base T (6/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCACGACGATATGTTCCTG  >  NC_000913/2525093‑2525525
                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
tCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTTGAGATATCGGCGTGGATTATCTGGGCTGGACCAGCCGGACTACCGAAAAAACGGCGTCCAAAACACCAAAAGAAGATTTTGGGTAACATTAAATTTAAGgggggcccaac                                                                                                                                                                                                                  >  2:143884/1‑222 (MQ=255)
tCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATTTGGGCTGGACCAGGCGGACTAACGGAAAATCGTCGACCAAATCACACAAAGGAGCTTTTGGTTCTCTTGAATTTGAAGGTGgccgttaa                                                                                                                                                                                                                  >  1:292069/1‑223 (MQ=255)
tCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATTTGGGCTGGACCAGCCGGGCTTACGGAAAAACGTCACCCAAATAACAACAAAAAAGTTTTGGGTTATTTTAAATTTGAgggggggccaaa                                                                                                                                                                                                                   >  2:672386/1‑223 (MQ=255)
tCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGGAAAACCGTCGACCAAATTACCCCAAGAAGATATTGGGTTTCTTGAATTTAAAGggggggtcaac                                                                                                                                                                                                                  >  1:113253/1‑222 (MQ=255)
tCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATATGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGGTTTTGGTTATCTTGAATTTTAAGGTGGCGCCAAc                                                                                                                                                                                                                  >  1:625892/1‑229 (MQ=255)
tCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCAATATCTGGACTCGACCCGCCGGAGTTCCGGAAAATCGACGACCCAATTACACAAAGATGAATTTGGTATTTTTGAATTTTAAGGTgtgcgcaac                                                                                                                                                                                                                  >  1:737044/1‑221 (MQ=255)
  tactaATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTCTCTGGACTGGACCCTCCGGACTACCTAAACATCGTCGACCCAATCCCAAAAAAGTGCTTTTGGGTTTCCTGAATTTGAAGGAGGGCCTAAAtt                                                                                                                                                                                                                >  2:476045/1‑229 (MQ=255)
       atatCGATCGGCATAGGGAATTATTCTATGAAAAAACTTCTTTTAACTCTGACACTTTCCTCTATAGTAGCGATACCGGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGtt                                                                                                                                                                                                           <  2:2998/229‑1 (MQ=255)
            gATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTTCAGATATCGGCGTGCATTTTCTGGAATGGACCAGCCGGACTACCGCAAAAACGGCGACCCAATCAAACCAAGAAGAGTTTGGTTTACTTGGATTTTAATGGGGGCGCAACATTTGTTTGGga                                                                                                                                                                                                      >  1:708508/1‑229 (MQ=255)
                ggCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAGCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGaaa                                                                                                                                                                                                  <  1:672387/229‑1 (MQ=255)
                      gggAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGCTTTTCAGATATCGGCGTGCACTTATTTGACTTGAGCAGCCGGATTAACGAAAAATTGTCGAACAAATCAACACAAAGTAGTTTTGGGTGACTTGAATTTGAAGTAGGCGCCAACTTTACATTGGGAGGAAttgaa                                                                                                                                                                                             >  1:582051/1‑224 (MQ=255)
                      gggAATTATGCTATAAAAAAACTTCTTTTAACTCTGACATTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAAGTTAAAGGCGGTTTTGCAGATAACGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACg                                                                                                                                                                                            <  2:215896/229‑1 (MQ=255)
                            tATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATTCAGAATTTAAAGGCGGTTTTGCAGAAATCGGCGTGCATTATCTTGACTTGACCAACCCTGCATCCGAAAAAACTGCTGACAAATACCCCAAAGATGATTTTGGTTATATTTAATTTTGAAGGTGCGCTAACTTTAATTTGGGGGAAATGTAAggggttt                                                                                                                                                                                      >  2:532701/1‑225 (MQ=255)
                            tATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTTCAGAAATCGGCGGGCCTTATTTTGCCATGAACCGCCCGCGTACCCGAAAATAGTCCGCCAAATCACACAAAAAGTATTTTGGTTATCTTGATATTGTAATGGGCGCTAAACTTAGATTGGGGGAAAtatcaagggttt                                                                                                                                                                                      >  2:454824/1‑218 (MQ=255)
                            tATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGGCTGGACCCGCCCGACTTCCGGAAAATCGTTCACCCAATCAAACAAAGAAGATTTTTGTTATCTTGAATTTGAAGGTGGCCCCAAACTTAATTTGGGGTAAATGTTCGGGtttt                                                                                                                                                                                      >  1:236367/1‑227 (MQ=255)
                            tATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATCTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCACCCGGCCTACCGAAAAATCGTCGACCCAATCACACAAAGATGATTTTGGGTATCTTGAATTTTAAGGTGGCGCTCACTTTAGTTTGGGAGAAATggtcggtgtct                                                                                                                                                                                      >  1:640963/1‑219 (MQ=255)
                                cttttaaaaaaCATTTTTTAACTCTGACACTTTCCTCTATATTAGCGATTCCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATAGCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGa                                                                                                                                                                                  <  2:625891/224‑1 (MQ=255)
                                       aaatctttttttAACCTTGAAACTTTCCTCTATTTTGGCGATTCCTTGTGTGTCTTCTGCAGAATTTAAAGGCGGTTTTGACAATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGaa                                                                                                                                                                          <  2:113253/221‑1 (MQ=255)
                                       aaatCATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATTTCTGTTGTATCTCATGCAGAATTTAAAGGCTGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGCGCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGaa                                                                                                                                                                          <  2:369570/229‑1 (MQ=255)
                                       aaaatctCTTTTAACTCTGACTCTTTACTCTATTATAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGACGTTTTGCAGATATCGGTGTGAATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTTGGGAGAAATGTACGGGTTCTTCGACTGGGaa                                                                                                                                                                          <  1:349310/223‑1 (MQ=255)
                                          aCATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTTGACTGGAACAGCCGGACTACCGAAAAATCGTCGACCAAAACCACCAAAGATGATTTTGGTTATCTTGAAATTGAAAGGGGCGCTCACCTTAGTTGGGGGGAAATGTAAGGGTTTTTCGGATGGGaaaaa                                                                                                                                                                       >  1:173375/1‑228 (MQ=255)
                                                    aCTCTGTCCTTTTCCTTTATATTAGCGATACCTGTTGCATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATAGCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACg                                                                                                                                                           <  2:277592/229‑1 (MQ=255)
                                                         gACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATACTGGTGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGGGGCGCTTACTTTAGTTTGGGAGAAATTTACGGGTTCTTCGACTGGGCAAACTTCTATAACGGTCGc                                                                                                                                                      >  1:100458/1‑229 (MQ=255)
                                                             ttttCCTCTATCTTAGCGATTCCGGTTTTTCTTCATGCAGAATTTAAAGGCGGTGTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTAGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCGTCGACTGGGAAAACTTCTATAACGGTCGCCata                                                                                                                                                  <  2:104076/228‑1 (MQ=255)
                                                                          tAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCCGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACCCAAAGATGATTTTGGTTATCTTTAATTTTAAGGTGGCGGTAACTTTAGTTTGGGAGAAATGTAAGGGTTCTTCGCCTTGGAAAAATTCTATAACGGTCGCCATAATAAAACAGGGAg                                                                                                                                     >  1:230245/1‑229 (MQ=255)
                                                                          tAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTTGACTTGAACAGCCGGACTTCCGGAAAATCGTCGCCCCAATCACACCAAAGTGATTTTTGTTATTTTCAAATTGAAGGGGGCGCTACCTTTCGTTTGGGGGAAATGTACCGGTTTTTCACCCGGGCAAACCTCTTTCAAGGGGGGCCATAAAAACCACGTaa                                                                                                                                     >  2:471374/1‑228 (MQ=255)
                                                                          tAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGGCTACCGGAAAATCGTCGACCAAATCACACAAAGATGATTTTTGGTATCTTGGATTTGAAAGGGGCGCTAACTTTAGTTGGGGGGGAATGTACGGGGTTTTTGACTGGGGAAACTTTTATAAAGGGCGCCCtaataactgtcctt                                                                                                                                      >  1:531617/1‑220 (MQ=255)
                                                                            gCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCAATATCTGGACTGGAACAGCCGGGCTACCCAAAAAACGTCGGACCAATTACACAAAAATAGATTTTGTTATTTTGAAGTTTAAGGGGGGCCCAACTTTAAATGGGGGGAAATTTAAGGGGTTTTTGGATGGGGAAAATTTTTTAACAGTCGCCCAAAAAAAcaaagtgagc                                                                                                                                   >  2:705703/1‑220 (MQ=255)
                                                                              tattccTGTTTGATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGag                                                                                                                                 <  2:442739/225‑1 (MQ=255)
                                                                              gATACCTGTTTGATCTCATTCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGag                                                                                                                                 <  2:486062/229‑1 (MQ=255)
                                                                              cattccTGTTTTATTTCATGCAGAATTTAAAGGCGGTTTTGCAGAAATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGag                                                                                                                                 <  2:230245/225‑1 (MQ=255)
                                                                                 aCCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGGTATCGGCGTGCGTTATCTGGGATGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGGTATCTTGAATTTGAAGGGGGGGCCAACTTTAGTTGGGAGGTAATGTACGGGGACTTAGCATGGGTAAAACTTTATAACAGGCGGCCAAAAAAAACAAGGAGGGagcaa                                                                                                                             >  2:167683/1‑228 (MQ=255)
                                                                                 aCCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTTGGTATCTTTAATTTTAAAGGGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGGTCTTTGACTGGGGAAACTTCTATAACGGTGGCCCTAATAAAACAGGTacggggcag                                                                                                                              >  1:371191/1‑221 (MQ=255)
                                                                                        gtATCTCATTCAGAATTTAAAGGCCGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTGTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTta                                                                                                                        <  1:104047/227‑1 (MQ=255)
                                                                                           tctcATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTTTGCTGGACCAGCCGGACTACCGAACAATTGTGGGCCAAATCACACAAAGGTTATTTTGGTTTTTTTGAATTTGAAGGTGGCGCTAACTTTATTTGGGGAGGAATGTATGGGGTCTTTGACTTGGAAAACTTCTATAACGGTGGCCCTAAAAAACTAGGTTGCGGGCAGCGGTTTACCt                                                                                                                   >  1:640092/1‑229 (MQ=255)
                                                                                                     aaTTTAAAGGCGGTTTTGCAGATATCGGCGGGCATTATCTGGGCTGGACCAGCCGGGCTACCGGAAAAACGTCGGCCAAATCACACAAAGATGATTTTTGTTATCTTTAAATTGAAGGGGGCGCTAACATTAGTTGGGGGGGAAATTACGGGGTTTTTGAATGGGGAAAATTTTATAAAGGGCGCCATAAATAACAAGGGAGGCTGGAAGGGTAAACCTTTaaaaaaaa                                                                                                         >  2:266825/1‑226 (MQ=255)
                                                                                                     aaTTTAAAGGCGGTTTTGCAGATATCGGCGGGACTTTTCTTGACTGGAACAGCCCGACTTCCGGAAAATCGTCGACCAAATTACACAAAAATGAATTTGGTTTATTTGAATTTGAAGGGGGGGCTCACTTTCGATGGGGTGACAAGTACGGGTTCTTCGACTAGGAGAACCTCTAAAACGGTGGCCCTAAAAAAACACGTAGGCAGCGAGGGTTTAACTTTaaaaaaaa                                                                                                         >  2:173162/1‑226 (MQ=255)
                                                                                                        ttAAAGGCGGTTTTGCAAATATCTGCGTTCATTATCTGGACTGGACCAGCCGGACTTCCGAAAAATCGTCGACCAAATCACACAAAGATGATTGTGGTTAGCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCGTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCaa                                                                                                      <  2:236367/229‑1 (MQ=255)
                                                                                                        ttAAAGGCGGTTTTGAAAAAATCGGCGGGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCaa                                                                                                      <  2:708507/229‑1 (MQ=255)
                                                                                                        ttAAAGGAGGTTGTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCaa                                                                                                      <  1:471374/229‑1 (MQ=255)
                                                                                                                  ttttGCAGGTATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGATAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTAc                                                                                            <  1:476045/229‑1 (MQ=255)
                                                                                                                  ttttGCAGATAACGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTAc                                                                                            <  1:532701/229‑1 (MQ=255)
                                                                                                                              cGGCGTGCATTATCTGGACTGAAACGGGCGGGCTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATTTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGTGTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGacac                                                                                 <  1:719498/228‑1 (MQ=255)
                                                                                                                                    gCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGa                                                                             >  2:588278/1‑226 (MQ=255)
                                                                                                                                     cATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGAtt                                                                           <  1:454824/227‑1 (MQ=255)
                                                                                                                                                     ccccggcgagcTACCGAAATATCGTCGACGAAATCACACAAGGGTGATTTTGGTTATGTTGAATTTGAAGGGGGAGCTAACTTTAGTGGGGGAGAAATGTACGGGTTCGTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATATCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATg                                                          <  2:356770/219‑1 (MQ=255)
                                                                                                                                                     acCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTTCAt                                                           >  1:218048/1‑228 (MQ=255)
                                                                                                                                                     acCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGACATAATAAACCAGGAAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTCCCTGGGGGAAACCGGATTTAATCTCTACTTACAt                                                           >  1:729055/1‑228 (MQ=255)
                                                                                                                                                         gCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGAATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTa                                                      >  1:249830/1‑229 (MQ=255)
                                                                                                                                                            ggACTTCCGAAAAATCGGCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACCCCGGATTTAATCGCTACTTGCATGCGTATg                                                    <  1:585534/228‑1 (MQ=255)
                                                                                                                                                             gACTACCGAAAAATCGTCGACCAAATCAAACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGTTTTAATCTCTACTTGCATGCGTATGGc                                                  <  1:705704/229‑1 (MQ=255)
                                                                                                                                                                  ccGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTa                                             >  1:373424/1‑229 (MQ=255)
                                                                                                                                                                                     aTCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAACCCAGGTAGCGAGCAGCGTTATACCTTTAAAACTACCAACCGTATTTACCTGGCAGACACCGGCTTTAATCTCTACTGTCTGCCGAATGGCCCCTACGGTTCTGCGAATagctag                          >  2:34927/1‑226 (MQ=255)
                                                                                                                                                                                     aTCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTg                          <  1:588279/229‑1 (MQ=255)
                                                                                                                                                                                                aGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCac                 >  1:583725/1‑229 (MQ=255)
                                                                                                                                                                                                  atgatTTTGGTTAGCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATCCCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCacga               <  1:143884/229‑1 (MQ=255)
                                                                                                                                                                                                   tttgttttGGTTTTCTTGCATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATTCCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCacga               <  2:173375/225‑1 (MQ=255)
                                                                                                                                                                                                     tgttttGGTTTCTTTGAATTTGAAGGTGGCGTTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATCCCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCacgacg             <  2:583724/227‑1 (MQ=255)
                                                                                                                                                                                                        tttGGTTTTTTTGAATTTGAAGGGGGCCATAACTTTAGTTGGGGAGAAATGTCCGGGTTCTTCGACTGGGAAAACTTCTATAACGTTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATTCCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCACGACGata          <  2:729054/229‑1 (MQ=255)
                                                                                                                                                                                                        tttGGTTATCTTGATTTTGAAGGTGGCGCTAAGTTTAATTGGGGATAAATGTACGGGGTTTTCGCCTGGGAAAGCTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATTCCTTTAAAAATACCAACCGTATTTCCCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCACGACGata          <  2:292069/229‑1 (MQ=255)
                                                                                                                                                                                                        tttGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCACGACGata          <  1:173162/229‑1 (MQ=255)
                                                                                                                                                                                                            gTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAACCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGCATCCCGTGAATTTCCACGCCGATATGtt      >  1:496255/1‑229 (MQ=255)
                                                                                                                                                                                                             ttATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAATCTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGacac                                                                                 <  2:498862/153‑1 (MQ=255)
                                                                                                                                                                                                             ttATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGacac                                                                                 >  1:498862/1‑153 (MQ=255)
                                                                                                                                                                                                                 cTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGGGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGACTTTCCACGACGCTATGTTCCTg  >  1:272990/1‑228 (MQ=255)
                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
TCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTGCATTATCTGGACTGGACCAGCCGGACTACCGAAAAATCGTCGACCAAATCACACAAAGATGATTTTGGTTATCTTGAATTTGAAGGTGGCGCTAACTTTAGTTGGGGAGAAATGTACGGGTTCTTCGACTGGGAAAACTTCTATAACGGTCGCCATAATAAACCAGGTAGCGAGCAGCGTTATACCTTTAAAAATACCAACCGTATTTACCTGGGTGACACCGGATTTAATCTCTACTTGCATGCGTATGGCACCTACGGTTCTGCGAATCGCGTGAATTTCCACGACGATATGTTCCTG  >  NC_000913/2525093‑2525525

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