Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,790 T→G 6.2% E139A (GAG→GCG)  ecpE ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,7900TG 6.2% 435.4 / 12.7 151E139A (GAG→GCG) ecpEECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base T (78/62);  new base G (3/7);  total (82/69)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 153 total aligned reads displayed.

CTTCGGTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCC  >  NC_000913/304641‑304939
                                                                                                                                                     |                                                                                                                                                     
ctTCGGTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTGACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCt                                                                                                                                                       >  2:1080803/1‑150 (MQ=255)
    ggTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGt                                                                                                                                                   >  1:863994/1‑150 (MQ=255)
     gTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTg                                                                                                                                                  >  1:1466728/1‑150 (MQ=255)
         aaTCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGa                                                                                                                                              <  1:1768089/150‑1 (MQ=255)
          aTCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGAc                                                                                                                                             <  2:829656/150‑1 (MQ=255)
            cacaTCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCa                                                                                                                                           >  1:1004991/1‑150 (MQ=255)
             acaTCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACTTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCAc                                                                                                                                          >  1:1160169/1‑150 (MQ=255)
               aTCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCAcgc                                                                                                                                        >  1:418983/1‑150 (MQ=255)
                tCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCAcgcc                                                                                                                                       <  2:489192/150‑1 (MQ=255)
                 ccTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGACCCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGCCGTGCTGACCccgccg                                                                                                                                      >  2:131110/1‑150 (MQ=255)
                  cTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACAAGTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACAACCGGCTCCGTGCTGACCAcgccgc                                                                                                                                     >  1:49811/1‑150 (MQ=255)
                   tGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCCCCACCGGCTCCGTGCTGACCACgccgcc                                                                                                                                    >  2:929249/1‑150 (MQ=255)
                        tcaTCAGTAGCTTAAACCATGTGCTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACTTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtgg                                                                                                                               <  2:1090801/148‑3 (MQ=255)
                          aTCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCCCCTTGTCGAAGGACCATTTAAACTGAACCTGACGGGGTCGTACTACCAGATTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGGTggg                                                                                                                             >  1:1617913/1‑150 (MQ=255)
                          aTCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGCTCCTGTCACCTTGTCGAAGGACCATTTACACTGAACCTGACGTGGGCGTACTACCAGAATGGTATCCAGCACCACCACCGGCGGCGTGGTCAGCACGGAGCCGGtggtg                                                                                                                             <  1:398609/150‑5 (MQ=255)
                              gTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGc                                                                                                                         >  2:213093/1‑150 (MQ=255)
                                aGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGt                                                                                                                       >  1:467815/1‑150 (MQ=255)
                                aGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGt                                                                                                                       >  1:245982/1‑150 (MQ=255)
                                 gcataccccATGGGGTGCCAGGGTTACTTACCGTGCCTGTCACCTTGTCGAAGGACCATTTAAACGGAACCTGACTTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTc                                                                                                                      <  2:270292/143‑1 (MQ=255)
                                 gCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGta                                                                                                                      >  1:149036/1‑149 (MQ=255)
                                    tAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCGGCGTCTTg                                                                                                                   >  1:319152/1‑150 (MQ=255)
                                     aaaCCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGt                                                                                                                  <  2:119286/150‑1 (MQ=255)
                                              gtTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACg                                                                                                         >  2:1824850/1‑150 (MQ=255)
                                              gtTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTACGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACg                                                                                                         >  2:403062/1‑150 (MQ=255)
                                                   cAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTgg                                                                                                    <  1:1074519/150‑1 (MQ=255)
                                                   cAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTgg                                                                                                    >  2:1531189/1‑150 (MQ=255)
                                                   cAGTCTTAGTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCTTCTTGTATTGTTACGAGTgg                                                                                                    <  1:1089963/150‑1 (MQ=255)
                                                     gtgtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTgggg                                                                                                  <  2:30329/150‑1 (MQ=255)
                                                       gtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGAc                                                                                                <  2:153695/150‑1 (MQ=255)
                                                       gtTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGGCTTGTATGGTTACGAGTGGGGAc                                                                                                <  1:1822296/150‑1 (MQ=255)
                                                         tacttacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCt                                                                                              >  1:498260/1‑150 (MQ=255)
                                                             tacCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTATTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTcgcg                                                                                          >  1:1176993/1‑150 (MQ=255)
                                                              acCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGa                                                                                         <  2:1822916/150‑1 (MQ=255)
                                                                cGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGaaa                                                                                       >  1:1804277/1‑150 (MQ=255)
                                                                    ccTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATgag                                                                                   <  2:1391339/150‑1 (MQ=255)
                                                                      tGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGAc                                                                                 <  1:902381/150‑1 (MQ=255)
                                                                         cACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCAGCCGGTTGGGCTGCGGCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCg                                                                              >  2:1084088/1‑150 (MQ=255)
                                                                          aCCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCCCCACCACCGGCTCCGAGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCgg                                                                             >  1:99219/1‑150 (MQ=255)
                                                                           ccTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGTTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgt                                                                            <  2:1417159/150‑1 (MQ=255)
                                                                             ttGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCTGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGgtag                                                                          <  2:752044/150‑1 (MQ=255)
                                                                                 cGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAgc                                                                      >  1:1746345/1‑150 (MQ=255)
                                                                                    aGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCt                                                                   <  2:1234533/150‑1 (MQ=255)
                                                                                     ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGCCCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                  >  1:407883/1‑150 (MQ=255)
                                                                                     ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                  >  1:1678450/1‑150 (MQ=255)
                                                                                     ggACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                                  <  1:1571296/150‑1 (MQ=255)
                                                                                       aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCTCgcgc                                                                >  2:1450621/1‑150 (MQ=255)
                                                                                       aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                                >  1:1716419/1‑150 (MQ=255)
                                                                                       aCCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                                >  1:922634/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)
                                                                                               aaCTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTATCGCTCGCGGTTATcc                                                        >  1:483392/1‑150 (MQ=255)
                                                                                                 cTGAACCTGATGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCGGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAAGGAGACCCGGTCGGAGCGCTGGCGGTTTTCCCg                                                      >  2:1099992/1‑150 (MQ=255)
                                                                                                 cTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCg                                                      >  1:580642/1‑150 (MQ=255)
                                                                                                   gacccTGACGGGGTCGTACTACCAGACTGGAATCCATCACCACCACCCGCTCCGGGCGGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTCACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTg                                                    <  1:1668138/147‑1 (MQ=255)
                                                                                                      ccTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAc                                                 >  1:606887/1‑150 (MQ=255)
                                                                                                       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)
                                                                                                              gggcGTACTACCAGACGGGGATCCATCACCACCACCGGCGCCGTGCGGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                         <  1:131110/147‑1 (MQ=255)
                                                                                                              ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGTTCCGTGCTGACCTCGCCGCCGGTTGGGCTGCGTCTTTTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTTGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                         <  2:1669634/150‑1 (MQ=255)
                                                                                                              ggTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAg                                         >  1:316502/1‑150 (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)
                                                                                                                      tacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAttt                                 >  1:849377/1‑150 (MQ=255)
                                                                                                                       acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGCCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTa                                >  1:418450/1‑150 (MQ=255)
                                                                                                                       acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTa                                >  1:1246290/1‑150 (MQ=255)
                                                                                                                        ccAGAATGGTATGCATCTGCAGCACGGGCTCCGTGCTGACCACGCCGGCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGATTTGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATATTAAAATTTaa                               <  2:1784596/150‑1 (MQ=255)
                                                                                                                         cAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaa                              <  2:1717335/150‑1 (MQ=255)
                                                                                                                          aGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaaa                             <  1:1297283/150‑1 (MQ=255)
                                                                                                                            aaTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATa                           <  2:1683254/150‑1 (MQ=255)
                                                                                                                             aTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                          >  2:400048/1‑150 (MQ=255)
                                                                                                                             aTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGCCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                          >  2:374767/1‑150 (MQ=255)
                                                                                                                              tGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACt                         <  1:705145/150‑1 (MQ=255)
                                                                                                                                 tATCCATCACCACCTCCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCTGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGc                      <  1:1354976/150‑1 (MQ=255)
                                                                                                                                  aTCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGct                     <  1:621289/150‑1 (MQ=255)
                                                                                                                                    ccctCACCACCACCGGCTCCGTGCTGACCAAGCCGCCGGTTGGGCTGCGTCTTGTATTGTTACGAGTGGGGACCTCGCGACATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGctct                   <  2:303005/147‑1 (MQ=255)
                                                                                                                                        cagcccaaccGGCGGCGTGGTCAGCACGGAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCTTGATAGTTAATTTTAAAATTCTCGCGCTcccc               >  1:284360/8‑150 (MQ=25)
                                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTATAATTTAAAATACTCGCTCTCACCAGcc           >  2:1127370/1‑150 (MQ=255)
                                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc           >  1:1254380/1‑150 (MQ=255)
                                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc           <  1:671815/150‑1 (MQ=255)
                                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc           <  1:658667/150‑1 (MQ=255)
                                                                                                                                            accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATAATAGTAAAATTTAAAATACTCGCTCTCACCAGcc           >  2:710220/1‑150 (MQ=255)
                                                                                                                                             ccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCt          <  1:929249/150‑1 (MQ=255)
                                                                                                                                             ccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCt          >  2:1009381/1‑150 (MQ=255)
                                                                                                                                             caaccGGCGGCGTGGTCAGCACGGAGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGCGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGAAAAATTTAAAATACTCCCTCTCCCCTGCCt          >  1:1138819/3‑150 (MQ=255)
                                                                                                                                              caccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCAGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTATCACCAGCATg         >  2:644118/1‑150 (MQ=255)
                                                                                                                                               aacgtcTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTTTGGTTACGAGTGGGGACCTCGCGCAATGAGTCCCGGTAGTAGGGCGCGCGGTTCTCCCGTGGACCAAGATAGTAAAATTTAAAATACACGCTCTCCCCAGCCTGc        >  2:75017/6‑150 (MQ=255)
                                                                                                                                                ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCa       <  2:1225339/150‑1 (MQ=255)
                                                                                                                                                ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCa       <  2:243785/150‑1 (MQ=255)
                                                                                                                                                    cTCCGTGCTGAGCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGCGCCCCGGTAGTAGCGCTCGCGGTTATCCCGGGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGc   >  2:1647655/1‑150 (MQ=255)
                                                                                                                                                    cTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGc   <  2:498260/150‑1 (MQ=255)
                                                                                                                                                     tCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGcc  >  1:365265/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTCGGTCGAATCACATCCTGGCTTAATCAGTAGCTTAAACCATGTGTTGCCAGTGTTACTTACCGTTCCTGTCACCTTGTCGAAGGACCATTTAAACTGAACCTGACGTGGTCGTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTGCAACGCC  >  NC_000913/304641‑304939

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