Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,143 A→T 5.1% I15N (ATC→AAC)  lptE ← LPS assembly OM complex LptDE, lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913672,1430AT 5.1% 261.3 / 4.8 81I15N (ATC→AAC) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base T (2/2);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

ATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCC  >  NC_000913/671995‑672291
                                                                                                                                                    |                                                                                                                                                     
aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                                                       >  2:304430/1‑149 (MQ=255)
aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                                                       <  1:48063/149‑1 (MQ=255)
 tCAAGCCACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAt                                                                                                                                                      <  1:80193/149‑1 (MQ=255)
   aaGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTa                                                                                                                                                    >  3:82564/1‑149 (MQ=255)
    aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                                                   >  1:506359/1‑149 (MQ=255)
    aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                                                   >  2:630336/1‑149 (MQ=255)
    aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                                                   <  1:125494/149‑1 (MQ=255)
     gCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑c                                                                                                                                                 <  3:234250/149‑1 (MQ=255)
     gCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑c                                                                                                                                                 >  4:443319/1‑149 (MQ=255)
      cAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑ca                                                                                                                                                >  3:275648/1‑149 (MQ=255)
             acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAAT‑CACCGccgg                                                                                                                                         <  2:33850/149‑3 (MQ=255)
             acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCag                                                                                                                                         >  1:508223/1‑149 (MQ=255)
             acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTG‑TTAATCACCGccgg                                                                                                                                         <  3:593981/149‑3 (MQ=255)
              cacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCaga                                                                                                                                        <  3:443319/149‑1 (MQ=255)
               acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCagag                                                                                                                                       >  2:593829/1‑149 (MQ=255)
               acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTG‑TTAATCACCGccgggt                                                                                                                                       <  1:266077/149‑5 (MQ=255)
                 cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGAt                                                                                                                                     >  3:656981/1‑149 (MQ=255)
                  aTTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTGTTAAT‑CACCGccgggt                                                                                                                                       >  4:251721/1‑142 (MQ=255)
                  aTTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTG‑TTAATCACCGccgggt                                                                                                                                       <  3:251721/146‑5 (MQ=255)
                  aTTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATa                                                                                                                                    >  2:81449/1‑149 (MQ=255)
                   ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaa                                                                                                                                   <  1:168157/149‑1 (MQ=255)
                         cGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAAcaacaa                                                                                                                             <  4:86958/149‑1 (MQ=255)
                         cGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAAcaacaa                                                                                                                             <  2:235648/149‑1 (MQ=255)
                         cGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAAcaacaa                                                                                                                             <  1:397307/149‑1 (MQ=255)
                             aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGTAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                                                         <  2:601755/149‑1 (MQ=255)
                             aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                                                         <  4:331956/149‑1 (MQ=255)
                                       caccgcacGGCTTAATGGCCCGTTCGGATCGTCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                                                               <  4:103329/149‑1 (MQ=255)
                                        accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCg                                                                                                              <  1:625446/149‑1 (MQ=255)
                                        accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCg                                                                                                              >  3:123149/1‑149 (MQ=255)
                                           gcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCAc                                                                                                           >  1:586662/1‑149 (MQ=255)
                                            cacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACg                                                                                                          <  1:3733/149‑1 (MQ=255)
                                                     aTGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTccc                                                                                                 <  4:335602/149‑1 (MQ=255)
                                                      tGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcg                                                                                                >  3:404992/1‑149 (MQ=255)
                                                        gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgcg                                                                                              >  2:38464/1‑149 (MQ=255)
                                                        gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgcg                                                                                              >  3:436509/1‑149 (MQ=255)
                                                         cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCgcgc                                                                                             <  2:337897/149‑1 (MQ=255)
                                                          ccGTTCGGATCGCCTGAGTCCAGGCTCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGCTATCGCACGCTTCGTCCCGCGCt                                                                                            >  3:575289/1‑149 (MQ=255)
                                                           cGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCtt                                                                                           <  2:135148/149‑1 (MQ=255)
                                                            gTTCGGATCGCCTGAGTCCAGGATCATGGCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                                                          <  4:329042/149‑1 (MQ=255)
                                                              tCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGc                                                                                        >  3:281699/1‑149 (MQ=255)
                                                                 gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                                                     >  1:98214/1‑149 (MQ=255)
                                                                 gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                                                     >  2:617624/1‑149 (MQ=255)
                                                                       cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                                                               >  1:645958/1‑149 (MQ=255)
                                                                            tCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTg                                                                          >  3:384532/1‑149 (MQ=255)
                                                                              cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                                                        <  4:656981/149‑1 (MQ=255)
                                                                              cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                                                        <  2:586662/149‑1 (MQ=255)
                                                                                ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                                                      <  1:62534/149‑1 (MQ=255)
                                                                                ggATCATGACCTTCATAGTGGAACGAACCTGCGTGGTACAACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACTACAATGTTGCCAGATATCGCACGCTACCTCCCGCGCTTAGCCAACGACCAGATTgagga                                                                      <  3:248024/149‑1 (MQ=255)
                                                                                 gATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgaggag                                                                     <  1:630336/149‑1 (MQ=255)
                                                                                 gATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgaggag                                                                     >  3:335142/1‑149 (MQ=255)
                                                                                    cATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGttt                                                                  <  4:288399/149‑1 (MQ=255)
                                                                                      tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                                                                <  1:142848/149‑1 (MQ=255)
                                                                                       gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGCCCAGATTCAGGAGTCTAcc                                                               <  4:431572/149‑1 (MQ=255)
                                                                                        aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCt                                                              <  4:351658/149‑1 (MQ=255)
                                                                                         ccTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTg                                                             <  3:667666/149‑1 (MQ=255)
                                                                                          cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg                                                            >  1:105227/1‑149 (MQ=255)
                                                                                           ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGt                                                           >  4:343247/1‑149 (MQ=255)
                                                                                            tCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCGCGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGt                                                                    >  2:424131/1‑139 (MQ=255)
                                                                                            tCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGt                                                                    <  1:424131/139‑1 (MQ=255)
                                                                                                   ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGGGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                                                   >  3:608851/1‑149 (MQ=255)
                                                                                                    gtaggaaCCTGCGTGGTATCCCGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATATTGCCAGATATCACACGCTTCCTCCCGCGCTTAGCCAACGACCAGACTGAGGAGCTTACCTGGTACGTAAATc                                                  >  4:468399/3‑149 (MQ=255)
                                                                                                     aaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCa                                                 <  4:281699/149‑1 (MQ=255)
                                                                                                        gaaACTGCGTGGTATCACGCAGATGCCAGCCACACCCGGGGGTGATTAG‑CACCGCGAGAGATAACAGCAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTATCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                                              <  4:621842/149‑1 (MQ=255)
                                                                                                          atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                            >  3:357591/3‑149 (MQ=255)
                                                                                                             tGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGAcc                                                                                >  3:170244/1‑110 (MQ=255)
                                                                                                             tGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGAcc                                                                                <  4:170244/110‑1 (MQ=255)
                                                                                                              gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTa                                        <  4:82564/149‑1 (MQ=255)
                                                                                                               cGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAc                                       <  2:464952/149‑1 (MQ=255)
                                                                                                                 tGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAg                                     <  4:123149/149‑1 (MQ=255)
                                                                                                                 tGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAg                                     >  1:488097/1‑149 (MQ=255)
                                                                                                                      tCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                                                               >  4:207957/1‑102 (MQ=255)
                                                                                                                      tCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                                                               <  3:207957/102‑1 (MQ=255)
                                                                                                                      tCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACg                                >  1:391031/1‑149 (MQ=255)
                                                                                                                        aCGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGcc                              >  3:601365/1‑149 (MQ=255)
                                                                                                                                 ccagccaCACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGata                     <  4:615521/149‑1 (MQ=255)
                                                                                                                                    gccaCACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATAttt                  <  4:436509/149‑1 (MQ=255)
                                                                                                                                    gccaCACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATAttt                  <  3:303229/149‑1 (MQ=255)
                                                                                                                                     ccaCACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATAtttt                 <  2:391031/149‑1 (MQ=255)
                                                                                                                                        cacCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCt              <  2:98214/149‑1 (MQ=255)
                                                                                                                                         acCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTa             >  3:621710/1‑149 (MQ=255)
                                                                                                                                           ctggcggTGTTAAT‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc                                 <  4:600359/125‑1 (MQ=255)
                                                                                                                                           ctggcggTGTTAAT‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc                                 >  3:600359/3‑127 (MQ=255)
                                                                                                                                           ccGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTAcc           <  2:119076/149‑1 (MQ=255)
                                                                                                                                                    aTTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTcc  >  2:96878/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCC  >  NC_000913/671995‑672291

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑