Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,806 C→A 6.1% G134C (GGC→TGC)  ecpE ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,8060CA 6.1% 410.2 / 12.1 141G134C (GGC→TGC) ecpEECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (71/60);  new base A (4/5);  total (75/66)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 153 total aligned reads displayed.

TCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGA  >  NC_000913/304657‑304955
                                                                                                                                                     |                                                                                                                                                     
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)
   tGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACgccgcc                                                                                                                                                    <  2:382219/150‑1 (MQ=255)
     gCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCgg                                                                                                                                                  <  1:1093145/150‑1 (MQ=255)
        tcaTCAGTAGCTTAAACCATGTGCTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACTTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtgg                                                                                                                                               <  2:1090801/148‑3 (MQ=255)
        taataaGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtgg                                                                                                                                               <  2:442162/150‑3 (MQ=255)
        tAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtgg                                                                                                                                               <  1:1265056/150‑3 (MQ=255)
          aTCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTggg                                                                                                                                             >  1:1320427/1‑150 (MQ=255)
          aTCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGCTCCTGTCACCTTGTCGAAGGACCATTTACACTGAACCTGACGTGGGCGTACTACCAGAATGGTATCCAGCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtggtg                                                                                                                                             <  1:398609/150‑5 (MQ=255)
             aGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACAACCGGCGGCGTGGTCAGCACGGAGCCGGtggtggtg                                                                                                                                          <  2:152867/150‑8 (MQ=38)
              gTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGc                                                                                                                                         >  2:213093/1‑150 (MQ=255)
                aGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGt                                                                                                                                       >  1:467815/1‑150 (MQ=255)
                 gCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGta                                                                                                                                      >  1:149036/1‑149 (MQ=255)
                 gCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTc                                                                                                                                      >  1:694225/1‑150 (MQ=255)
                    tAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCGGCGTCTTg                                                                                                                                   >  1:319152/1‑150 (MQ=255)
                     aaaCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGt                                                                                                                                  <  2:119286/150‑1 (MQ=255)
                       aCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTAt                                                                                                                                <  2:1089089/150‑1 (MQ=255)
                       aCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTAt                                                                                                                                >  1:24506/1‑150 (MQ=255)
                                   cAGTCTTAGTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCTTCTTGTATTGTTACGAGTgg                                                                                                                    <  1:1089963/150‑1 (MQ=255)
                                       gtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGAc                                                                                                                <  2:153695/150‑1 (MQ=255)
                                       gtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGGCTTGTATGGTTACGAGTGGGGAc                                                                                                                <  1:1822296/150‑1 (MQ=255)
                                          acttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTc                                                                                                             <  2:1044181/150‑1 (MQ=255)
                                             tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTATTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                                          >  1:1176993/1‑150 (MQ=255)
                                              acCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGa                                                                                                         <  2:1812736/150‑1 (MQ=255)
                                              acCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGa                                                                                                         <  2:1822916/150‑1 (MQ=255)
                                                    ccTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgag                                                                                                   <  2:1391339/150‑1 (MQ=255)
                                                     cTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgaga                                                                                                  >  1:1571767/1‑150 (MQ=255)
                                                      tGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGAc                                                                                                 <  1:902381/150‑1 (MQ=255)
                                                       gTCCCCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAACTGAGAcc                                                                                                >  1:874873/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)
                                                                tCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGtagtag                                                                                       <  1:1237382/150‑1 (MQ=255)
                                                                 cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                                      >  1:1746345/1‑150 (MQ=255)
                                                                    aGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCt                                                                                   <  2:1234533/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)
                                                                       aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCTCgcgc                                                                                >  2:1450621/1‑150 (MQ=255)
                                                                       aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                                                >  1:922634/1‑150 (MQ=255)
                                                                       aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                                                >  1:1716419/1‑150 (MQ=255)
                                                                         cATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCgg                                                                              >  1:772347/1‑150 (MQ=255)
                                                                         cATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCgg                                                                              >  1:1722318/1‑150 (MQ=255)
                                                                            ttAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTGGCGCTCGCGGTTa                                                                           <  2:26945/150‑1 (MQ=255)
                                                                                 cTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCg                                                                      >  1:580642/1‑150 (MQ=255)
                                                                                   gacccTGACGGGGTCGTACTACCAGACTGGAATCCATCACCACCACCCGCTCCGGGCGGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTCACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTg                                                                    <  1:1668138/147‑1 (MQ=255)
                                                                                      ccTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAc                                                                 >  1:606887/1‑150 (MQ=255)
                                                                                       cTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAcc                                                                <  2:1780802/150‑1 (MQ=255)
                                                                                       cTGACGTGGTCGATCTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAcc                                                                <  2:734977/150‑1 (MQ=255)
                                                                                          aCGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATg                                                             <  2:149036/150‑1 (MQ=255)
                                                                                          aCGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATg                                                             <  2:694225/150‑1 (MQ=255)
                                                                                             tGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATa                                                          >  1:1057571/1‑150 (MQ=255)
                                                                                              gggcGTACTACCAGACGGGGATCCATCACCACCACCGGCGCCGTGCGGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                                         <  1:131110/147‑1 (MQ=255)
                                                                                              ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGTTCCGTGCTGACCTCGCCGCCGGTTGGGCTGCGTCTTTTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTTGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                                         <  2:1669634/150‑1 (MQ=255)
                                                                                              ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                                         <  1:526607/150‑1 (MQ=255)
                                                                                                tCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTa                                                       >  1:1035707/1‑150 (MQ=255)
                                                                                                  gTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTAGCGGTTATCCCGTGGACCATGATAGTaaa                                                     >  2:348088/1‑150 (MQ=255)
                                                                                                   tactacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTaaaa                                                    >  2:330988/1‑150 (MQ=255)
                                                                                                   tactacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTaaaa                                                    <  2:863994/150‑1 (MQ=255)
                                                                                                    actacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAt                                                   <  1:1373190/150‑1 (MQ=255)
                                                                                                     ctacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAtt                                                  >  1:1655983/1‑150 (MQ=255)
                                                                                                     ctacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAtt                                                  >  1:1108227/1‑150 (MQ=255)
                                                                                                       acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGCCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTa                                                >  1:418450/1‑150 (MQ=255)
                                                                                                        ccAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaa                                               >  1:405122/1‑150 (MQ=255)
                                                                                                         cAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaa                                              >  2:1717524/1‑150 (MQ=255)
                                                                                                            aaTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATa                                           <  2:1683254/150‑1 (MQ=255)
                                                                                                             aTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                                          >  2:400048/1‑150 (MQ=255)
                                                                                                              tGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACt                                         <  1:705145/150‑1 (MQ=255)
                                                                                                                   tCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGctc                                    >  1:991810/1‑150 (MQ=255)
                                                                                                                    ccctCACCACCACCGGCTCCGTGCTGACCAAGCCGCCGGTTGGGCTGCGTCTTGTATTGTTACGAGTGGGGACCTCGCGACATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGctct                                   <  2:303005/147‑1 (MQ=255)
                                                                                                                       tCACCACCACCTGCTCCTTTCTTGCCACGCCGCCGGTTGGGCTGCGTATTGTATGGTTACGAGTGGGGACCTCTCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCAc                                <  2:307825/150‑1 (MQ=255)
                                                                                                                        cagcccaaccGGCGGCGTGGTCAGCACGGAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCTTGATAGTTAATTTTAAAATTCTCGCGCTcccc                               >  1:284360/8‑150 (MQ=25)
                                                                                                                         accaccaccGTCTCCGTGCTGACCACGCCGCGTGTTGGGCTGCGTCTTGTATGGTTTCGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCa                              <  2:278401/150‑1 (MQ=255)
                                                                                                                         accaccaccGGCTCCGTGCTGACCACGCCGCGGGTTGGGGTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAATATTTAAAATACTCGCTCTCACCa                              <  2:1544188/150‑1 (MQ=255)
                                                                                                                           caccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGc                            >  1:637656/1‑150 (MQ=255)
                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTATAATTTAAAATACTCGCTCTCACCAGcc                           >  2:1127370/1‑150 (MQ=255)
                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc                           <  1:658667/150‑1 (MQ=255)
                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc                           >  1:1254380/1‑150 (MQ=255)
                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATAATAGTAAAATTTAAAATACTCGCTCTCACCAGcc                           >  2:710220/1‑150 (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)
                                                                                                                                ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCa                       >  2:54940/1‑150 (MQ=255)
                                                                                                                                ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCa                       <  2:243785/150‑1 (MQ=255)
                                                                                                                                  ggCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAAc                     <  1:523364/150‑1 (MQ=255)
                                                                                                                                     tCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGcc                  >  1:365265/1‑150 (MQ=255)
                                                                                                                                       ggTGCTGACCACGCCGCCGGTTGGGCTGTGTCTTGTATGGATACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAg                <  2:810493/149‑1 (MQ=255)
                                                                                                                                            tGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGc           >  1:936662/1‑150 (MQ=255)
                                                                                                                                              aCCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCgg         <  1:54940/150‑1 (MQ=255)
                                                                                                                                                 aCGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCggggg      <  1:1798326/150‑1 (MQ=255)
                                                                                                                                                  cgccgcCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCgggggg     <  2:540696/150‑1 (MQ=255)
                                                                                                                                                    ccgccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCg   <  1:617249/150‑1 (MQ=255)
                                                                                                                                                    ccgccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCg   <  2:245982/150‑1 (MQ=255)
                                                                                                                                                     cgccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGa  <  2:1466728/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGA  >  NC_000913/304657‑304955

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: