Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,452 C→A 9.3% G1279V (GGT→GTT)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,648,4520CA 9.3% 179.1 / 6.1 61G1279V (GGT→GTT) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base A (3/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTA  >  NC_000913/2648156‑2648749
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                          
atatCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAACCCg                                                                                                                                                                                                                                                                                                        >  2:556795/1‑301 (MQ=255)
 tatCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATCGCGCAGCAGACGCTCATTACCCCgg                                                                                                                                                                                                                                                                                                       >  2:505000/1‑301 (MQ=255)
 tatCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGAGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAACAGCGCAGCAAACGCTCATTACCCCgg                                                                                                                                                                                                                                                                                                       >  1:188129/1‑301 (MQ=255)
  atCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCACAGAAGACGCTCATCACCCCGGt                                                                                                                                                                                                                                                                                                      >  1:35936/1‑301 (MQ=255)
         catcaCTATTACGCGGCGTTTTCAGCGCCAGCGCACTCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTTACGGGCCAACACCAGCGCAGCGCAAGACTGTACGGCGAATTTACTGGCCTTGAGATTATCCGCGTCCGGGATCGCCATACTGCCCGGATCTTGACAATCACGCAGCAGACACTCCTTACCCCGCTTCATG‑Gc                                                                                                                                                                                                                                                                                              >  1:247542/1‑301 (MQ=255)
         catcaCTATTACACGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCGGCGGTAAACCAGAAGCGGCATCTGCACGATGCCCCCAGATTTCACGCAGTGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCGGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGATATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTACATAGCGCAGCAGACACTCATTACCCCGGTTAATg‑gg                                                                                                                                                                                                                                                                                              >  2:512324/1‑300 (MQ=255)
               tattaggggtggTTGTCTGCGCCTGCGCAATCGCTTCGTCACCACGCGTCGCATCACCCAGGGTTTGCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTGTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                                                                                                                                                                        <  1:112326/290‑2 (MQ=255)
                  tACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctgc                                                                                                                                                                                                                                                                                     <  2:624296/301‑5 (MQ=255)
                      cggcgTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACg                                                                                                                                                                                                                                                                                 <  2:35987/301‑1 (MQ=255)
                      cggcgTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACg                                                                                                                                                                                                                                                                                 <  2:574249/301‑1 (MQ=255)
                          tgtttCGGCGCCAGCGCAATCGCTTCTTCACCCCGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTGTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGt                                                                                                                                                                                                                                                                             <  1:238671/299‑1 (MQ=255)
                                    ccagcgcAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGACCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGttt                                                                                                                                                                                                                                                                   >  2:65067/1‑300 (MQ=255)
                                     ccgcgcGATCGCGTCTTCACCATGCTGCGCATCACCCAGTGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCg                                                                                                                                                                                                                                                                  <  1:64999/299‑1 (MQ=255)
                                               gCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGa                                                                                                                                                                                                                                                        >  2:229641/1‑301 (MQ=255)
                                                   tttgacctcGCGTCGCATCACCCTTGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                                                                                                                                                                    <  2:291804/297‑1 (MQ=255)
                                                            gcgTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                                                                                                                                                                           <  2:84473/301‑1 (MQ=255)
                                                                       cccATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCg                                                                                                                                                                                                                                >  2:469546/1‑301 (MQ=255)
                                                                                  cAGCGCAACGCCACGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCCCGATGCTCCCAGCTTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGACTTTGAGATTATCCGCGTACGGGATCGAAATCATGCCCGGAACTTGTAAATAGCGCAGCAGACGCCCATTACCCCGGTAAATG‑GCCTATGTCGGCCCCCTCTAACCCTCATCGCCTACCCCGACCAAGaacacaaacacgcacacccacaaccatcac                                                                                                                                                                                                                     >  1:33517/1‑273 (MQ=255)
                                                                                     cgcAACGCCACGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTCNCGCGCATTTACTGGCTTTCCGATTATCCCCGTACGGGATCGCNNNNATGCCCGGCTCNNGTACATAGCCCCGCANNNCCTCATTACNNNNNTCAATG‑GCATCTGTCCGCAAGCTGTAAACCCTTTCGCCTGCGCGAACCCGAACACCCAACACCTACACCATCACACAGTActct                                                                                                                                                                                                                  >  1:402200/1‑301 (MQ=255)
                                                                                            ccAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                                                                                                                                                                           >  2:197299/1‑301 (MQ=255)
                                                                                                  tGCAGCAGCGGTAAACCAGAATCGGCAGTTGCACGATGCTCCCGGATTTCACGCAGCGCACCCAGCGTAGCCTTTTGCTGACTGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAAGTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtttt                                                                                                                                                                                                     <  1:456754/301‑1 (MQ=255)
                                                                                                        agcGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                                                                                                                                                                               >  2:629021/1‑301 (MQ=255)
                                                                                                               aaCCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GAGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGCCCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGCCACCGTTTTTATCCCACAgcgc                                                                                                                                                                                        >  1:392244/1‑301 (MQ=255)
                                                                                                                          gCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGt                                                                                                                                                                             <  2:247809/301‑1 (MQ=255)
                                                                                                                                   ggATGCTCCCGGATTTCCCGCAGGGCACCCAGCGGTGCGTTTGGCGGACGGGCCAACACCGTCGTAGCGTAAGAGTGTACGGCGAATTGACTGGCTTTGAGATTATCCGCGTACGGGAACGACATGATGCCCGGATCTGGTAAATAGCGCAGCAGACGATCATGACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGtt                                                                                                                                                                    <  1:504328/300‑1 (MQ=255)
                                                                                                                                      tGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa                                                                                                                                                                 <  2:143573/301‑1 (MQ=255)
                                                                                                                                        cTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAtt                                                                                                                                                               >  2:632673/1‑301 (MQ=255)
                                                                                                                                          cccAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAACACGTTGCAtttt                                                                                                                                                             >  1:305054/1‑300 (MQ=255)
                                                                                                                                           ccgggttttACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGGGATTATCGGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCACCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                                                                                                                                                            <  1:100070/296‑1 (MQ=255)
                                                                                                                                           ccAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                                                                                                                                                            >  2:132533/1‑301 (MQ=255)
                                                                                                                                                   cACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAACTCCATCCCGTAAGCCGTCAGCCCGTACTCTTCGTCCCCGATTTTATCCCACAGCGCAAAGCCGCCGTTATCCCGTTGCCTTTGCAGCAAAcc                                                                                                                                                    >  1:379948/1‑300 (MQ=255)
                                                                                                                                                               ccggtggggcgttttGCTGAGGGGCGAACGCCAGCGCAGGGTAAGACTGTACGGGGAATTTACTGGGTTGGAGATTAGCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGat                                                                                                                                        <  1:584238/290‑1 (MQ=255)
                                                                                                                                                                                     gggcaacctaagcggagggTGAGACTGTACGGGGAATTGTCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAAGAGCGCAGCAGACGCTCATTACGCCGGTTAATG‑GCGTCTGTCGGCAAGCTGTAACCCTGTGCGCCTGCGCGGACCAGGAAATCCATCAGGTAAGCCGTCAGCCAGTACTCGTCGTCACCGTTTTTAACCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTtc                                                                                                                  <  1:511652/291‑1 (MQ=255)
                                                                                                                                                                                     ggCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTtc                                                                                                                  <  2:305458/301‑1 (MQ=255)
                                                                                                                                                                                      gCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGATCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACATCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATAACCACCGATGCACGGCGattat                                                                                                                 >  2:216372/1‑299 (MQ=255)
                                                                                                                                                                                      ccccacacCAGTGAGGCGTAAGTCTGTACGGCGAATTTACTGGCTGTGAGATTATCCGCGTACGGGATCGACATCATGCGCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTtct                                                                                                                 <  1:197056/297‑1 (MQ=255)
                                                                                                                                                                                        cacctccggggcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCtca                                                                                                               <  1:631796/291‑1 (MQ=255)
                                                                                                                                                                                        cAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GAGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCtca                                                                                                               <  2:392758/301‑1 (MQ=255)
                                                                                                                                                                                           cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCtcatca                                                                                                            >  2:432451/1‑301 (MQ=255)
                                                                                                                                                                                            cccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCAc                                                                                                           <  2:288841/300‑1 (MQ=255)
                                                                                                                                                                                              cggggcggcgtaggcttttACGGCGAATTTACTGGCTTTGTGATTATCCGCGTAGGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTg                                                                                                         <  1:556081/283‑1 (MQ=255)
                                                                                                                                                                                                                   ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTCTCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCAACACTGACCCCTTTGAAGCCCAAcccc                                                                                    >  1:342667/1‑301 (MQ=255)
                                                                                                                                                                                                                        ttttACTGTCTTTGAGATTATCCGCGTACGGGATTGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAg                                                                               <  2:187646/300‑1 (MQ=255)
                                                                                                                                                                                                                             gttgCTTGGTGATTATCCGGGTAGGGGATCGACAGCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTggg                                                                          <  1:363632/298‑1 (MQ=255)
                                                                                                                                                                                                                             cTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAActgcaggcgttggg                                                                          <  2:188371/301‑14 (MQ=255)
                                                                                                                                                                                                                             cTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTggg                                                                          >  2:546430/1‑301 (MQ=255)
                                                                                                                                                                                                                               ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCt                                                                        >  2:421114/1‑301 (MQ=255)
                                                                                                                                                                                                                                       atttgtCCGCTTACGGGTTGGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAActgcaggcgttgggcatcaaaggc                                                                <  1:431887/296‑24 (MQ=255)
                                                                                                                                                                                                                                        tttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGta                                                               <  2:178802/300‑1 (MQ=255)
                                                                                                                                                                                                                                                  tACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAACCGGGAAATGCCGATATCGACCGCTGCCCGGCGTTTCTCATCCCTGTCGCCTTTGCTGCCCAACCCCCGCAGTTCGGCTGCATTCGtcaaacatcacg                                                     >  1:364720/1‑290 (MQ=255)
                                                                                                                                                                                                                                                  tACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACg                                                     >  2:458974/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         tCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCCCATCCCTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAg                                              >  1:18467/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          cGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAga                                             >  2:141911/1‑300 (MQ=255)
                                                                                                                                                                                                                                                          cGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGACTTTGATACCCAACGCCTGCAGTTGGGCTGCGATGGTATAAAATGAAGGAAACAgg                                             >  1:122792/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                  tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGc                                     <  2:380444/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gcggACGCTCATTACCCCGGTTAATG‑GCTTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACgg             <  1:545757/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gcagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACgg             >  2:627956/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            agACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGaa           >  2:243740/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATa          >  2:210760/6‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                     ttACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGTCTGAGCGGACCATGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTa  <  2:402714/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     ttAACCGGGGTAATGAGCGTCTGTCGGCCCGTTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCCCGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCtttt   <  1:243473/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                          
ATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTA  >  NC_000913/2648156‑2648749

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