Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,110 G→A 7.7% T26I (ACC→ATC)  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,1100GA 7.7% 369.0 / 15.4 117T26I (ACC→ATC) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base G (56/52);  new base A (6/3);  total (62/55)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 117 total aligned reads displayed.

AGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAG  >  NC_000913/671964‑672256
                                                                                                                                                  |                                                                                                                                                   
aGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTa                                                                                                                                                   <  1:742670/149‑1 (MQ=255)
    tGGAACGTCCTTATGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCAc                                                                                                                                               <  3:749144/149‑1 (MQ=255)
    tGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCAc                                                                                                                                               >  4:567182/1‑149 (MQ=255)
      gAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGc                                                                                                                                             >  4:330771/1‑149 (MQ=255)
       aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                            <  1:160691/149‑1 (MQ=255)
       aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                            <  2:711281/149‑1 (MQ=255)
        aCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                           <  4:306980/149‑1 (MQ=255)
              ttACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                     >  3:581391/1‑149 (MQ=255)
                  gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                 >  3:843445/1‑149 (MQ=255)
                  gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                 >  1:898540/1‑149 (MQ=255)
                  gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                 <  4:116540/149‑1 (MQ=255)
                       gTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAccc                                                                                                                            >  1:327321/1‑149 (MQ=255)
                       gTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAccc                                                                                                                            <  4:89560/149‑1 (MQ=255)
                        tttCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcg                                                                                                                           >  3:75972/1‑149 (MQ=255)
                        tttCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAggt                                                                                                                                         >  2:626387/1‑133 (MQ=255)
                        tttCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAggt                                                                                                                                         <  1:626387/135‑3 (MQ=255)
                          tCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACTGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCGTCATAGTGGAAGGAACCTGCGTGGTATCACGCGGATGCCAGCCACACCcggc                                                                                                                         <  4:555751/149‑1 (MQ=255)
                            tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                       >  4:311319/1‑149 (MQ=255)
                            tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                       <  1:621649/149‑1 (MQ=255)
                            tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                       >  1:350675/1‑149 (MQ=255)
                            tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                       >  3:399522/1‑149 (MQ=255)
                             ttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                      <  1:149275/149‑1 (MQ=255)
                             ttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                      <  3:614008/149‑1 (MQ=255)
                              tagcAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTg                                                                                                                     >  4:856027/4‑149 (MQ=255)
                              tATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTg                                                                                                                     <  2:373667/149‑1 (MQ=255)
                               aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                    <  4:679337/149‑1 (MQ=255)
                               aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                    <  2:322086/149‑1 (MQ=255)
                                tCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAt                                                                                                                   <  3:621451/149‑1 (MQ=255)
                                   aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGCCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGGCACACCCGGCGGTGATTaa                                                                                                                <  2:56348/149‑1 (MQ=255)
                                   aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                >  2:671151/1‑149 (MQ=255)
                                   aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                <  4:655708/149‑1 (MQ=255)
                                   aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                >  2:85771/1‑149 (MQ=255)
                                     cAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑ca                                                                                                             <  1:671151/149‑1 (MQ=255)
                                      aaCTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑cac                                                                                                            >  2:381740/1‑149 (MQ=255)
                                         tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAAT‑CACCGc                                                                                                         <  3:276497/149‑1 (MQ=255)
                                         tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGc                                                                                                         >  1:115079/1‑149 (MQ=255)
                                             cgccATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCTCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTACGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCaga                                                                                                     <  2:170134/147‑1 (MQ=255)
                                              acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTG‑TTAATCACCGccgggt                                                                                                    <  4:639326/149‑5 (MQ=255)
                                              acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCagag                                                                                                    >  3:61177/1‑149 (MQ=255)
                                              acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCagag                                                                                                    >  4:608952/1‑149 (MQ=255)
                                                cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGAt                                                                                                  >  4:476368/1‑149 (MQ=255)
                                                  ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaa                                                                                                <  3:370036/149‑1 (MQ=255)
                                                   tCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaac                                                                                               <  3:476368/149‑1 (MQ=255)
                                                     aGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCGAGAGATaacaa                                                                                             <  3:515380/149‑1 (MQ=255)
                                                      gACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaacaac                                                                                            <  4:14116/149‑1 (MQ=255)
                                                         gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAAt                                                                                         <  1:561890/149‑1 (MQ=255)
                                                         gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAAt                                                                                         <  1:159469/149‑1 (MQ=255)
                                                          tAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                           >  3:14843/1‑99 (MQ=255)
                                                          tAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                           <  4:14843/99‑1 (MQ=255)
                                                            aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                      <  2:84070/149‑1 (MQ=255)
                                                            aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                      <  1:203078/149‑1 (MQ=255)
                                                              tGGGTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGc                                                                                    >  1:766854/1‑149 (MQ=255)
                                                                      caccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                            <  3:330771/149‑1 (MQ=255)
                                                                       accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGACAGATATCg                                                                           >  3:496810/1‑149 (MQ=255)
                                                                         cgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCa                                                                         <  1:471130/149‑1 (MQ=255)
                                                                          gcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCAc                                                                        >  1:349977/1‑149 (MQ=255)
                                                                           cacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaacaa                                                                                             <  3:291532/127‑1 (MQ=255)
                                                                           cacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaacaa                                                                                             >  4:291532/1‑127 (MQ=255)
                                                                               gCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                   >  1:137478/1‑149 (MQ=255)
                                                                               gCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                   >  1:850350/1‑149 (MQ=255)
                                                                                cTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATACCGCACGCTtcc                                                                  <  4:75972/149‑1 (MQ=255)
                                                                                 ttAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATATTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAAT‑CACCGccggg                                                                                                     <  4:506925/113‑4 (MQ=255)
                                                                                 ttAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAAT‑CACCGccggg                                                                                                     >  3:506925/1‑110 (MQ=255)
                                                                                  tAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcctc                                                                <  4:584272/149‑1 (MQ=255)
                                                                                        cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCgcgc                                                          <  3:776263/149‑1 (MQ=255)
                                                                                                gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                  >  4:114643/1‑149 (MQ=255)
                                                                                                 aTCGCCTTAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAAc                                                 >  2:477425/1‑149 (MQ=255)
                                                                                                       tGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAg                                           >  1:483792/1‑149 (MQ=255)
                                                                                                        gAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGa                                          >  1:231605/1‑149 (MQ=255)
                                                                                                        gAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGa                                          <  2:222823/149‑1 (MQ=255)
                                                                                                             cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                     >  1:52653/1‑149 (MQ=255)
                                                                                                             cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                     >  2:658367/1‑149 (MQ=255)
                                                                                                              aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                    >  4:801074/1‑149 (MQ=255)
                                                                                                                     tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                             >  3:680217/1‑149 (MQ=255)
                                                                                                                     tGACCTTCATAGTGGAAGAAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                             >  2:660988/1‑149 (MQ=255)
                                                                                                                       aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                      <  3:482113/90‑1 (MQ=255)
                                                                                                                       aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                      >  4:482113/1‑90 (MQ=255)
                                                                                                                        ccTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTg                          >  1:482965/1‑149 (MQ=255)
                                                                                                                         cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGa                                                                                <  4:375974/94‑1 (MQ=255)
                                                                                                                         cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGa                                                                                >  3:375974/1‑94 (MQ=255)
                                                                                                                          ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGt                        >  4:31559/1‑149 (MQ=255)
                                                                                                                          ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGt                        >  3:771138/1‑149 (MQ=255)
                                                                                                                            cATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTAc                      >  3:405727/1‑149 (MQ=255)
                                                                                                                                gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTaa                  >  1:362188/1‑149 (MQ=255)
                                                                                                                                gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTaa                  >  1:699698/1‑149 (MQ=255)
                                                                                                                                gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAAATAAAAACAATGTTGCCAGATATCGCACGCTTCCTCCCACGCTTAGCGAACGGACAGATTTAGGAGTTTACCTGGTACGTaa                  >  3:463079/1‑149 (MQ=255)
                                                                                                                                  ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                <  2:349977/149‑1 (MQ=255)
                                                                                                                                       gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt           <  4:149986/149‑1 (MQ=255)
                                                                                                                                       gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt           >  1:323424/1‑149 (MQ=255)
                                                                                                                                       gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt           <  4:581391/149‑1 (MQ=255)
                                                                                                                                         atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa         >  4:554436/3‑149 (MQ=255)
                                                                                                                                         aCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa         >  3:636999/1‑149 (MQ=255)
                                                                                                                                         aCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa         >  2:646100/1‑149 (MQ=255)
                                                                                                                                          tcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTAc        >  3:747775/2‑149 (MQ=255)
                                                                                                                                          ccTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTAc        <  4:602784/149‑1 (MQ=255)
                                                                                                                                           cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                       >  2:790896/1‑101 (MQ=255)
                                                                                                                                           cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACg       >  1:410003/1‑149 (MQ=255)
                                                                                                                                             cgttGGTACGACCCAAATGCCAGCAATAACCGGCGGTGATAAC‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGATTAGCCAACGACCAGATTGAGTAGTTTACCTGGTACGTAAAGTGATTTACGTa     <  2:461837/146‑1 (MQ=255)
                                                                                                                                               gTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACa   >  2:488284/1‑149 (MQ=255)
                                                                                                                                                tGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAg  <  3:567182/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAG  >  NC_000913/671964‑672256

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