Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,641,817 T→G 8.6% N1277T (AAC→ACC)  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
*CP0092732,641,8170TG 8.6% 152.9 / 5.7 72N1277T (AAC→ACC) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base G (3/5);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

ACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGA  >  CP009273/2641672‑2641965
                                                                                                                                                 |                                                                                                                                                     
aCCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGAACGACATCCTGCCCGGTGCTTGTACATAGCGCAGCAGACGCCCATTACCCCGGGTaa                                                                                                                                                    >  1:161125/1‑149 (MQ=255)
 cccAGCGTAGCCTTTTGCTGACAGGCCAACCCCCGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAAt                                                                                                                                                   <  1:883791/149‑1 (MQ=255)
 cccAGCGGAGCCTTTTGCTGGCGGGGCAGCGCCAGCGCAGCGTAAGACTGTGCGGCGAATTTACTGGCTTTGAGATTATCTGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAAt                                                                                                                                                   <  2:507531/149‑1 (MQ=255)
 cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAAt                                                                                                                                                   <  2:826198/149‑1 (MQ=255)
 cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAAt                                                                                                                                                   <  1:603447/149‑1 (MQ=255)
 cccACCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAAt                                                                                                                                                   <  1:39224/149‑1 (MQ=255)
  ccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGCCATCATGCCCGGATCTTGTAAATAGCGCAGCAGCCGCTCATTACCCCGGTTAATg                                                                                                                                                  >  2:156620/1‑149 (MQ=255)
  ccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGGTAATg                                                                                                                                                  >  1:501197/1‑149 (MQ=255)
  ccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCACAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                  >  1:829485/1‑149 (MQ=255)
   cAGCGGAGCCTTTGTGTAACTGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg‑g                                                                                                                                                <  2:303064/149‑1 (MQ=255)
    aGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAAAAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑Gc                                                                                                                                               >  1:256976/1‑149 (MQ=255)
     gCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTCATG‑GCg                                                                                                                                              >  1:1021119/1‑149 (MQ=255)
      cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCCTACGGGATCGACATCATGCCCGGATCTTGTAAATGTCGCAGCAGACGCCCATTAACCCGGTCAACG‑CCGt                                                                                                                                             >  2:264810/1‑149 (MQ=255)
         acccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTg                                                                                                                                          <  2:818385/147‑1 (MQ=255)
         aGCCTTTTGCCGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATGTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTg                                                                                                                                          <  1:416042/149‑1 (MQ=255)
            cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGAAAGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctg                                                                                                                                       <  2:6038/149‑4 (MQ=38)
             ttttGCTGACGTGCCAACACCAGCGCATCGTATGACTGTACGGCGTGTTTTCGACCGTAGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctgc                                                                                                                                      <  1:346862/149‑5 (MQ=19)
             ttttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGCACGTGTAAATAGCGCACCCGCCGCTCATGACCCCGGTTCAGG‑GCGTCTGTccg                                                                                                                                      >  1:656339/1‑147 (MQ=255)
             ttttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                                                                      >  1:694502/1‑149 (MQ=255)
             ttttGCTGACGGGCCAACACCAGCGCAGCGGAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctgc                                                                                                                                      <  1:514876/149‑5 (MQ=38)
                  cTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGCGACCCACCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGc                                                                                                                                 <  1:904143/149‑1 (MQ=255)
                        ggCCAACACGAGCGCAGCGTAAGACTGGCCGTGGAATTTACTGGCTTTGAGATTATCCGCCTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAAc                                                                                                                           <  2:508168/149‑1 (MQ=255)
                             acacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACGTCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGTCACGCTGTAACCCGGt                                                                                                                      >  2:831520/1‑149 (MQ=255)
                                ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGGTAATG‑GGGTCTGTCGGCACGCTGTAACCCTGTTCg                                                                                                                   >  1:477138/1‑149 (MQ=255)
                                 cagcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                                  <  2:1033011/149‑1 (MQ=255)
                                  agcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGcc                                                                                                                 >  1:676931/1‑149 (MQ=255)
                                   gcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCt                                                                                                                >  2:634957/1‑149 (MQ=255)
                                    agcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCCGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACCTTGTAACCCTGTTCGCCTg                                                                                                               >  1:368038/2‑149 (MQ=255)
                                           tAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGAc                                                                                                        <  1:357621/149‑1 (MQ=255)
                                            aaGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGAcc                                                                                                       <  1:971303/149‑1 (MQ=255)
                                                   tACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGGCCAGGAAAt                                                                                                >  2:1114742/1‑149 (MQ=255)
                                                    aCGGCGAATTGACGGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGGTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATc                                                                                               <  2:161125/149‑1 (MQ=255)
                                                      ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCACCAGACGCGCATAACCCCGGTTAATG‑ACGTCTGTCGGCACGCAGTAACCCTGTTCGCCTGCGCGACCCAGGaaaccac                                                                                             >  2:314857/1‑145 (MQ=255)
                                                         gAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCa                                                                                          <  1:58699/149‑1 (MQ=255)
                                                         gAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCa                                                                                          <  1:271618/149‑1 (MQ=255)
                                                         gAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACACGGAAATCCATCa                                                                                          >  2:452837/1‑149 (MQ=255)
                                                                 cggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGcc                                                                                  <  1:421322/148‑1 (MQ=255)
                                                                   gCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGt                                                                                >  1:194722/1‑149 (MQ=255)
                                                                    ctgtGAGAGTACCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAAGAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGGCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                                                                               <  1:84982/146‑1 (MQ=255)
                                                                    cTTTGGGATTATCCGCGTACGGGATCGCCATCATGCCAGGATTTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                                                                               <  1:230908/149‑1 (MQ=255)
                                                                    cTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                                                                               <  1:634957/149‑1 (MQ=255)
                                                                    cTTTGAGAGTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                                                                               <  1:674367/149‑1 (MQ=255)
                                                                     tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGTGCATCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAACTCCATCACGTAAGCCGTCa                                                                              <  1:901782/149‑1 (MQ=255)
                                                                     tttGAGATTATCCGCGTACGGGATCGACAGCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCTGGTTAATG‑GCGTCCGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCa                                                                              <  1:70968/149‑1 (MQ=255)
                                                                            ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTa                                                                       <  2:51795/149‑1 (MQ=255)
                                                                              aTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAct                                                                     <  1:985607/149‑1 (MQ=255)
                                                                              aTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAct                                                                     <  1:609483/149‑1 (MQ=255)
                                                                                    gTACTGGATCGCCCGCATGCTCGGATCTTGGAAATAGCGACGCAGACGCTCATCACCCCGGTTAATG‑GCGGCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGt                                                               <  1:469662/149‑1 (MQ=255)
                                                                                       cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                            <  2:694502/149‑1 (MQ=255)
                                                                                        gggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCGTCACGTAAGCAGTCAGCCAGTACTCTTCGTCagc                                                           >  1:4873/1‑147 (MQ=255)
                                                                                          gATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCAAGTAAGCCGTCAACCAGTACTCTTCGTCACCGt                                                         >  2:1101288/1‑149 (MQ=255)
                                                                                           aTCGACATCGTGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTGAATG‑GCGTCTGTCGGCACGCTGTAACGCGGGTCGCCTGCGCGGACCAGGAAAACCATCACGTAAGCCGTCACGCAGTACTCTTCGTCACGGtt                                                        >  2:312161/1‑149 (MQ=255)
                                                                                           aTCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtt                                                        >  2:608736/1‑149 (MQ=255)
                                                                                           aTCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGCCACCGtt                                                        >  2:333937/1‑149 (MQ=255)
                                                                                            tCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGttt                                                       <  1:547235/149‑1 (MQ=255)
                                                                                               aCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTACTG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa                                                    <  2:1167482/149‑1 (MQ=255)
                                                                                                    atGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCGTCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTCTTATCCca                                               >  1:587318/1‑149 (MQ=255)
                                                                                                     tGCCCGGATCTTGTAAATAGCGCAGCAGAAGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTCTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcac                                              >  1:760003/1‑149 (MQ=255)
                                                                                                         cGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCCGTACTCTTCGTCCCCGTTTTTATCCCACAgcg                                          >  2:1141659/1‑149 (MQ=255)
                                                                                                          ggATCTTGTAAATAGCGCAGCAGACTCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTCAGCCGTCAGGCAGTACTCTTCGTCACCGTTTTTATCCCACAgcgc                                         >  1:693294/1‑149 (MQ=255)
                                                                                                             tctgcgaatttgcgcAGCAGACGCTCATTCCCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCGTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCaaa                                      <  1:410419/138‑1 (MQ=255)
                                                                                                             ccTTGTAAACAGCGCAGCAGTCGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCGGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCaaa                                      <  1:1005412/148‑1 (MQ=255)
                                                                                                               ttGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCATGTTCGCCTGCGCGGACCAGGAAATCGATCACGTAAGCCGCCAGCCAGTACTCTACGGCACCGGGTTGATCCCACAGCGCACAgc                                    >  2:107965/1‑149 (MQ=255)
                                                                                                                tGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAACTCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGCTTTTATCCCACAGCGCAAAgcc                                   >  2:827062/1‑149 (MQ=255)
                                                                                                                 gTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGAGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgccg                                  >  1:170194/1‑149 (MQ=255)
                                                                                                                  tAAATAGCGCAGCAGATGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCCGTAAACCTGTTCGACTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAGTCTTCGTCACCGGTTTTATCACACACCGCACAgcggc                                 >  2:873836/1‑149 (MQ=255)
                                                                                                                         cgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGACTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCGGACAGTACTCTTCGTCACCGTTTTTATCCCACAACGCAAGACCGCCGGTAtg                          >  1:505244/1‑148 (MQ=255)
                                                                                                                          gcggcagCCGCTCATTACCCCGGTTCATG‑CGGTCTGGACGCACGCTGTAACCCTGTTCCCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCa                         <  1:1067403/146‑1 (MQ=255)
                                                                                                                                gACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                   >  2:760912/1‑149 (MQ=255)
                                                                                                                                        gtcccccGGTTACCG‑GCGTCCGTCGGCAGCCTGTACCCCTGTTCGTCTGCGCCGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcag           <  1:174767/146‑1 (MQ=255)
                                                                                                                                              cGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGAACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACg     >  2:1046201/1‑149 (MQ=255)
                                                                                                                                                gTTGATG‑GCGTCTGTCGGCACGCTGTAGCCCTGTGCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACggg   <  1:360993/149‑1 (MQ=255)
                                                                                                                                                 ttAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGa  >  2:549209/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGA  >  CP009273/2641672‑2641965

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