Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,836,192 C→T 9.0% G71S (GGC→AGC)  hycI ← protease involved in processing C‑terminal end of HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,836,1920CT 9.0% 1009.5 / 59.3 300G71S (GGC→AGC) hycIprotease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base C (145/128);  new base T (18/9);  total (163/137)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 308 total aligned reads displayed.

AAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGCCAGCA  >  CP009273/2836043‑2836341
                                                                                                                                                     |                                                                                                                                                     
  aGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTgg                                                                                                                                                     >  2:2353698/1‑150 (MQ=255)
   gCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCggg                                                                                                                                                    <  2:3074689/150‑1 (MQ=255)
     cgacgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGt                                                                                                                                                  >  1:480382/1‑150 (MQ=255)
     cgacgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGt                                                                                                                                                  >  1:2082609/1‑150 (MQ=255)
     cgacgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGt                                                                                                                                                  >  1:2791471/1‑150 (MQ=255)
       acgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGttt                                                                                                                                                <  2:1494028/150‑1 (MQ=255)
           tatCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGcc                                                                                                                                            >  2:3380116/1‑150 (MQ=255)
           tatCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGcc                                                                                                                                            >  2:3353826/1‑150 (MQ=255)
           tatCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGcc                                                                                                                                            <  2:2921059/150‑1 (MQ=255)
            atCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGccc                                                                                                                                           >  1:864624/1‑150 (MQ=255)
            atCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGccc                                                                                                                                           >  2:4593038/1‑150 (MQ=255)
                     aaTGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGt                                                                                                                                  >  2:3129019/1‑150 (MQ=255)
                      aTGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCGACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTg                                                                                                                                 >  2:4143071/1‑150 (MQ=255)
                      aTGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCGACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTg                                                                                                                                 >  2:545622/1‑150 (MQ=255)
                      aTGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTg                                                                                                                                 <  1:2548709/150‑1 (MQ=255)
                      aTGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTg                                                                                                                                 <  1:2478520/150‑1 (MQ=255)
                        gCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGc                                                                                                                               <  1:109112/150‑1 (MQ=255)
                         ccGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCg                                                                                                                              >  2:1903208/1‑150 (MQ=255)
                          cGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGt                                                                                                                             >  1:2999508/1‑150 (MQ=255)
                            aGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCg                                                                                                                           >  2:3524275/1‑150 (MQ=255)
                             ggAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGa                                                                                                                          >  2:1673966/1‑150 (MQ=255)
                                 aaTCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAAt                                                                                                                      >  2:3248511/1‑150 (MQ=255)
                                 aaTCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAAt                                                                                                                      >  2:1587432/1‑150 (MQ=255)
                                  aTCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATc                                                                                                                     >  1:4187388/1‑150 (MQ=255)
                                    cACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATcag                                                                                                                   >  2:2480499/1‑150 (MQ=255)
                                     aCTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATcagc                                                                                                                  >  2:2827429/1‑150 (MQ=255)
                                      cTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATcagca                                                                                                                 <  2:2055210/150‑1 (MQ=255)
                                       ttCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATcagcag                                                                                                                <  2:434771/150‑1 (MQ=255)
                                           ccAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGAAAATCAGCAGTCgt                                                                                                            <  2:1237271/150‑1 (MQ=255)
                                             aaTATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCgtgt                                                                                                          >  1:2418457/1‑150 (MQ=255)
                                                atCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCgg                                                                                                       >  2:2205422/1‑150 (MQ=255)
                                                 tcttctTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCggg                                                                                                      >  1:781047/1‑150 (MQ=255)
                                                     cttTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCa                                                                                                  <  1:1172764/150‑1 (MQ=255)
                                                       ttCAACTGGTCGATAAGGTAATTCATCGGCATGTTATGTGTAGTCATCATAAACATCTCGGCGATATCTTGCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCatt                                                                                                <  1:1819018/150‑3 (MQ=255)
                                                       ttCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGt                                                                                                <  1:2700899/150‑1 (MQ=255)
                                                       ttCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGt                                                                                                >  2:141450/1‑150 (MQ=255)
                                                        tCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGtt                                                                                               <  1:1640645/150‑1 (MQ=255)
                                                          aaCTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCa                                                                                             <  2:646942/150‑1 (MQ=255)
                                                            cTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACg                                                                                           >  2:2421686/1‑150 (MQ=255)
                                                            cTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACg                                                                                           >  2:3835969/1‑150 (MQ=255)
                                                              ggTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGa                                                                                         >  2:829081/1‑150 (MQ=255)
                                                                tCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATa                                                                                       <  2:4063543/150‑1 (MQ=255)
                                                                 cGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAg                                                                                      <  2:2132050/150‑1 (MQ=255)
                                                                 cGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAg                                                                                      <  2:194651/150‑1 (MQ=255)
                                                                  gATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGc                                                                                     >  2:4366340/1‑150 (MQ=255)
                                                                  gATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGc                                                                                     >  1:1474574/1‑150 (MQ=255)
                                                                   aTAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcg                                                                                    >  1:4227774/1‑150 (MQ=255)
                                                                    tAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcga                                                                                   >  1:2001127/1‑150 (MQ=255)
                                                                     aaGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATTATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcgac                                                                                  <  1:4453771/150‑1 (MQ=255)
                                                                     aaGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcgac                                                                                  <  1:1838859/150‑1 (MQ=255)
                                                                      aGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGAGATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcgacg                                                                                 <  1:1311914/150‑1 (MQ=255)
                                                                       ggTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcgacga                                                                                >  2:3916590/1‑150 (MQ=255)
                                                                        gTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGAt                                                                               >  1:826265/1‑150 (MQ=255)
                                                                         tAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATg                                                                              >  2:74146/1‑150 (MQ=255)
                                                                          aaTTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGt                                                                             >  1:4461831/1‑150 (MQ=255)
                                                                           aTTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTc                                                                            <  2:3104631/150‑1 (MQ=255)
                                                                            ttCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCg                                                                           >  1:1787396/1‑150 (MQ=255)
                                                                               aaCGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGttt                                                                        >  2:4475455/1‑150 (MQ=255)
                                                                                  ggCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTcc                                                                     <  1:535570/150‑1 (MQ=255)
                                                                                   gCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCg                                                                    <  2:2597414/150‑1 (MQ=255)
                                                                                   gCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCg                                                                    >  1:3272979/1‑150 (MQ=255)
                                                                                     aTGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGt                                                                  >  1:2336682/1‑150 (MQ=255)
                                                                                      tGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTg                                                                 <  1:319778/150‑1 (MQ=255)
                                                                                       gTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTgc                                                                <  1:287939/150‑1 (MQ=255)
                                                                                        ttATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTgcg                                                               >  2:3736259/1‑150 (MQ=255)
                                                                                        ttATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTgcg                                                               >  2:2249463/1‑150 (MQ=255)
                                                                                        ttATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTgcg                                                               >  2:3991794/1‑150 (MQ=255)
                                                                                          aTGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCt                                                             >  2:363871/1‑150 (MQ=255)
                                                                                          aTGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCt                                                             >  1:1305135/1‑150 (MQ=255)
                                                                                           tGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTa                                                            >  1:4056427/1‑150 (MQ=255)
                                                                                              gTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAccg                                                         >  1:1379426/1‑150 (MQ=255)
                                                                                                 gTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAccgccg                                                      >  1:1584766/1‑150 (MQ=255)
                                                                                                   catcatAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTc                                                    <  1:4018339/150‑1 (MQ=255)
                                                                                                   catcatAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTc                                                    <  2:477183/150‑1 (MQ=255)
                                                                                                     tcatAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCaa                                                  <  1:2282262/150‑1 (MQ=255)
                                                                                                       atAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATc                                                >  1:2195648/1‑150 (MQ=255)
                                                                                                            cATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATcaccac                                           <  2:2562448/150‑1 (MQ=255)
                                                                                                             aTCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCaccacc                                          >  1:3675631/1‑150 (MQ=255)
                                                                                                                tcGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAg                                       >  2:1782407/1‑150 (MQ=255)
                                                                                                                   gCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTa                                    <  2:3029337/150‑1 (MQ=255)
                                                                                                                     gATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTAcc                                  >  2:2624289/1‑150 (MQ=255)
                                                                                                                     gATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTAcc                                  >  1:3427619/1‑150 (MQ=255)
                                                                                                                         tcatcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTc                              <  1:3454256/150‑1 (MQ=255)
                                                                                                                         tcatcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTc                              >  2:1044662/1‑150 (MQ=255)
                                                                                                                           atcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcgg                            <  1:4526404/150‑1 (MQ=255)
                                                                                                                           atcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcgg                            <  1:3796053/150‑1 (MQ=255)
                                                                                                                               gggTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggcgcg                        >  2:670185/1‑150 (MQ=255)
                                                                                                                               gggTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggcgcg                        >  2:2283537/1‑150 (MQ=255)
                                                                                                                                  tCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggcgcggcg                     <  1:3383051/150‑1 (MQ=255)
                                                                                                                                   cGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCgcggcgg                    >  2:1653339/1‑150 (MQ=255)
                                                                                                                                   cGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCgcggcgg                    >  1:3759411/1‑150 (MQ=255)
                                                                                                                                    gatgatGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCgcggcggc                   >  1:1725218/1‑150 (MQ=255)
                                                                                                                                     atgatgCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTAGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCgcggcggcg                  <  1:4112670/150‑1 (MQ=255)
                                                                                                                                        atgCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCAc               >  2:1686497/1‑150 (MQ=255)
                                                                                                                                                ctcGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGc       >  2:2702945/1‑150 (MQ=255)
                                                                                                                                                 tcGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGcc      <  1:3105698/150‑1 (MQ=255)
                                                                                                                                                 tcGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGcc      >  2:1720871/1‑150 (MQ=255)
                                                                                                                                                 tcGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGcc      <  1:1482372/150‑1 (MQ=255)
                                                                                                                                                    ccGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGCcagc   >  1:655701/1‑150 (MQ=255)
                                                                                                                                                     cGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGCcagca  <  1:3946460/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGCCAGCA  >  CP009273/2836043‑2836341

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