Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,110 G→A 5.4% 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 5.4% 800.1 / 23.9 243T26I (ACC→ATC) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base G (111/117);  new base A (6/7);  total (119/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-02

Only 100 of 244 total aligned reads displayed.

CGGCCGGGGATCAACACGGTCGCATTAACCGTCATGATCATCTGATACTCTGCTGTTTGACCGTTACGGAATACCGATGCGGTATCTTTCGCGATGCTCACTTTACCCAAACGCAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCT  >  NC_000913/671849‑672370
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                    
cGGCCGGGGTTCAAAATGGTTGCATTAAACGTCATGTTCATCTGATTCTCCGCTGCTTGACCTTTACGGCATCCCCATGCTGTATCTTTCGTTATGCTACCTTTAACCAAACGCAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAGGCAAGTCGACACCATTCAGAC‑GTAACTGGTTACGCTCCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATg                                                                                                                                                                                                                                                         <  4:66207/275‑1 (MQ=255)
                                    atcatcTGATACTCTGCTGTTTGACCGTTACGGAATACCGATGCGGTATCTTTCGCGATGCTCACTTTACCCAAACGCAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCCGGATCCTGACCTTCATAGTGGAAGGAACCTGCGTGTTATCACTCAGATGCCAGCCACACACGGCGGTGATTAACACCGcacgag                                                                                                                                                                                                                      >  1:10518/1‑269 (MQ=255)
                                         cTGATACTCTGCTGTTTGACCGTTACGGAATACCGATGCGGTATCTTTCGCGATGCTCACTTTACCCAAACGCAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATAACCTTCATAGGGGAAGGAACCTGCGGGGTATCACGGAGATGCCACCCACACCCGGAGGTGATTAACACCCCCCGAGaaaaac                                                                                                                                                                                                                >  4:200212/1‑273 (MQ=255)
                                                      gTTTGACCGTTACGGAATACCGATGCGGTATCTTTCGCGATGCTCACTTTACCCAAACGCAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGAACATGACCTTCATAGTGGAAGGAAACTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAAAACCTCCAGACATAACAACAATGTTGCCa                                                                                                                                                                                                    >  4:125030/1‑274 (MQ=255)
                                                                                        tCGCGATGCTCACTTTACCCAAACGCAAGGATGGAACGTACTTACTCGTGGTTTCTTGAGCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACATCAGTGTTGCCAGATATCGCACGATTCCTCCCGCGCTTAGCCAACGa                                                                                                                                                                 <  2:53019/275‑1 (MQ=255)
                                                                                                tctaCTTTAACCAAACGCAAGGATGGAACGTCCTTACGCGTGGTTTATTTTTCAAGCAACTCCAAACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTg                                                                                                                                                         <  4:10578/272‑1 (MQ=255)
                                                                                                ggcACTGTAGCCAAAGTAAAGGGTGGAAAGTCCTTACTGGTGGTTTCTTTATCAAGCAAATCGACAGCCGTC‑GACAGCAACTGGTTAATGACCGTACGGCTTAAGGGCCCGTTCGGATAGCCTGAGTCCAGGATCATGAGCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTg                                                                                                                                                         <  4:323495/273‑1 (MQ=255)
                                                                                                             aaCGCAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATACCAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCCCTTAGCCAACGACCAGCTTGAGGAGTTTACCTg                                                                                                                                            >  4:173559/1‑275 (MQ=255)
                                                                                                                 cAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCa                                                                                                                                                                                                                         <  2:102352/194‑1 (MQ=255)
                                                                                                                  aaGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTCCGTGGt                                                                                                                                                                                                                                                                     <  6:1045535/149‑1 (MQ=255)
                                                                                                                    ggATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTAt                                                                                                                                                                                                                                                                   <  8:1014405/149‑1 (MQ=255)
                                                                                                                     gATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATc                                                                                                                                                                                                                                                                  <  5:1342620/149‑1 (MQ=255)
                                                                                                                          aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCACACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCa                                                                                                                                                                                                                                                             <  7:743992/149‑1 (MQ=255)
                                                                                                                           aCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAg                                                                                                                                                                                                                                                            >  5:602001/1‑149 (MQ=255)
                                                                                                                             gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                                                                                                                                          <  8:276180/149‑1 (MQ=255)
                                                                                                                             gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                                                                                                                                          <  6:814432/149‑1 (MQ=255)
                                                                                                                              tCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCGAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATg                                                                                                                                                                                                                                                         >  8:186467/1‑149 (MQ=255)
                                                                                                                                cTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcc                                                                                                                                                                                                                                                       >  5:394101/1‑149 (MQ=255)
                                                                                                                                 ttACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCCCGCAGATgcca                                                                                                                                                                                                                                                      <  6:587183/149‑1 (MQ=255)
                                                                                                                                     gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                                                                                                                                  >  7:24069/1‑149 (MQ=255)
                                                                                                                                     gcgTGGTTTCTTCATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCTCCTGAGTCCATGATCATGACCTTCATAGTGGACGGAACCTGCGTGGTCTCACGCAGATgccaacca                                                                                                                                                                                                                                                  >  8:915762/1‑149 (MQ=255)
                                                                                                                                         ggTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAGGGTGCCAGATCACGCACGCTTCCTCCCGCCCTTAGCCAACGACCAGATTGAGGCGTTTANNNGGTACGTAAATAACTTTaactacctgccc                                                                                                                >  3:73798/1‑264 (MQ=255)
                                                                                                                                          gTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAccc                                                                                                                                                                                                                                             >  6:1152181/1‑149 (MQ=255)
                                                                                                                                           tttCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcg                                                                                                                                                                                                                                            >  5:679976/1‑149 (MQ=255)
                                                                                                                                            ttCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcgg                                                                                                                                                                                                                                           <  8:1326248/149‑1 (MQ=255)
                                                                                                                                              cTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACt                                                                                                                              <  2:120565/256‑1 (MQ=255)
                                                                                                                                                   tCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAt                                                                                                                                                                                                                                    >  6:313581/1‑149 (MQ=255)
                                                                                                                                                    cAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                                                                                                                                   <  5:52941/149‑1 (MQ=255)
                                                                                                                                                     aaGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTa                                                                                                                                                                                                                                  <  5:1010534/149‑1 (MQ=255)
                                                                                                                                                      aGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTgttaat                                                                                                                                                                                                                                 >  8:1181570/1‑144 (MQ=255)
                                                                                                                                                      aGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                                                                                                                                 >  5:315384/1‑149 (MQ=255)
                                                                                                                                                       gCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                                                                                                                                        <  5:121350/141‑1 (MQ=255)
                                                                                                                                                       gCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAac                                                                                                                                                                                                                                <  7:712995/149‑1 (MQ=255)
                                                                                                                                                         aaCTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAacac                                                                                                                                                                                                                              <  8:635335/149‑1 (MQ=255)
                                                                                                                                                         aaCTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAacac                                                                                                                                                                                                                              >  7:822771/1‑149 (MQ=255)
                                                                                                                                                           cTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCg                                                                                                                                                                                                                            >  7:459513/1‑149 (MQ=255)
                                                                                                                                                            tCGACACCATTCAGAC‑GTAACTGGTTACGGACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGc                                                                                                                                                                                                                           >  7:1251585/1‑149 (MQ=255)
                                                                                                                                                            tCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGc                                                                                                                                                                                                                           >  6:1039411/1‑149 (MQ=255)
                                                                                                                                                            tCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGc                                                                                                                                                                                                                           >  7:783087/1‑149 (MQ=255)
                                                                                                                                                            tCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGc                                                                                                                                                                                                                           >  8:431812/1‑149 (MQ=255)
                                                                                                                                                              gACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCa                                                                                                                                                                                                                         >  7:186241/1‑149 (MQ=255)
                                                                                                                                                               acacCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCag                                                                                                                                                                                                                        >  8:626016/1‑149 (MQ=255)
                                                                                                                                                                cacCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCaga                                                                                                                                                                                                                       <  7:1018722/149‑1 (MQ=255)
                                                                                                                                                                 acCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTGTTAATCACCGccgggt                                                                                                                                                                                                                      <  7:1181570/149‑5 (MQ=255)
                                                                                                                                                                 acCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                                                                                                                                   <  7:1168525/136‑1 (MQ=255)
                                                                                                                                                                 acCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCagag                                                                                                                                                                                                                      >  8:828653/1‑149 (MQ=255)
                                                                                                                                                                  ccATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAgaga                                                                                                                                                                                                                     <  8:1251585/149‑1 (MQ=255)
                                                                                                                                                                   cATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGAt                                                                                                                                                                                                                    <  6:1127534/149‑1 (MQ=255)
                                                                                                                                                                     ttCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAaca                                                                                                                                                                                                                               >  5:875655/1‑136 (MQ=255)
                                                                                                                                                                     ttCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAaca                                                                                                                                                                                                                               <  6:875655/136‑1 (MQ=255)
                                                                                                                                                                     ttCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaa                                                                                                                                                                                                                  >  5:506763/1‑149 (MQ=255)
                                                                                                                                                                      tCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaac                                                                                                                                                                                                                 <  6:983414/149‑1 (MQ=255)
                                                                                                                                                                         gAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaac                                                                                                                                                                                                              <  5:61806/149‑1 (MQ=255)
                                                                                                                                                                             gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAt                                                                                                                                                                                                           >  8:1377105/1‑149 (MQ=255)
                                                                                                                                                                             gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAt                                                                                                                                                                                                           >  6:1081649/1‑149 (MQ=255)
                                                                                                                                                                               aaCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAgg                                                                                                                                                                                                                                                           >  8:498333/1‑98 (MQ=255)
                                                                                                                                                                                aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATA‑CACGCAggtt                                                                                                                                                                                                                                                         <  6:1432379/99‑4 (MQ=255)
                                                                                                                                                                                  tGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGc                                                                                                                                                                                                      <  7:946482/149‑1 (MQ=255)
                                                                                                                                                                                  tGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGc                                                                                                                                                                                                      <  8:223964/149‑1 (MQ=255)
                                                                                                                                                                                  tGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGc                                                                                                                                                                                                      >  5:1063839/1‑149 (MQ=255)
                                                                                                                                                                                   ggTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGcc                                                                                                                                                                                                     <  6:425549/149‑1 (MQ=255)
                                                                                                                                                                                    gTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc                                                                                                                <  4:248364/233‑1 (MQ=255)
                                                                                                                                                                                                    cTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATc                                                                                                                                                                                                                                                                  >  7:387595/1‑71 (MQ=255)
                                                                                                                                                                                                      tAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcctc                                                                                                                                                                                  <  8:897955/149‑1 (MQ=255)
                                                                                                                                                                                                        aTGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACTCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTccc                                                                                                                                                                                >  7:138348/1‑149 (MQ=255)
                                                                                                                                                                                                           gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgcg                                                                                                                                                                             >  7:67676/1‑149 (MQ=255)
                                                                                                                                                                                                            cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGACATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCATTGc                                              >  2:226630/1‑275 (MQ=255)
                                                                                                                                                                                                             ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGa                                                                                                                                                                                                  >  5:1049346/1‑126 (MQ=255)
                                                                                                                                                                                                             ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCt                                                                                                                                                                           <  8:129606/149‑1 (MQ=255)
                                                                                                                                                                                                               gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                                                                                                                                       <  8:1201897/87‑1 (MQ=255)
                                                                                                                                                                                                               gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                                                                                                                                       >  7:1201897/1‑87 (MQ=255)
                                                                                                                                                                                                               gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                                                                                                                                         <  5:1273840/149‑1 (MQ=255)
                                                                                                                                                                                                                  cgtaTCGCCTGAGTCCAGGATCATGACCTTCATAGTTGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGa                                                                                       <  2:199628/225‑1 (MQ=255)
                                                                                                                                                                                                                   ggATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCa                                                                                                                                                                     >  6:1111268/1‑149 (MQ=255)
                                                                                                                                                                                                                        gCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGAc                                                                                                                                                                >  7:867078/1‑149 (MQ=255)
                                                                                                                                                                                                                          cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                                                                                                                                              >  7:59610/1‑149 (MQ=255)
                                                                                                                                                                                                                              gTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGAtt                                                                                                                                                          <  8:1106055/149‑1 (MQ=255)
                                                                                                                                                                                                                               tCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGtt                                                                                                                                                                                                        <  7:1318287/102‑1 (MQ=255)
                                                                                                                                                                                                                               tCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTg                                                                                                                                                         >  5:1322317/1‑149 (MQ=255)
                                                                                                                                                                                                                                ccAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAATGACCAGATTga                                                                                                                                                        <  5:1097680/149‑1 (MQ=255)
                                                                                                                                                                                                                                ccAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTga                                                                                                                                                        <  6:1271718/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                                                                                                                                       <  7:1359476/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                                                                                                                                       <  8:799673/149‑1 (MQ=255)
                                                                                                                                                                                                                                  aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATc                                                                                                                                                                                                                                                                  >  5:728013/1‑41 (MQ=255)
                                                                                                                                                                                                                                  aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                                                                                                                                      >  7:1069142/1‑149 (MQ=255)
                                                                                                                                                                                                                                  aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                                                                                                                                      <  8:568270/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                                                                                                                                     >  7:268233/1‑149 (MQ=255)
                                                                                                                                                                                                                                   ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                                                                                                                                     <  5:243991/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGttt                                                                                                                                                 >  8:668992/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aTGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTga                                                                                                                                                        >  6:487123/1‑141 (MQ=255)
                                                                                                                                                                                                                                         tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                                                                                                                                               <  8:487211/149‑1 (MQ=255)
                                                                                                                                                                                                                                         tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                                                                                                                                               <  6:679976/149‑1 (MQ=255)
                                                                                                                                                                                                                                         tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                                                                                                                                               >  5:1410153/1‑149 (MQ=255)
                                                                                                                                                                                                                                             cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg                                                                                                                                           <  8:337592/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGt                                                                                                                                          <  2:155081/143‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                                                                                                                                  <  5:1472476/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       gaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATc                                                                                                                                 <  8:763433/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                                                                                                                             <  7:632750/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACttt                                                                                                                            <  7:533221/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             atcTGCGTGATACCTCGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                                                                                                           >  7:110736/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                    
CGGCCGGGGATCAACACGGTCGCATTAACCGTCATGATCATCTGATACTCTGCTGTTTGACCGTTACGGAATACCGATGCGGTATCTTTCGCGATGCTCACTTTACCCAAACGCAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAC‑GTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCT  >  NC_000913/671849‑672370

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