Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,800 C→G 5.6% V136L (GTC→CTC)  ecpE ← putative fimbrial chaperone EcpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,8000CG 5.6% 352.4 / 6.4 103V136L (GTC→CTC) ecpEputative fimbrial chaperone EcpE
Reads supporting (aligned to +/- strand):  ref base C (51/46);  new base G (4/2);  total (55/48)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

TCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCG  >  NC_000913/304652‑304947
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                         >  1:1563314/1‑150 (MQ=255)
                           aaCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTa                                                                                                                         <  2:810504/150‑1 (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:1190526/1‑150 (MQ=255)
                                              tacttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCt                                                                                                      >  2:1026479/1‑150 (MQ=255)
                                              tacttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCt                                                                                                      >  2:484163/1‑150 (MQ=255)
                                              tacttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCt                                                                                                      >  2:517637/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                                                                                                  <  2:845715/150‑1 (MQ=255)
                                                  tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                                  <  2:709323/150‑1 (MQ=255)
                                                  tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                                  >  1:795231/1‑150 (MQ=255)
                                                  tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                                  >  2:312221/1‑150 (MQ=255)
                                                   acCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGa                                                                                                 >  1:1223035/1‑150 (MQ=255)
                                                     cGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGAcc                                                                                                       <  1:1621495/142‑1 (MQ=255)
                                                     cGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGAcc                                                                                                       >  2:1621495/1‑142 (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)
                                                                ccTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTTCGTCTTGTATGGTTACGACTGGGGACAAACCGAAATAAGAcacagt                                                                                    >  2:693054/1‑145 (MQ=255)
                                                                ccTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgt                                                                                    >  2:429169/1‑150 (MQ=255)
                                                                ccTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCAGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgt                                                                                    >  2:513207/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                                                                              >  1:1137575/1‑150 (MQ=255)
                                                                      cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                              >  2:79282/1‑150 (MQ=255)
                                                                      cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                              >  1:870575/1‑150 (MQ=255)
                                                                      cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                              >  1:700901/1‑150 (MQ=255)
                                                                          ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCTCCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                          >  1:134823/1‑150 (MQ=255)
                                                                            aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                                        <  1:1518348/150‑1 (MQ=255)
                                                                             ccATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgcg                                                                       <  1:1337303/150‑1 (MQ=255)
                                                                              cATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCgg                                                                      >  1:163961/1‑150 (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)
                                                                                                   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                                       <  2:1399069/150‑1 (MQ=255)
                                                                                                             ccAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaa                                       <  1:1115221/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                      >  2:132262/1‑150 (MQ=255)
                                                                                                                              accaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCa                      >  1:247082/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)
                                                                                                                                             gTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGc       <  1:1190526/150‑1 (MQ=255)
                                                                                                                                              tGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCt      >  2:894885/1‑150 (MQ=255)
                                                                                                                                               gCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTg     >  1:604283/1‑150 (MQ=255)
                                                                                                                                                  gACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCg  <  1:76034/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCG  >  NC_000913/304652‑304947

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: