Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,800 C→G 5.2% V136L (GTC→CTC)  ecpE ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,8000CG 5.2% 492.7 / 12.7 149V136L (GTC→CTC) ecpEECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (79/62);  new base G (3/5);  total (82/67)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 150 total aligned reads displayed.

ATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGG  >  NC_000913/304651‑304948
                                                                                                                                                     |                                                                                                                                                    
aTCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGAc                                                                                                                                                      <  2:829656/150‑1 (MQ=255)
  cacaTCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCa                                                                                                                                                    <  1:351232/150‑1 (MQ=255)
  cacaTCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCa                                                                                                                                                    >  1:1004991/1‑150 (MQ=255)
   acaTCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACTTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCAc                                                                                                                                                   >  1:1160169/1‑150 (MQ=255)
      tCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCAcgcc                                                                                                                                                <  2:489192/150‑1 (MQ=255)
       ccTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGACCCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGCCGTGCTGACCccgccg                                                                                                                                               >  2:131110/1‑150 (MQ=255)
        cTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACAAGTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACAACCGGCTCCGTGCTGACCAcgccgc                                                                                                                                              >  1:49811/1‑150 (MQ=255)
         tGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCCCCACCGGCTCCGTGCTGACCACgccgcc                                                                                                                                             >  2:929249/1‑150 (MQ=255)
                aTCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCCCCTTGTCGAAGGACCATTTAAACTGAACCTGACGGGGTCGTACTACCAGATTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGGTggg                                                                                                                                      >  1:1617913/1‑150 (MQ=255)
                aTCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGCTCCTGTCACCTTGTCGAAGGACCATTTACACTGAACCTGACGTGGGCGTACTACCAGAATGGTATCCAGCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtggtg                                                                                                                                      <  1:398609/150‑5 (MQ=255)
                    gTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGc                                                                                                                                  >  2:213093/1‑150 (MQ=255)
                      aGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGt                                                                                                                                >  1:467815/1‑150 (MQ=255)
                      aGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGt                                                                                                                                >  1:245982/1‑150 (MQ=255)
                             aCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTAt                                                                                                                         <  2:1089089/150‑1 (MQ=255)
                             aCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTAt                                                                                                                         >  1:24506/1‑150 (MQ=255)
                                    gtTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACg                                                                                                                  >  2:1824850/1‑150 (MQ=255)
                                    gtTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTACGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACg                                                                                                                  >  2:403062/1‑150 (MQ=255)
                                      tGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAg                                                                                                                >  2:1193508/1‑150 (MQ=255)
                                         cAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTgg                                                                                                             <  1:1074519/150‑1 (MQ=255)
                                         cAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTgg                                                                                                             >  2:1531189/1‑150 (MQ=255)
                                         cAGTCTTAGTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCTTCTTGTATTGTTACGAGTgg                                                                                                             <  1:1089963/150‑1 (MQ=255)
                                           gtgtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTgggg                                                                                                           <  2:30329/150‑1 (MQ=255)
                                             gtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGAc                                                                                                         <  2:153695/150‑1 (MQ=255)
                                                   tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTATTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                                   >  1:1176993/1‑150 (MQ=255)
                                                    acCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGa                                                                                                  <  2:1812736/150‑1 (MQ=255)
                                                      cGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGaaa                                                                                                >  1:1804277/1‑150 (MQ=255)
                                                       gTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAAt                                                                                               <  1:1550980/150‑1 (MQ=255)
                                                          ccTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgag                                                                                            <  2:1398955/150‑1 (MQ=255)
                                                            tGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGAc                                                                                          <  1:902381/150‑1 (MQ=255)
                                                               cACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCAGCCGGTTGGGCTGCGGCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCg                                                                                       >  2:1084088/1‑150 (MQ=255)
                                                                aCCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCCCCACCACCGGCTCCGAGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCgg                                                                                      >  1:99219/1‑150 (MQ=255)
                                                                 ccTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGTTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgt                                                                                     <  2:1417159/150‑1 (MQ=255)
                                                                  cTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgta                                                                                    >  1:314009/1‑150 (MQ=255)
                                                                   ttGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCTGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgtag                                                                                   <  2:752044/150‑1 (MQ=255)
                                                                   ttGTCGAAGGACCATTTAAACTCAACCTGCCGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCGGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgtag                                                                                   <  1:496998/150‑1 (MQ=255)
                                                                      tCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGtagtag                                                                                >  1:777552/1‑150 (MQ=255)
                                                                       cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                               >  1:1746345/1‑150 (MQ=255)
                                                                          aGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGCCGTGCTGACCACGCCGCCGCTTGGGCTGCGGCTTGTATGGTTCCGAGTTGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCt                                                                            <  2:907248/150‑1 (MQ=255)
                                                                           ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                           <  1:1571296/150‑1 (MQ=255)
                                                                           ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                           >  1:1678450/1‑150 (MQ=255)
                                                                           ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGGATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                           >  1:414530/1‑150 (MQ=255)
                                                                             aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                                         >  1:1716419/1‑150 (MQ=255)
                                                                               cATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCgg                                                                       >  1:1722318/1‑150 (MQ=255)
                                                                                  ttAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTGGCGCTCGCGGTTa                                                                    <  2:26945/150‑1 (MQ=255)
                                                                                       cTGAACCTGATGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCGGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAAGGAGACCCGGTCGGAGCGCTGGCGGTTTTCCCg                                                               >  2:1099992/1‑150 (MQ=255)
                                                                                       cTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCg                                                               >  1:580642/1‑150 (MQ=255)
                                                                                         gacccTGACGGGGTCGTACTACCAGACTGGAATCCATCACCACCACCCGCTCCGGGCGGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTCACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTg                                                             <  1:1668138/147‑1 (MQ=255)
                                                                                            ccTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAc                                                          >  1:606887/1‑150 (MQ=255)
                                                                                             cTGACGTGGTCGATCTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAcc                                                         <  2:734977/150‑1 (MQ=255)
                                                                                               gACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCAt                                                       <  2:24506/150‑1 (MQ=255)
                                                                                                aCGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATg                                                      <  2:149036/150‑1 (MQ=255)
                                                                                                aCGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATg                                                      <  2:694225/150‑1 (MQ=255)
                                                                                                    gggcGTACTACCAGACGGGGATCCATCACCACCACCGGCGCCGTGCGGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                                  <  1:131110/147‑1 (MQ=255)
                                                                                                    ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGTTCCGTGCTGACCTCGCCGCCGGTTGGGCTGCGTCTTTTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTTGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                                  <  2:1669634/150‑1 (MQ=255)
                                                                                                    ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                                  <  1:526607/150‑1 (MQ=255)
                                                                                                    ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                                  >  1:316502/1‑150 (MQ=255)
                                                                                                      tCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTa                                                >  1:1035707/1‑150 (MQ=255)
                                                                                                        gTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTAGCGGTTATCCCGTGGACCATGATAGTaaa                                              >  2:348088/1‑150 (MQ=255)
                                                                                                         tactacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTaaaa                                             >  2:330988/1‑150 (MQ=255)
                                                                                                          actacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAt                                            <  1:1373190/150‑1 (MQ=255)
                                                                                                           ctacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAtt                                           >  1:1108227/1‑150 (MQ=255)
                                                                                                           ctacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAtt                                           >  1:1655983/1‑150 (MQ=255)
                                                                                                            tacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAttt                                          >  1:849377/1‑150 (MQ=255)
                                                                                                             acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGCCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTa                                         >  1:418450/1‑150 (MQ=255)
                                                                                                             acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTa                                         >  1:433114/1‑150 (MQ=255)
                                                                                                             acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTa                                         >  1:1246290/1‑150 (MQ=255)
                                                                                                              ccAGAATGGTATGCATCTGCAGCACGGGCTCCGTGCTGACCACGCCGGCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGATTTGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATATTAAAATTTaa                                        <  2:1784596/150‑1 (MQ=255)
                                                                                                              ccAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaa                                        >  1:405122/1‑150 (MQ=255)
                                                                                                               cAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaa                                       >  2:1717524/1‑150 (MQ=255)
                                                                                                               cAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaa                                       <  2:1717335/150‑1 (MQ=255)
                                                                                                                aGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaaa                                      >  2:1249406/1‑150 (MQ=255)
                                                                                                                aGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaaa                                      <  1:1297283/150‑1 (MQ=255)
                                                                                                                   aTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                                   >  2:400048/1‑150 (MQ=255)
                                                                                                                   aTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGCCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                                   >  2:374767/1‑150 (MQ=255)
                                                                                                                    tGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACt                                  <  1:705145/150‑1 (MQ=255)
                                                                                                                       tATCCATCACCACCTCCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCTGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGc                               <  1:1354976/150‑1 (MQ=255)
                                                                                                                        aTCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGct                              <  1:621289/150‑1 (MQ=255)
                                                                                                                         tCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGctc                             <  2:1459604/150‑1 (MQ=255)
                                                                                                                          ccctCACCACCACCGGCTCCGTGCTGACCAAGCCGCCGGTTGGGCTGCGTCTTGTATTGTTACGAGTGGGGACCTCGCGACATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGctct                            <  2:303005/147‑1 (MQ=255)
                                                                                                                            atcaccatcaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCa                          <  1:917822/150‑1 (MQ=255)
                                                                                                                               accaccaccGTCTCCGTGCTGACCACGCCGCGTGTTGGGCTGCGTCTTGTATGGTTTCGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCa                       <  2:278401/150‑1 (MQ=255)
                                                                                                                                  accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTATAATTTAAAATACTCGCTCTCACCAGcc                    >  2:1127370/1‑150 (MQ=255)
                                                                                                                                  accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc                    <  1:658667/150‑1 (MQ=255)
                                                                                                                                  accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc                    >  1:1254380/1‑150 (MQ=255)
                                                                                                                                  accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc                    <  1:671815/150‑1 (MQ=255)
                                                                                                                                  accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATAATAGTAAAATTTAAAATACTCGCTCTCACCAGcc                    >  2:710220/1‑150 (MQ=255)
                                                                                                                                   ccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCt                   <  1:929249/150‑1 (MQ=255)
                                                                                                                                   ccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCt                   >  2:1009381/1‑150 (MQ=255)
                                                                                                                                   caaccGGCGGCGTGGTCAGCACGGAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGCGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGAAAAATTTAAAATACTCCCTCTCCCCTGCCt                   >  1:1138819/3‑150 (MQ=255)
                                                                                                                                    caccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCAGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTATCACCAGCATg                  >  2:644118/1‑150 (MQ=255)
                                                                                                                                     aacgtcTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTTTGGTTACGAGTGGGGACCTCGCGCAATGAGTCCCGGTAGTAGGGCGCGCGGTTCTCCCGTGGACCAAGATAGTAAAATTTAAAATACACGCTCTCCCCAGCCTGc                 >  2:75017/6‑150 (MQ=255)
                                                                                                                                      ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCa                <  2:1225339/150‑1 (MQ=255)
                                                                                                                                        ggCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAAc              <  1:523364/150‑1 (MQ=255)
                                                                                                                                          cTCCGTGCTGAGCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGCGCCCCGGTAGTAGCGCTCGCGGTTATCCCGGGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGc            >  2:1647655/1‑150 (MQ=255)
                                                                                                                                          cTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGc            <  2:498260/150‑1 (MQ=255)
                                                                                                                                           tCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGcc           >  1:365265/1‑150 (MQ=255)
                                                                                                                                             ggTGCTGACCACGCCGCCGGTTGGGCTGTGTCTTGTATGGATACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAg         <  2:810493/149‑1 (MQ=255)
                                                                                                                                              gTGCTGACCACGCCGCCGGTTGGACTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGc        >  1:752211/1‑150 (MQ=255)
                                                                                                                                                  tGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGc    >  1:936662/1‑150 (MQ=255)
                                                                                                                                                    aCCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCgg  <  1:54940/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGG  >  NC_000913/304651‑304948

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