Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,113 T→C 7.8% D25G (GAT→GGT)  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,1130TC 7.8% 391.8 / 14.2 123D25G (GAT→GGT) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base T (61/52);  new base C (6/4);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 124 total aligned reads displayed.

TGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACG  >  NC_000913/671968‑672261
                                                                                                                                                 |                                                                                                                                                     
   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                                                                                                                                      <  4:116540/149‑1 (MQ=255)
              gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                      >  1:898540/1‑149 (MQ=255)
                   gTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAccc                                                                                                                                 <  4:89560/149‑1 (MQ=255)
                   gTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAccc                                                                                                                                 >  1:327321/1‑149 (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                                                                                                                            >  3:399522/1‑149 (MQ=255)
                        tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                            >  4:311319/1‑149 (MQ=255)
                        tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                            <  1:621649/149‑1 (MQ=255)
                        tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                            >  1:350675/1‑149 (MQ=255)
                         ttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                           <  3:614008/149‑1 (MQ=255)
                         ttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                           <  1:149275/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)
                           aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACACACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                         >  4:657807/1‑149 (MQ=255)
                            tCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAt                                                                                                                        <  3:621451/149‑1 (MQ=255)
                               aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGCCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGGCACACCCGGCGGTGATTaa                                                                                                                     <  2:56348/149‑1 (MQ=255)
                               aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                     >  2:85771/1‑149 (MQ=255)
                               aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                     <  4:655708/149‑1 (MQ=255)
                               aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                     >  2:671151/1‑149 (MQ=255)
                                  aaCTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAacac                                                                                                                 >  2:381740/1‑149 (MQ=255)
                                     tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGc                                                                                                              >  1:115079/1‑149 (MQ=255)
                                        acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCag                                                                                                           >  1:222823/1‑149 (MQ=255)
                                         cgccATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCTCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTACGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCaga                                                                                                          <  2:170134/147‑1 (MQ=255)
                                          acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTGTTAATCACCGccgggt                                                                                                         <  4:639326/149‑5 (MQ=255)
                                          acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCagag                                                                                                         >  3:61177/1‑149 (MQ=255)
                                            cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGAt                                                                                                       >  4:476368/1‑149 (MQ=255)
                                              ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaa                                                                                                     <  3:311319/149‑1 (MQ=255)
                                              ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaa                                                                                                     <  3:370036/149‑1 (MQ=255)
                                               tCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaac                                                                                                    <  3:476368/149‑1 (MQ=255)
                                                 aGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCGAGAGATaacaa                                                                                                  <  3:515380/149‑1 (MQ=255)
                                                  gACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaac                                                                                                 <  4:14116/149‑1 (MQ=255)
                                                    cGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAAcaacaa                                                                                               >  4:40278/1‑149 (MQ=255)
                                                     gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAt                                                                                              <  1:159469/149‑1 (MQ=255)
                                                     gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAt                                                                                              <  1:561890/149‑1 (MQ=255)
                                                      tAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                                <  4:14843/99‑1 (MQ=255)
                                                      tAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                                >  3:14843/1‑99 (MQ=255)
                                                        aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGtt                                                                                           <  2:84070/149‑1 (MQ=255)
                                                        aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGtt                                                                                           <  4:576849/149‑1 (MQ=255)
                                                          tGGGTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGc                                                                                         >  1:766854/1‑149 (MQ=255)
                                                            gTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCa                                                                                       >  1:574590/1‑149 (MQ=255)
                                                                  caccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                                 <  3:330771/149‑1 (MQ=255)
                                                                   accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGACAGATATCg                                                                                >  3:496810/1‑149 (MQ=255)
                                                                       cacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaa                                                                                                  >  4:291532/1‑127 (MQ=255)
                                                                       cacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaa                                                                                                  <  3:291532/127‑1 (MQ=255)
                                                                           gCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                        >  1:137478/1‑149 (MQ=255)
                                                                             ttAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATATTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAATCACCGccggg                                                                                                          <  4:506925/113‑4 (MQ=255)
                                                                             ttAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAATCACCGccggg                                                                                                          >  3:506925/1‑110 (MQ=255)
                                                                              tAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcctc                                                                     <  4:584272/149‑1 (MQ=255)
                                                                                    cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCgcgc                                                               <  3:776263/149‑1 (MQ=255)
                                                                                            gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                       >  4:114643/1‑149 (MQ=255)
                                                                                             aTCGCCTTAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAAc                                                      >  2:477425/1‑149 (MQ=255)
                                                                                                   tGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAg                                                >  1:483792/1‑149 (MQ=255)
                                                                                                    gAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGa                                               >  1:231605/1‑149 (MQ=255)
                                                                                                    gAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGa                                               <  2:222823/149‑1 (MQ=255)
                                                                                                         cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                          >  2:658367/1‑149 (MQ=255)
                                                                                                         cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                          <  1:381740/149‑1 (MQ=255)
                                                                                                          aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                         >  4:801074/1‑149 (MQ=255)
                                                                                                            gATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgaggag                                       >  3:803706/1‑149 (MQ=255)
                                                                                                              tCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                     <  3:657807/149‑1 (MQ=255)
                                                                                                                 tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                                  >  3:680217/1‑149 (MQ=255)
                                                                                                                 tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                                  <  4:184902/149‑1 (MQ=255)
                                                                                                                 tGACCTTCATAGTGGAAGAAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                                  >  2:660988/1‑149 (MQ=255)
                                                                                                                  gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAcc                                 >  3:529674/1‑149 (MQ=255)
                                                                                                                   aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGtt                                                                                           >  4:482113/1‑90 (MQ=255)
                                                                                                                    ccTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTg                               >  1:482965/1‑149 (MQ=255)
                                                                                                                     cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGa                                                                                     <  4:375974/94‑1 (MQ=255)
                                                                                                                      ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGt                             >  3:771138/1‑149 (MQ=255)
                                                                                                                      ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGt                             >  4:31559/1‑149 (MQ=255)
                                                                                                                        cATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTAc                           >  3:405727/1‑149 (MQ=255)
                                                                                                                            gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTaa                       >  1:699698/1‑149 (MQ=255)
                                                                                                                            gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTaa                       >  1:362188/1‑149 (MQ=255)
                                                                                                                            gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAAATAAAAACAATGTTGCCAGATATCGCACGCTTCCTCCCACGCTTAGCGAACGGACAGATTTAGGAGTTTACCTGGTACGTaa                       >  3:463079/1‑149 (MQ=255)
                                                                                                                             tGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTaaa                      >  2:576854/1‑149 (MQ=255)
                                                                                                                              ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                     <  2:349977/149‑1 (MQ=255)
                                                                                                                                   gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                >  1:323424/1‑149 (MQ=255)
                                                                                                                                   gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                <  4:581391/149‑1 (MQ=255)
                                                                                                                                     atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa              >  3:419197/3‑149 (MQ=255)
                                                                                                                                     atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa              >  4:554436/3‑149 (MQ=255)
                                                                                                                                      tcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTAc             >  3:747775/2‑149 (MQ=255)
                                                                                                                                      ccTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTAc             <  4:602784/149‑1 (MQ=255)
                                                                                                                                       cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                            <  1:790896/101‑1 (MQ=255)
                                                                                                                                       cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                            >  2:790896/1‑101 (MQ=255)
                                                                                                                                       cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACg            >  1:410003/1‑149 (MQ=255)
                                                                                                                                         cgttGGTACGACCCAAATGCCAGCAATAACCGGCGGTGATAACCACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGATTAGCCAACGACCAGATTGAGTAGTTTACCTGGTACGTAAAGTGATTTACGTa          <  2:461837/146‑1 (MQ=255)
                                                                                                                                               tATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTaa    >  1:775592/1‑149 (MQ=255)
                                                                                                                                               tATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTaa    <  2:50793/149‑1 (MQ=255)
                                                                                                                                               tATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTaa    >  3:714662/1‑149 (MQ=255)
                                                                                                                                               tATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTaa    >  4:814475/1‑149 (MQ=255)
                                                                                                                                                aTCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc   >  4:563182/1‑149 (MQ=255)
                                                                                                                                                aTCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc   >  1:516465/1‑149 (MQ=255)
                                                                                                                                                 tCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACg  >  1:485289/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACG  >  NC_000913/671968‑672261

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: