Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,486,869 T→A L218M (TTG→ATG)  lptF → lipopolysaccharide transport system protein LptF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,486,8690TA97.3% 249.7 / NA 75L218M (TTG→ATG) lptFlipopolysaccharide transport system protein LptF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (40/33);  minor base C (2/0);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGCGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGTCGAT  >  NC_000913/4486624‑4487113
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
aagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGGGACGGCGCCCAGGTCGTCACTCTCAGCCAGGGAACGCGCTTCGAAGGCACTGCAATGTTa                                                                                                                                                                                                                                                 <  2:182749/251‑1 (MQ=255)
       cccTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGAtt                                                                                                                                                                                                                                          <  2:271047/251‑1 (MQ=255)
       cccTGGCATGGAGGCGCTGGCGCAAGGGCAATTTTAGCAAGCGACTACTGGCAGCTCGGTGCTGTTAATCGAAAGCGTTGACGGCAGCGCTTTCAAAGATGTGTTATTCGCGCAACTTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACCTTTAACCCAGCGGCGCGAAGGCTCCCCGGACGGCACTATCCACTTGGGAACGCGCTTCGAAGGCACTGCGCTGTTACGTGAtt                                                                                                                                                                                                                                          >  2:361930/1‑251 (MQ=255)
             cATGGCGGAGCAGGCGCAAGGGCCATTCCAGCCAGCGAATAAAGGCAGCTCGGCGCTGCTCATCGAAAGCGTTGCCGGCAGCGATTTCAAAACTGGGTTCCTCGCCCAAATTCGACCAAAAGGTAATGAACGTCCTTCTGTGGTGGTGGGCGATTCCGGAAATTTAACCCAGATCCCCGACGGCCCCCAGGACGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGGTTCGTGCTTTcccc                                                                                                                                                                                                                                     >  2:342821/1‑248 (MQ=255)
                ggcggcGCTGGCGCCATGGCAATTCCAGCAAGCGAATAATGGCAGATCGGTGCTGTTCATCGAAAGAGTTGCCGGCAGAGCTTTCAAAGCTGTGTTCCTCGCGCAACTTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCCCCCAGGTCGTCCCTCTAAACCAGGGAACGCGCTTCGAAGGCATTGCACTGTTACGTGCTTTCCGCAtt                                                                                                                                                                                                                                  >  2:56957/1‑250 (MQ=255)
                ggcggcGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCCCTGCAATGTTACGTGATTTCCGCATTa                                                                                                                                                                                                                                 <  2:29863/251‑1 (MQ=255)
                 gcggcgCTGGCGCAAGGGCAAGTCCAGCAAGCGACTAATGGCAGTTCGGTGCTGTTCACCGAAAGCGTTGACGGCAGCGATTTCACAGATGTGTTCCGTGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCAACGGCTCCCAGGTCGTCACTATCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTAc                                                                                                                                                                                                                                <  2:66941/251‑1 (MQ=255)
                     cgcTGGCGCACGGGCAAGTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTACTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTATGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGAc                                                                                                                                                                                                                            >  2:57468/1‑251 (MQ=255)
                     cgcTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGAc                                                                                                                                                                                                                            <  1:235926/251‑1 (MQ=255)
                           cgcAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTcc                                                                                                                                                                                                                        >  2:260062/1‑249 (MQ=255)
                           cgcAAGGGCAATTCCAGCAAGAGACTACTGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTcc                                                                                                                                                                                                                        >  2:201731/1‑249 (MQ=255)
                            gcAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGCCTTCCAgg                                                                                                                                                                                                                     >  2:310885/1‑251 (MQ=255)
                                  gCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATc                                                                                                                                                                                                               <  2:311837/251‑1 (MQ=255)
                                  gCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATc                                                                                                                                                                                                               <  1:353715/251‑1 (MQ=255)
                                         cAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCTCTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGAt                                                                                                                                                                                                        <  2:385376/251‑1 (MQ=255)
                                                 gACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTc                                                                                                                                                                                                >  1:90895/1‑251 (MQ=255)
                                                 gACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTc                                                                                                                                                                                                >  2:231552/1‑251 (MQ=255)
                                                 gACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTc                                                                                                                                                                                                >  1:352001/1‑251 (MQ=255)
                                                    tAATGGCAGCTCGGGGCTGTTCATCGAAAGCGTTGACGGCAGCGCTTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAACGGAAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCCCACTCAAACAGGGAACGCGCTTCGAAGGCACTGCCATGTTACGTGCTTTCCGCATTACGGACTTCCCGGATTATCAGGCGCTCATTGGTCAcc                                                                                                                                                                                             >  1:392088/1‑251 (MQ=255)
                                                         gCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCg                                                                                                                                                                                        >  1:8235/1‑251 (MQ=255)
                                                               cGGTGCTGCTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGCTGTGTTCCTCGCGCAAATTCGACCAAAAGGGAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGAAATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGgcg                                                                                                                                                                                  >  1:277964/1‑251 (MQ=255)
                                                                 gTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATAAGGCGATCATTGGTCACCAGGCGGTGGCGCt                                                                                                                                                                                >  1:345521/1‑251 (MQ=255)
                                                                        tCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCCAATTCGACCAAAAGGTAATGCACGTCCTTATGTGGCGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGAcccc                                                                                                                                                                         >  1:185784/1‑250 (MQ=255)
                                                                            cGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACg                                                                                                                                                                     >  1:72451/1‑251 (MQ=255)
                                                                                     tGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGAc                                                                                                                                                             <  1:176615/250‑1 (MQ=255)
                                                                                        cGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGAGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGACTCCCAGGTCGTCAATCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGACCGATACCGAcc                                                                                                                                                            >  1:308275/1‑248 (MQ=255)
                                                                                          gCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATg                                                                                                                                                       >  1:133787/1‑251 (MQ=255)
                                                                                                aTTTCAAAGATGTGTTCCTCGCGCAAATTCGACAAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCTTCACTCTCAACCAGGGAACGCGCTTCGTAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTAGCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGGCCCGTACGATACCGACCAGATGGACATg                                                                                                                                                 <  2:348451/251‑1 (MQ=255)
                                                                                                 tttCAAAGCTGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCATTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATgc                                                                                                                                                >  1:153307/1‑251 (MQ=255)
                                                                                                       aGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACAt                                                                                                                                          >  1:122100/1‑251 (MQ=255)
                                                                                                         aTGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACTTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGGCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATtg                                                                                                                                        <  1:10917/251‑1 (MQ=255)
                                                                                                           gtgtTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATtgtg                                                                                                                                      <  1:310885/251‑1 (MQ=255)
                                                                                                             gtTCCTCGCGCAAATGCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGAAATTTCACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCGTTACGGTCTTCCATGATTATCAGGCGAACATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATAGACATGCGCACAGTGTGGa                                                                                                                                    <  2:77290/251‑1 (MQ=255)
                                                                                                              ttCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGCGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCCCTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGAGGCGCTCGACCCGAACGATACCGACCATATGGACATGCGCACATTGTGGa                                                                                                                                    >  1:248470/1‑250 (MQ=255)
                                                                                                              ttCCTCGAGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCTGTGGCGCTCGACCCGAACGCTACCGACCAGATGGACATGCGCCCATTGTGGa                                                                                                                                    >  2:251347/1‑250 (MQ=255)
                                                                                                                ccTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGTACTGCCAGGATTATCAGGCGATCATTGTTCACCAGGCGGTGGCGCTCTACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAaca                                                                                                                                 <  1:174841/251‑1 (MQ=255)
                                                                                                                    gcgcAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGa                                                                                                                             <  1:253313/251‑1 (MQ=255)
                                                                                                                          atgcGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGGTGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCGCTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGATCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGa                                                                                                                       <  2:199385/248‑1 (MQ=255)
                                                                                                                             cGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCg                                                                                                                    <  2:96900/251‑1 (MQ=255)
                                                                                                                              gACCAAAAGGTAATGCACGTCCTTCTGGGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGt                                                                                                                   <  2:217980/251‑1 (MQ=255)
                                                                                                                               aCCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTg                                                                                                                  <  1:377933/251‑1 (MQ=255)
                                                                                                                                    aaGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGCTCGTGCTcgc                                                                                                             >  1:393312/1‑251 (MQ=255)
                                                                                                                                        tAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCgcgccg                                                                                                         >  2:110113/1‑249 (MQ=255)
                                                                                                                                         aaTGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCgcgc                                                                                                           >  1:43375/1‑248 (MQ=255)
                                                                                                                                            gCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAgaac                                                                                                      <  1:231552/250‑1 (MQ=255)
                                                                                                                                             cACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAgaactg                                                                                                    >  1:172191/1‑251 (MQ=255)
                                                                                                                                              ccGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAgaactga                                                                                                   <  1:281973/250‑1 (MQ=255)
                                                                                                                                                      cTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCCGAACTGAACTggtgt                                                                                           >  2:112847/1‑248 (MQ=255)
                                                                                                                                                          ggttgtggCCGATTCCTGACATTTAACCCAGCTGCGCGACGGCTCCCAGGGCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGGGTATCa                                                                                       <  2:269794/247‑1 (MQ=255)
                                                                                                                                                           gtggtggCCGATTCCGTACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACTGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGGGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCAc                                                                                      <  1:361930/251‑1 (MQ=255)
                                                                                                                                                                     aTTCCTGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTc                                                                            <  2:209402/251‑1 (MQ=255)
                                                                                                                                                                     aTTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCATGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGAt                                                                                                                                                                   <  2:358326/164‑1 (MQ=255)
                                                                                                                                                                     aTTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGAt                                                                                                                                                                   >  1:358326/1‑164 (MQ=255)
                                                                                                                                                                      ttCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGCGTATCACGTTGGTATTcc                                                                           >  1:389976/1‑250 (MQ=255)
                                                                                                                                                                           gACATTTACCCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGTAGAACGGAACTGGTGTATCACGTGGGTATTCACCgtg                                                                      <  2:394311/251‑1 (MQ=255)
                                                                                                                                                                           gACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCgtg                                                                      <  2:133787/251‑1 (MQ=255)
                                                                                                                                                                           gACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCgtg                                                                      >  1:97666/1‑251 (MQ=255)
                                                                                                                                                                            aCATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGACCTGAACTGGTGTATCACGGTGGTATTCACCgtgt                                                                     >  2:253551/1‑251 (MQ=255)
                                                                                                                                                                                ttAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTat                                                                 >  1:46359/1‑251 (MQ=255)
                                                                                                                                                                                      cAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGc                                                           <  2:149063/251‑1 (MQ=255)
                                                                                                                                                                                            cgcgACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGTCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGATCTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTAt                                                     <  2:317860/251‑1 (MQ=255)
                                                                                                                                                                                            cgcgACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGCCACCGATCGTGATCGCGCAGCACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTAt                                                     >  1:362214/1‑251 (MQ=255)
                                                                                                                                                                                                 cGGCTCCCAGGGCGTCACTCTCAACCAGGGAAGGCGCTTCGAAGGCACTGCCATGTTACGAGATTTACGCATTACGGCCTTCCAGGATTATCAGGCGATCATTGGTCACAAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCAAGTTGGTATTAACCGCGTTTATGATGGCACTTATGGTCg                                                >  1:129538/1‑251 (MQ=255)
                                                                                                                                                                                                       ccAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTGCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTAGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGc                                          >  2:328879/1‑251 (MQ=255)
                                                                                                                                                                                                           gtcgtcACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAg                                      <  1:347064/251‑1 (MQ=255)
                                                                                                                                                                                                             tgtcACTCTCAACCAGGGAACGCGCTTGGAAGCCACTGCAATGTTACGTGATATCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGCACGATACCGACCAGCTGGACCTGCGCACATTGTGGAACACTGACACCGGTCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCAATTATGGTAGTACCgacgcgcg                                    >  2:20930/2‑244 (MQ=255)
                                                                                                                                                                                                              gtcACTCTCAGCCAGGGAACGCTCTTCCACGGCACTGCTATGTTAAGTGATTTCCGCATTACGCACTTCCAGGCTTATCAGGCGATCATTGGTCACCAGGCGACGGCGCTCGACCCGAAAGATACCGCCAAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTCTCACGATGGTATTCACCGTGTTTATGGTGGCACTTATGGTCGTACCGCTGAGCGt                                   >  2:161471/1‑251 (MQ=255)
                                                                                                                                                                                                                cACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTgg                                 <  1:372920/251‑1 (MQ=255)
                                                                                                                                                                                                                   tctcAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGtt                               >  1:398169/1‑250 (MQ=255)
                                                                                                                                                                                                                     tcAACCAGGGAACGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGCTTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGTTCGTACCGCTGAGCGTGGTTAAc                            >  2:178313/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaCGCGCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAgg                  >  2:121789/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cgcgCTTCGAAGGCACTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGa                <  2:394941/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cgcTTCGCAGGCACTGCAATGTTACGTGCTTTCCGCATTACGGACTTCCAGGATTCTCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGCCCCGAACGATACCGCCCAGCTGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGAcg              >  2:146963/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gAAGGCACTGCAATGTTACGTGATTTCCGCAATACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGAGGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGTGTACt        <  2:182062/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCTGCAATGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGTGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGTCGAt  >  2:313259/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
AAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTGGCCGATTCCGGACATTTAACCCAGCTGCGCGACGGCTCCCAGGTCGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCGACCAGATGGACATGCGCACATTGTGGAACACTGACACCGATCGTGCTCGCGCAGAACTGAACTGGCGTATCACGTTGGTATTCACCGTGTTTATGATGGCACTTATGGTCGTACCGCTGAGCGTGGTTAACCCACGTCAGGGACGCGTACTGTCGAT  >  NC_000913/4486624‑4487113

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