Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,836,192 C→T 10.9% 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,1920CT10.9% 1151.3 / 91.9 358G71S (GGC→AGC) hycIprotease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base C (187/132);  new base T (23/16);  total (210/148)
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.89e-01

Only 100 of 360 total aligned reads displayed.

AAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTCCGCCAGCA  >  CP009273/2836043‑2836341
                                                                                                                                                     |                                                                                                                                                     
 aaGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCg                                                                                                                                                      >  1:74226/1‑150 (MQ=255)
   gCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAATCATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCggg                                                                                                                                                    >  2:2304358/1‑150 (MQ=255)
      gacgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGtt                                                                                                                                                 <  1:3080370/150‑1 (MQ=255)
       acgaTATCCGGCTGACTGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGttt                                                                                                                                                <  1:4634556/150‑1 (MQ=255)
       acgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGttt                                                                                                                                                <  1:2936670/150‑1 (MQ=255)
       acgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGttt                                                                                                                                                <  1:1545093/150‑1 (MQ=255)
        cgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTa                                                                                                                                               >  2:1051360/1‑150 (MQ=255)
         gaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAg                                                                                                                                              >  2:523667/1‑150 (MQ=255)
         gaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAg                                                                                                                                              <  2:3648892/150‑1 (MQ=255)
           tatCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGcc                                                                                                                                            >  2:210102/1‑150 (MQ=255)
                ggCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCata                                                                                                                                       >  1:2116477/1‑150 (MQ=255)
                ggCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCata                                                                                                                                       <  2:1633851/150‑1 (MQ=255)
                ggCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCata                                                                                                                                       >  1:1371540/1‑150 (MQ=255)
                 gCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCatat                                                                                                                                      >  2:3124491/1‑150 (MQ=255)
                    gAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCg                                                                                                                                   >  2:1253965/1‑150 (MQ=255)
                       tGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTgg                                                                                                                                >  2:2847309/1‑150 (MQ=255)
                        gCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGc                                                                                                                               <  1:4551370/150‑1 (MQ=255)
                        gCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGc                                                                                                                               >  2:2960816/1‑150 (MQ=255)
                              gAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGAc                                                                                                                         <  2:1551839/150‑1 (MQ=255)
                                 aaTCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAAt                                                                                                                      <  2:500450/150‑1 (MQ=255)
                                 aaTCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAAt                                                                                                                      >  1:1847206/1‑150 (MQ=255)
                                  aTCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATc                                                                                                                     <  1:3962522/150‑1 (MQ=255)
                                     aCTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATcagc                                                                                                                  >  1:1361373/1‑150 (MQ=255)
                                      cTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATcagca                                                                                                                 >  2:1726120/1‑150 (MQ=255)
                                        tCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGt                                                                                                               <  2:2058265/150‑1 (MQ=255)
                                        tCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGt                                                                                                               <  2:2174109/150‑1 (MQ=255)
                                           ccAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCgt                                                                                                            >  1:1955229/1‑150 (MQ=255)
                                              atatCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTc                                                                                                         <  2:4530671/150‑1 (MQ=255)
                                              atatCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTc                                                                                                         >  1:1881186/1‑150 (MQ=255)
                                              atatCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTc                                                                                                         >  1:4303666/1‑150 (MQ=255)
                                              atatCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTc                                                                                                         >  1:3886483/1‑150 (MQ=255)
                                              atatCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTc                                                                                                         >  1:38434/1‑150 (MQ=255)
                                               tatCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCg                                                                                                        >  1:3222725/1‑150 (MQ=255)
                                                atCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCgg                                                                                                       >  2:3046252/1‑150 (MQ=255)
                                                atCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCgg                                                                                                       >  2:2460845/1‑150 (MQ=255)
                                                   ttcttTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGgcg                                                                                                    >  1:2222948/1‑150 (MQ=255)
                                                    tcttTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGgcgc                                                                                                   >  2:3748170/1‑150 (MQ=255)
                                                     cttTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCa                                                                                                  >  2:3571053/1‑150 (MQ=255)
                                                        tCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGtt                                                                                               <  1:88918/150‑1 (MQ=255)
                                                         cAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTc                                                                                              <  1:3728014/150‑1 (MQ=255)
                                                         cAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTc                                                                                              >  1:2558237/1‑150 (MQ=255)
                                                          aaCTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCa                                                                                             >  2:2731914/1‑150 (MQ=255)
                                                          aaCTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCa                                                                                             <  2:4160567/150‑1 (MQ=255)
                                                             tGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCCTCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATAACCGTGGCGTCGACACTCAGCAGTCGTGTCGGGCGCAGTTCACgg                                                                                          >  2:467030/1‑150 (MQ=255)
                                                                 cGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAg                                                                                      >  2:1186549/1‑150 (MQ=255)
                                                                  gATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGc                                                                                     >  1:2873263/1‑150 (MQ=255)
                                                                     aaGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcgac                                                                                  >  1:84833/1‑150 (MQ=255)
                                                                       ggTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGcgacga                                                                                >  1:3751279/1‑150 (MQ=255)
                                                                       ggTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTAACGGATAGcgacga                                                                                >  2:332794/1‑150 (MQ=255)
                                                                        gTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGAt                                                                               >  1:1265106/1‑150 (MQ=255)
                                                                        gTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGAt                                                                               <  1:4104394/150‑1 (MQ=255)
                                                                         tAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATg                                                                              >  2:1504170/1‑150 (MQ=255)
                                                                          aaTTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGt                                                                             >  1:2509897/1‑150 (MQ=255)
                                                                            ttCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCg                                                                           <  2:3199362/150‑1 (MQ=255)
                                                                            ttCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCg                                                                           >  1:2528268/1‑150 (MQ=255)
                                                                             tCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGt                                                                          >  2:1932654/1‑150 (MQ=255)
                                                                              cAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGtt                                                                         >  2:4447382/1‑150 (MQ=255)
                                                                               aaCGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGttt                                                                        <  1:951227/150‑1 (MQ=255)
                                                                                aCGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGtttt                                                                       <  1:642948/150‑1 (MQ=255)
                                                                                aCGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGtttt                                                                       >  2:12990/1‑150 (MQ=255)
                                                                                    cATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCgg                                                                   >  2:882721/1‑150 (MQ=255)
                                                                                    cATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCgg                                                                   >  2:2230523/1‑150 (MQ=255)
                                                                                    cATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCgg                                                                   >  2:3717482/1‑150 (MQ=255)
                                                                                     aTGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGt                                                                  >  2:1894244/1‑150 (MQ=255)
                                                                                      tGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTg                                                                 >  2:3967409/1‑150 (MQ=255)
                                                                                        ttATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTgcg                                                               <  1:990860/150‑1 (MQ=255)
                                                                                          aTGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCt                                                             <  2:1392700/150‑1 (MQ=255)
                                                                                          aTGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCt                                                             >  2:3437745/1‑150 (MQ=255)
                                                                                          aTGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCt                                                             >  1:267756/1‑150 (MQ=255)
                                                                                           tGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTa                                                            <  1:2426091/150‑1 (MQ=255)
                                                                                             ggTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATTGCGACGATGTCGTTTTCCGGTGCGCTAcc                                                          >  1:1041897/1‑150 (MQ=255)
                                                                                             ggTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc                                                          <  2:267756/150‑1 (MQ=255)
                                                                                                aGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAccgcc                                                       >  1:1008540/1‑150 (MQ=255)
                                                                                                aGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAccgcc                                                       <  1:2344764/150‑1 (MQ=255)
                                                                                                 gTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAccgccg                                                      <  2:2574701/150‑1 (MQ=255)
                                                                                                  tcatcaTAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCTGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGt                                                     >  1:4037204/1‑150 (MQ=255)
                                                                                                      catAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAAt                                                 <  1:2273235/150‑1 (MQ=255)
                                                                                                       atAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATc                                                <  1:521726/150‑1 (MQ=255)
                                                                                                                tcGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAg                                       <  1:3125790/150‑1 (MQ=255)
                                                                                                                  ggCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGtt                                     >  1:3975157/1‑150 (MQ=255)
                                                                                                                       tatCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCtt                                <  2:2399229/150‑1 (MQ=255)
                                                                                                                       tatCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCtt                                >  1:2733770/1‑150 (MQ=255)
                                                                                                                        atcatcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCttt                               <  2:3886483/150‑1 (MQ=255)
                                                                                                                        atcatcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCttt                               <  1:115356/150‑1 (MQ=255)
                                                                                                                         tcatcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTc                              <  2:428991/150‑1 (MQ=255)
                                                                                                                          catcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcg                             <  1:3676532/150‑1 (MQ=255)
                                                                                                                          catcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcg                             >  1:4513725/1‑150 (MQ=255)
                                                                                                                          catcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcg                             <  1:3534926/150‑1 (MQ=255)
                                                                                                                            tcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggc                           <  2:3517256/150‑1 (MQ=255)
                                                                                                                            tcCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggc                           <  2:4191221/150‑1 (MQ=255)
                                                                                                                              cGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggcgc                         >  1:4413680/1‑150 (MQ=255)
                                                                                                                               gggTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggcgcg                        >  1:2660866/1‑150 (MQ=255)
                                                                                                                               gggTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggcgcg                        <  2:4107995/150‑1 (MQ=255)
                                                                                                                               gggTCGATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTcggcgcg                        <  1:3244162/150‑1 (MQ=255)
                                                                                                                                    gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCgcggcggc                   >  1:1199604/1‑150 (MQ=255)
                                                                                                                                    gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCgcggcggc                   >  2:706249/1‑150 (MQ=255)
                                                                                                                                     atgatgCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCgcggcggcg                  >  1:1973961/1‑150 (MQ=255)
                                                                                                                                         tgCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACt              <  1:3282673/150‑1 (MQ=255)
                                                                                                                                             gATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTc          >  1:1694213/1‑150 (MQ=255)
                                                                                                                                              aTCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACCGCCGTCAATCACCACCCAGTTACCTTTCGGCGCGGCGGCGCACTTTTcc         <  1:2412059/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: