Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,791 C→G 5.5% E139Q (GAG→CAG)  ecpE ← putative fimbrial chaperone EcpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,7910CG 5.5% 349.8 / 7.0 104E139Q (GAG→CAG) ecpEputative fimbrial chaperone EcpE
Reads supporting (aligned to +/- strand):  ref base C (53/44);  new base G (3/3);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 105 total aligned reads displayed.

CGGTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCA  >  NC_000913/304644‑304940
                                                                                                                                                   |                                                                                                                                                     
cGGTCGAATCACCTCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCg                                                                                                                                                     >  2:548311/1‑150 (MQ=255)
 ggTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGt                                                                                                                                                    >  1:855180/1‑150 (MQ=255)
  gTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTg                                                                                                                                                   >  1:845715/1‑150 (MQ=255)
     gAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTg                                                                                                                                                <  2:172931/150‑1 (MQ=255)
        tCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGAcc                                                                                                                                             <  1:623706/150‑1 (MQ=255)
               cTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCAcgccgc                                                                                                                                      <  2:707863/150‑1 (MQ=255)
                    ttAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGtt                                                                                                                                 <  1:895429/150‑1 (MQ=255)
                     tAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtgg                                                                                                                                <  1:364412/150‑3 (MQ=255)
                       aTCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGGCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTggg                                                                                                                              <  1:932805/150‑1 (MQ=255)
                         cAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtggtggt                                                                                                                            <  2:630839/150‑7 (MQ=255)
                            tAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCg                                                                                                                         <  2:696928/150‑1 (MQ=255)
                              gCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTc                                                                                                                       >  1:1374084/1‑150 (MQ=255)
                              gCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTc                                                                                                                       >  1:766750/1‑150 (MQ=255)
                                   agccATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTa                                                                                                                  <  1:863400/148‑1 (MQ=255)
                                   aaCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTa                                                                                                                  <  2:810504/150‑1 (MQ=255)
                                   aaCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTa                                                                                                                  >  1:1563314/1‑150 (MQ=255)
                                     ccATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATg                                                                                                                >  1:669174/1‑150 (MQ=255)
                                       aTGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGt                                                                                                              <  1:799848/150‑1 (MQ=255)
                                       aTGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGACGCCGGTTGGGCTGCGTCTTGTATGGt                                                                                                              >  1:1524230/1‑150 (MQ=255)
                                           gtTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACg                                                                                                          <  2:739931/150‑1 (MQ=255)
                                             tGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAg                                                                                                        <  1:1210928/150‑1 (MQ=255)
                                               ccAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTg                                                                                                      <  2:140209/150‑1 (MQ=255)
                                                 aGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTggg                                                                                                    <  2:770497/150‑1 (MQ=255)
                                                  gtgtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTgggg                                                                                                   <  1:45936/150‑1 (MQ=255)
                                                      tacttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCt                                                                                               >  2:517637/1‑150 (MQ=255)
                                                      tacttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCt                                                                                               >  2:1026479/1‑150 (MQ=255)
                                                      tacttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCt                                                                                               >  2:484163/1‑150 (MQ=255)
                                                       acttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTc                                                                                              >  1:1545604/1‑150 (MQ=255)
                                                       acttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTc                                                                                              >  1:50339/1‑150 (MQ=255)
                                                       acttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACAACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTc                                                                                              >  2:1128297/1‑150 (MQ=255)
                                                        cttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcg                                                                                             >  1:367936/1‑150 (MQ=255)
                                                         ttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgc                                                                                            <  1:44722/150‑1 (MQ=255)
                                                          tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                           >  1:795231/1‑150 (MQ=255)
                                                          tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                           <  2:845715/150‑1 (MQ=255)
                                                          tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                           <  2:709323/150‑1 (MQ=255)
                                                          tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                           >  2:312221/1‑150 (MQ=255)
                                                           acCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGa                                                                                          >  1:1223035/1‑150 (MQ=255)
                                                             cGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGAcc                                                                                                >  2:1621495/1‑142 (MQ=255)
                                                             cGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGAcc                                                                                                <  1:1621495/142‑1 (MQ=255)
                                                             cGTTCCTGTCAACTTGTCGAAGGACAATTTAAACTGAAACTGAAGTGGTCGTACTACCAGAATGGTATCAATCACCACAACCGCCTCAGTGCAGACCACGACGCCGGGTGGGCTGAGTCTTGTATGGTTACGAGTGGGGACCTAGCGaaa                                                                                        >  2:1486793/1‑150 (MQ=255)
                                                              gTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAAt                                                                                       <  1:1477458/150‑1 (MQ=255)
                                                                 ccTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgag                                                                                    >  2:5216/1‑150 (MQ=255)
                                                                 ccTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgag                                                                                    >  2:1301537/1‑150 (MQ=255)
                                                                  cTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgaga                                                                                   >  1:666419/1‑150 (MQ=255)
                                                                  cTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgaga                                                                                   >  1:312398/1‑150 (MQ=255)
                                                                   tGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGAc                                                                                  <  2:762026/150‑1 (MQ=255)
                                                                      cACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCg                                                                               >  2:196871/1‑150 (MQ=255)
                                                                        ccTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTTCGTCTTGTATGGTTACGACTGGGGACAAACCGAAATAAGAcacagt                                                                             >  2:693054/1‑145 (MQ=255)
                                                                        ccTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgt                                                                             >  2:429169/1‑150 (MQ=255)
                                                                         cTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgta                                                                            <  2:1223035/150‑1 (MQ=255)
                                                                           tGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgtagt                                                                          >  1:1399069/1‑150 (MQ=255)
                                                                           tGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgtagt                                                                          >  1:1314540/1‑150 (MQ=255)
                                                                             tCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGACGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGtagtag                                                                        <  2:1524230/150‑1 (MQ=255)
                                                                              cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                       >  2:79282/1‑150 (MQ=255)
                                                                              cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                       >  1:1137575/1‑150 (MQ=255)
                                                                              cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                       >  1:870575/1‑150 (MQ=255)
                                                                              cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                       >  1:700901/1‑150 (MQ=255)
                                                                               gAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgcg                                                                      <  1:300567/150‑1 (MQ=255)
                                                                                  ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCTCCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                   >  1:134823/1‑150 (MQ=255)
                                                                                    aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                                 <  1:1518348/150‑1 (MQ=255)
                                                                                     ccATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgcg                                                                <  1:1337303/150‑1 (MQ=255)
                                                                                       aTTTAACCTGAACCTGCCGTGGTCGTTCTACCAGATTGGTATCGATCTCGTCCTGCGGCTGCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGt                                                              <  1:693054/150‑1 (MQ=255)
                                                                                       aTTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGAGCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGt                                                              >  2:1172147/1‑150 (MQ=255)
                                                                                           aaaCTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATc                                                          >  2:1623338/1‑150 (MQ=255)
                                                                                             aCTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATccc                                                        >  2:1288153/1‑150 (MQ=255)
                                                                                               tGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGt                                                      <  1:1301537/150‑1 (MQ=255)
                                                                                                  aCCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGATGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGa                                                   >  1:1139336/1‑150 (MQ=255)
                                                                                                   ccTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAc                                                  <  2:1314540/150‑1 (MQ=255)
                                                                                                    cTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTTGCGGTTATCCCGTGGAcc                                                 <  1:15533/150‑1 (MQ=255)
                                                                                                      gACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTTGACCAt                                               <  1:142108/150‑1 (MQ=255)
                                                                                                        cGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGa                                             >  2:76034/1‑150 (MQ=255)
                                                                                                         gTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGAt                                            <  1:79282/150‑1 (MQ=255)
                                                                                                         gTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGAt                                            <  1:879174/150‑1 (MQ=255)
                                                                                                           ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                          <  2:842720/150‑1 (MQ=255)
                                                                                                           ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                          <  1:429521/150‑1 (MQ=255)
                                                                                                                   tacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAttt                                  >  2:1480406/1‑150 (MQ=255)
                                                                                                                   tacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAttt                                  >  1:480482/1‑150 (MQ=255)
                                                                                                                   tacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTAATCCCGTGGACCATGATAGTAAAAttt                                  >  1:939805/1‑150 (MQ=255)
                                                                                                                    acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTa                                 >  1:175297/1‑150 (MQ=255)
                                                                                                                     ccAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaa                                <  1:1115221/150‑1 (MQ=255)
                                                                                                                     ccAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaa                                <  2:1399069/150‑1 (MQ=255)
                                                                                                                      cAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaa                               <  1:429169/150‑1 (MQ=255)
                                                                                                                        gAATGGTATCCATCACCTCCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAAt                             <  2:134823/150‑1 (MQ=255)
                                                                                                                          aTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                           <  1:905440/150‑1 (MQ=255)
                                                                                                                             gTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTTGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCg                        <  2:855180/150‑1 (MQ=255)
                                                                                                                                    tCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCAc                 >  2:757798/1‑150 (MQ=255)
                                                                                                                                      accaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCa               >  1:247082/1‑150 (MQ=255)
                                                                                                                                      accaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCa               >  2:132262/1‑150 (MQ=255)
                                                                                                                                      accaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCa               >  2:1104924/1‑150 (MQ=255)
                                                                                                                                       gcccaaccGGCGGCGTGGTCAGCACGGAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAg              >  2:1127445/6‑150 (MQ=255)
                                                                                                                                       ccaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAg              >  2:355510/1‑150 (MQ=255)
                                                                                                                                       ccaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAg              >  2:727909/1‑150 (MQ=255)
                                                                                                                                            accGGCGGCGTGGTCAGCACGGAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGc         >  1:722610/1‑150 (MQ=255)
                                                                                                                                             ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCa        >  1:958371/1‑150 (MQ=255)
                                                                                                                                             ccGGCGGCGTGGTCAGCACGGAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCa        >  2:1603832/1‑150 (MQ=255)
                                                                                                                                              cGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCaa       <  1:355510/150‑1 (MQ=255)
                                                                                                                                               ggCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAAc      <  1:1480406/150‑1 (MQ=255)
                                                                                                                                                gCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACg     <  1:508707/150‑1 (MQ=255)
                                                                                                                                                gCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACg     <  2:939805/150‑1 (MQ=255)
                                                                                                                                                   ccGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCCCACCAGCCTGCAACGCCa  <  1:280238/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCA  >  NC_000913/304644‑304940

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