Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,522,716 T→A 5.4% D117V (GAC→GTC)  yncD ← putative iron outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,522,7160TA 5.4% 905.0 / 28.5 254D117V (GAC→GTC) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (115/125);  new base A (10/4);  total (125/129)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-02

Only 100 of 256 total aligned reads displayed.

CGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTATCCAGTTCTGAAACCACCTGCGGTGCGGCACTGACAATCATAGTCTGTT  >  NC_000913/1522448‑1522989
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
cGTAAAACGCTTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtggt                                                                                                                                                                                                                                                                         <  1:91567/279‑1 (MQ=255)
                                             tgtgCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGGGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGACCTTCCACATTTTACACACTCCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCTGGCATGGTGGCGGTACAACCGTCCACATACAGGCGACTCGCGCGAATAGCGTa                                                                                                                                                                                                                            >  1:283838/1‑279 (MQ=255)
                                                                                            tGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                                                                                                                                                <  2:609028/212‑1 (MQ=255)
                                                                                                                          aaaTGGTTGGTGGTGTTTTTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTg                                                                                                                                                    <  2:92770/273‑1 (MQ=255)
                                                                                                                           aaTGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCgg                                                                                                                                                                                                                                                                                >  6:874445/1‑149 (MQ=255)
                                                                                                                             tggttggtGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGc                                                                                                                                                                                                                                                                              <  8:932187/149‑1 (MQ=255)
                                                                                                                               gttggtGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAt                                                                                                                                                                                                                                                                            <  7:241749/149‑1 (MQ=255)
                                                                                                                                 tggtggCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtgg                                                                                                                                                                                                                                                                          <  5:138353/149‑1 (MQ=255)
                                                                                                                                 tggtggCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTTTACCTGCAAACCAg                                                                                                                                            >  2:623998/1‑275 (MQ=255)
                                                                                                                                      gCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCg                                                                                                                                                                                                                                                                     <  7:352624/149‑1 (MQ=255)
                                                                                                                                         gttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGa                                                                                                                                                                                                                                                                  <  7:1263410/149‑1 (MQ=255)
                                                                                                                                          ttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGaa                                                                                                                                                                                                                                                                 <  7:1689495/149‑1 (MQ=255)
                                                                                                                                                 gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                                                                                                                                                          <  7:852088/149‑1 (MQ=255)
                                                                                                                                                 gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                                                                                                                                                          <  8:342561/149‑1 (MQ=255)
                                                                                                                                                  tCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCa                                                                                                                                                                                                                                                                             >  4:217362/1‑129 (MQ=255)
                                                                                                                                                   cTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCca                                                                                                                                                                                                                                                        <  5:1044255/149‑1 (MQ=255)
                                                                                                                                                   cTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACGGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCca                                                                                                                                                                                                                                                        <  7:300271/149‑1 (MQ=255)
                                                                                                                                                     gggTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCcaca                                                                                                                                                                                                                                                      >  5:682969/1‑149 (MQ=255)
                                                                                                                                                      ggtggtGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtggt                                                                                                                                                                                                                                                                         <  2:217365/129‑1 (MQ=255)
                                                                                                                                                         ggtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACa                                                                                                                                                                                                                                                  <  6:1604959/149‑1 (MQ=255)
                                                                                                                                                         ggtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACa                                                                                                                                                                                                                                                  <  8:126219/149‑1 (MQ=255)
                                                                                                                                                           tgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCggg                                                                                                                                                                                                                                                                   >  7:1411328/1‑130 (MQ=255)
                                                                                                                                                           tgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                                                                                                                                                <  8:707605/149‑1 (MQ=255)
                                                                                                                                                           tgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                                                                                                                                                <  8:1631595/149‑1 (MQ=255)
                                                                                                                                                           tgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAg                                                                                                                  >  2:306944/1‑275 (MQ=255)
                                                                                                                                                               attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATGCAGgcgaa                                                                                                                                                                                                                                            <  7:648051/149‑1 (MQ=255)
                                                                                                                                                               attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaa                                                                                                                                                                                                                                            >  6:122094/1‑149 (MQ=255)
                                                                                                                                                                ttcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCa                                                                                                                                                                                                                                                                             >  2:196364/1‑115 (MQ=255)
                                                                                                                                                                    ttACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtggt                                                                                                                                                                                                                                                                         <  4:196361/115‑1 (MQ=255)
                                                                                                                                                                     tACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgc                                                                                                                                                                                                                                      <  7:439218/149‑1 (MQ=255)
                                                                                                                                                                     tACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgc                                                                                                                                                                                                                                      <  7:66618/149‑1 (MQ=255)
                                                                                                                                                                        cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGTCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                                                                                                                                                   >  8:1191334/1‑149 (MQ=255)
                                                                                                                                                                        cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                                                                                                                                                          >  5:1433340/1‑126 (MQ=255)
                                                                                                                                                                        cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                                                                                                                                                   >  8:1224048/1‑149 (MQ=255)
                                                                                                                                                                        cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTGGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                                                                                                                                                   >  7:1454917/1‑149 (MQ=255)
                                                                                                                                                                           accaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATAACGGCGAATACCGCGAATACCGTAAGTGGAACGGTAGCCAAATCCGCGAATCGACAGCCGTAAATCCTGCCAATAGTTCTGCCGGTTTTTTACCTACAAACCAGGCACCCCGGTCAGTTATTCGCACAAATTATTGCTCGGTGtt                                                                                                   >  4:112764/1‑274 (MQ=255)
                                                                                                                                                                              aGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGt                                                                                                                                                                                                                             >  8:506172/1‑149 (MQ=255)
                                                                                                                                                                                         cATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCgg                                                                                                                                                                                                                  <  8:820990/149‑1 (MQ=255)
                                                                                                                                                                                            gtaggggctcgaAGGGGACACGGGGCCCTTCCCCATTTGTCACGCTGATTAATACGATGTGGGATGCTTGCCCTTGCCCGTCGGGCAGCGTGGCGGGAAATCCGGCCCCATCCAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAAGCCTGGACTGGACAGCTGTAAATCTTGCGCATAAGTCAGCCGGTTTTGTACCTGCAAACCAGGGACGCCGGTCAGTGCTGCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCTCCTCg                                                                                 <  4:430756/265‑1 (MQ=255)
                                                                                                                                                                                              aGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGcc                                                                                                                                                                                                             <  8:660043/149‑1 (MQ=255)
                                                                                                                                                                                               gggCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCa                                                                                                                                                                                                            >  5:443889/1‑149 (MQ=255)
                                                                                                                                                                                               gggCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGCGGAGCGGGAGCCa                                                                                                                                                                                                            >  8:1360311/1‑149 (MQ=255)
                                                                                                                                                                                                   agagAAGGGGCCAAGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATc                                                                                                                                                                                                        >  8:55307/1‑149 (MQ=255)
                                                                                                                                                                                                     agaANGGGNCACGCAGNACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCg                                                                                                                                                                                                                   <  3:40224/136‑1 (MQ=255)
                                                                                                                                                                                                      gaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTcc                                                                                                                                                                                                                                                         <  2:508875/97‑1 (MQ=255)
                                                                                                                                                                                                       aaGGGGTCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcgcg                                                                                                                                                                                                    <  7:1191334/149‑1 (MQ=255)
                                                                                                                                                                                                           ggCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATc                                                                                                                                                                                                <  8:1600681/149‑1 (MQ=255)
                                                                                                                                                                                                            gCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCggg                                                                                                                                                                                                                 <  5:1622472/131‑1 (MQ=255)
                                                                                                                                                                                                            gCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCg                                                                                                                                                                                               >  6:197163/1‑149 (MQ=255)
                                                                                                                                                                                                              cccGCAGCACTTCAACATTTTGAACACTGCTGAAACGGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCTCCTCGCCATCCACCACGCTTAc                                                                <  2:368030/272‑1 (MQ=255)
                                                                                                                                                                                                              cACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGAc                                                                                                                                                                                             >  8:264593/1‑149 (MQ=255)
                                                                                                                                                                                                                cGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACat                                                                                                                                                                                           <  8:1454917/149‑2 (MQ=255)
                                                                                                                                                                                                                cGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCg                                                                                                                                    >  2:48735/1‑204 (MQ=255)
                                                                                                                                                                                                                  cagcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCt                                                                                                                                                                                         <  6:443889/149‑1 (MQ=255)
                                                                                                                                                                                                                   agcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGTCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGt                                                                                                                                                                                                                         <  8:855940/116‑1 (MQ=255)
                                                                                                                                                                                                                   agcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTg                                                                                                                                                                                        >  5:209975/1‑149 (MQ=255)
                                                                                                                                                                                                                    gcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGt                                                                                                                                                                                       >  6:1126917/1‑149 (MQ=255)
                                                                                                                                                                                                                    gcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTGGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCa                                                                                                                                <  4:48735/204‑1 (MQ=255)
                                                                                                                                                                                                                      aCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTaa                                                                                                                                                                                     >  6:63208/1‑149 (MQ=255)
                                                                                                                                                                                                                       cTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTaaa                                                                                                                                                                                    <  8:1559128/149‑1 (MQ=255)
                                                                                                                                                                                                                        ttCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAAt                                                                                                                                                                                   >  6:238591/1‑149 (MQ=255)
                                                                                                                                                                                                                          ccACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCgg                                                                                                                                                                                                                                                                    >  6:801684/1‑66 (MQ=255)
                                                                                                                                                                                                                          ccACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCt                                                                                                                                                                                 >  7:794433/1‑149 (MQ=255)
                                                                                                                                                                                                                            acaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGGGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                                                                                                                                               <  5:783202/149‑1 (MQ=255)
                                                                                                                                                                                                                            acaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                                                                                                                                               >  6:149443/1‑149 (MQ=255)
                                                                                                                                                                                                                             caTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCTGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAACCAAATCCGCCAATCGACAGCGTTAAATCTGGAGCATAGTTCTGCAGGTTTGGTACCTGCAAACAAGGCAACGCGGTCAGTGAGTCGGCCAAGTTACTCACCGGTGTTGCGCGGCGCATCTCCTCGCCGTCCACCCCGCTTCATGCTCCCTGGGTagc                                                 >  4:66263/1‑272 (MQ=255)
                                                                                                                                                                                                                              aTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaatac                                                                                                                                                                                                                                >  5:815891/1‑98 (MQ=255)
                                                                                                                                                                                                                              aTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATc                                                                                                                                                                                                <  7:505530/130‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGtt                                                                                                                                                                      >  2:697012/1‑152 (MQ=255)
                                                                                                                                                                                                                                  tGCACACGGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGcc                                                                                                <  1:184836/222‑1 (MQ=255)
                                                                                                                                                                                                                                       acTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTcc                                                                                                                                                                                                                                                         <  6:1024557/64‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTg                                                                                                                                                                   >  5:644210/1‑149 (MQ=255)
                                                                                                                                                                                                                                        cTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCTCCTCGCCATCCACCACGCTTACTGCTGCTggg                                                      >  1:69385/1‑258 (MQ=255)
                                                                                                                                                                                                                                          gCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGcc                                                                                                                                                                 <  6:923873/149‑1 (MQ=255)
                                                                                                                                                                                                                                                aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                                                                                                                                                           <  6:1357242/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCGATGTTGGATGTTTGCCCTTGCCCGTCGGGTATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTg                                                                                                                                                          >  7:1050157/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 tCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTg                                                                                                                                                          <  5:874445/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                                                                                                                                                >  6:727512/1‑93 (MQ=255)
                                                                                                                                                                                                                                                   gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                                                                                                                                                <  5:727512/93‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTAc                                                                                                                                                       <  6:697077/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaatac                                                                                                                                                                                                                                         <  1:145553/66‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTAcc                                                                                                                                                      <  7:1601430/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGAc                                                                                                                     >  1:218398/1‑182 (MQ=255)
                                                                                                                                                                                                                                                      gTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCt                                                                                                                                                     <  7:928241/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCg                                                                                                                                                                                                                   >  5:757901/1‑86 (MQ=255)
                                                                                                                                                                                                                                                       ttGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTg                                                                                                                                                    >  8:1282658/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          gATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCaa                                                                                                                                                 >  5:1574983/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          gATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCaa                                                                                                                                                 >  5:1064885/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           aTGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCaaa                                                                                                                                                <  6:947402/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               ttgcccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACNTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCCCCTCCCCACCCACCCCGCTTACTGCTGCTCGGGTATcccggcnc                                          >  3:193322/1‑241 (MQ=255)
                                                                                                                                                                                                                                                                tgcccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAgg                                                                                                                                           <  7:1282658/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 caccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGttt                                                                                                                                                            <  7:1583075/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                 caccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGc                                                                                                                                          >  5:1548163/6‑149 (MQ=255)
                                                                                                                                                                                                                                                                  accatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCa                                                                                                                                         >  5:1184131/5‑149 (MQ=255)
                                                                                                                                                                                                                                                                    cttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTgc                                                                                                                                                                              >  8:474517/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                     ttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGc                                                                                                                                      >  5:1265346/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGAc                                                                                                                                                                                             <  6:473134/89‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTc                                                                                                                                 <  7:996438/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCa                                                                                                                                >  6:277125/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           gTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCAGCTCCTCGCCATccac                                                                         <  2:420402/204‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
CGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTATCCAGTTCTGAAACCACCTGCGGTGCGGCACTGACAATCATAGTCTGTT  >  NC_000913/1522448‑1522989

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