Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,601,404 C→T 10.7% L377F (CTC→TTC)  hyfF → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,601,4040CT10.7% 1091.1 / 90.5 328L377F (CTC→TTC) hyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (147/146);  new base T (15/20);  total (162/166)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

TGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTA  >  CP009273/2601255‑2601553
                                                                                                                                                     |                                                                                                                                                     
tGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGc                                                                                                                                                       <  2:1033000/150‑1 (MQ=255)
   aGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCa                                                                                                                                                    >  2:808376/1‑150 (MQ=255)
      tAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCaaaa                                                                                                                                                 >  2:2178600/1‑150 (MQ=255)
          cATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCAt                                                                                                                                             <  1:192322/150‑1 (MQ=255)
           aTTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATg                                                                                                                                            >  1:4830020/1‑150 (MQ=255)
            ttGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGc                                                                                                                                           <  2:4504119/150‑1 (MQ=255)
                cGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCAtt                                                                                                                                       <  1:3297682/150‑1 (MQ=255)
                 gggCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCAttt                                                                                                                                      >  2:2415007/1‑150 (MQ=255)
                  ggCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTa                                                                                                                                     <  1:3056685/150‑1 (MQ=255)
                   gCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTAc                                                                                                                                    <  2:3510477/150‑1 (MQ=255)
                   gCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAc                                                                                                                                    <  2:4672815/150‑1 (MQ=255)
                    ccGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAcc                                                                                                                                   <  1:1148337/150‑1 (MQ=255)
                      gCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAccgc                                                                                                                                 >  2:3648856/1‑150 (MQ=255)
                         gggAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGt                                                                                                                              >  1:4346270/1‑150 (MQ=255)
                         gggAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGt                                                                                                                              <  1:1318799/150‑1 (MQ=255)
                         gggAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGt                                                                                                                              <  1:1294438/150‑1 (MQ=255)
                          ggAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTg                                                                                                                             <  2:2766021/150‑1 (MQ=255)
                          ggAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTg                                                                                                                             <  1:4476115/150‑1 (MQ=255)
                             aTTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCACCGTCGTCTGTGGGATGCTCAACATCATGCCATTTACCGCCGTGCTg                                                                                                                          <  2:3143462/150‑1 (MQ=255)
                             aTTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTg                                                                                                                          >  1:3832681/1‑150 (MQ=255)
                                  tGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTgg                                                                                                                     >  2:2841185/1‑150 (MQ=255)
                                        gctgctGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGcggcgg                                                                                                               >  2:2135265/1‑150 (MQ=255)
                                         ctgctgCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGt                                                                                                              <  2:4656192/150‑1 (MQ=255)
                                         ctgctgCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGt                                                                                                              <  2:2466958/150‑1 (MQ=255)
                                          tgctgcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTg                                                                                                             <  2:1561438/150‑1 (MQ=255)
                                           gctgcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgc                                                                                                            <  2:3003625/150‑1 (MQ=255)
                                           gctgcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgc                                                                                                            <  1:1904048/150‑1 (MQ=255)
                                           gctgcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgc                                                                                                            <  2:4438344/150‑1 (MQ=255)
                                             tgcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGATTGGCGGCGGTgcgc                                                                                                          >  2:1275891/1‑150 (MQ=255)
                                               cacaTCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgcgctg                                                                                                        >  1:2371819/1‑150 (MQ=255)
                                               cacaTCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgcgctg                                                                                                        >  2:10951/1‑150 (MQ=255)
                                               cacaTCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgcgctg                                                                                                        <  1:3480114/150‑1 (MQ=255)
                                                  aTCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgcgctggcg                                                                                                     <  1:3517712/150‑1 (MQ=255)
                                                      tAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGcgctggcgctgg                                                                                                 >  1:3971487/1‑150 (MQ=255)
                                                       aaaCCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCgctggcgctggc                                                                                                >  2:1830291/1‑150 (MQ=255)
                                                       aaaCCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCgctggcgctggc                                                                                                <  2:3596619/150‑1 (MQ=255)
                                                          ccACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAgg                                                                                             >  2:3221644/1‑150 (MQ=255)
                                                             caGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGAt                                                                                          <  2:3832681/150‑1 (MQ=255)
                                                             caGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGAt                                                                                          <  2:543469/150‑1 (MQ=255)
                                                              aGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATg                                                                                         <  1:2267297/150‑1 (MQ=255)
                                                                 cTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgccg                                                                                      >  1:4596789/1‑150 (MQ=255)
                                                                 cTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgccg                                                                                      <  1:821129/150‑1 (MQ=255)
                                                                  tGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgccgc                                                                                     >  1:4202344/1‑150 (MQ=255)
                                                                    gCAAAAACGCTGCTGGTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGccc                                                                                   <  2:5126703/150‑1 (MQ=255)
                                                                     cAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCt                                                                                  >  2:344790/1‑150 (MQ=255)
                                                                      aaaaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCtt                                                                                 <  2:1216016/150‑1 (MQ=255)
                                                                       aaaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTc                                                                                >  1:4415224/1‑150 (MQ=255)
                                                                        aaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCGTCa                                                                               >  1:3227282/1‑150 (MQ=255)
                                                                         aaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCaa                                                                              >  1:5007222/1‑150 (MQ=255)
                                                                         aaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCaa                                                                              >  2:1354540/1‑150 (MQ=255)
                                                                             ctgctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAtt                                                                          <  1:4706045/150‑1 (MQ=255)
                                                                             ctgctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAtt                                                                          >  2:1936264/1‑150 (MQ=255)
                                                                              tgctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAttt                                                                         >  2:4771101/1‑150 (MQ=255)
                                                                               gctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAtttt                                                                        >  1:3928713/1‑150 (MQ=255)
                                                                                ctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAttttt                                                                       >  2:3360917/1‑150 (MQ=255)
                                                                                ctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAttttt                                                                       <  1:1277643/150‑1 (MQ=255)
                                                                                    tCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTa                                                                   <  2:3728367/150‑1 (MQ=255)
                                                                                      tGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGc                                                                 <  1:699585/150‑1 (MQ=255)
                                                                                       gCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCg                                                                <  1:1801584/150‑1 (MQ=255)
                                                                                          gTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAAt                                                             >  2:711680/1‑150 (MQ=255)
                                                                                          gTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAAt                                                             >  1:2967882/1‑150 (MQ=255)
                                                                                             ccGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTa                                                          <  1:4929792/150‑1 (MQ=255)
                                                                                               ggCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATg                                                        >  2:4280428/1‑150 (MQ=255)
                                                                                                 cAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGAc                                                      <  1:2007880/150‑1 (MQ=255)
                                                                                                 cAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGAc                                                      >  2:4979473/1‑150 (MQ=255)
                                                                                                  actGTACTGCTCAAGTACGGCACGCGCGATCTCACCGTCGTCTGGGGGATGTTCAAAATCATGCCATGTACCGCCGTGCTGTTTGGCGGCGGTGCGCGGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGAcc                                                     <  1:3607881/148‑1 (MQ=255)
                                                                                                     gTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCAtt                                                  <  2:1033605/150‑1 (MQ=255)
                                                                                                      tACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTa                                                 >  1:1021527/1‑150 (MQ=255)
                                                                                                         tGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAccg                                              <  2:3555598/150‑1 (MQ=255)
                                                                                                           cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAccgcc                                            >  2:2215929/1‑150 (MQ=255)
                                                                                                           cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAccgcc                                            >  1:4890945/1‑150 (MQ=255)
                                                                                                            tCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAccgccg                                           >  2:1719169/1‑150 (MQ=255)
                                                                                                            tCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAccgccg                                           <  2:3633248/150‑1 (MQ=255)
                                                                                                              aaGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGa                                         >  1:3432649/1‑150 (MQ=255)
                                                                                                               aGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGAc                                        >  1:5138666/1‑150 (MQ=255)
                                                                                                                  aCGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTgg                                     >  2:2962571/1‑150 (MQ=255)
                                                                                                                    ggCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCa                                   <  1:1643389/150‑1 (MQ=255)
                                                                                                                     gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCAc                                  >  1:3133494/1‑150 (MQ=255)
                                                                                                                     gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCAc                                  <  2:4127508/150‑1 (MQ=255)
                                                                                                                      cACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACg                                 <  1:1375836/150‑1 (MQ=255)
                                                                                                                       aCGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGt                                >  1:4983584/1‑150 (MQ=255)
                                                                                                                       aCGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGt                                >  1:4281253/1‑150 (MQ=255)
                                                                                                                        cgcgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTa                               <  2:3987067/150‑1 (MQ=255)
                                                                                                                          cgcgATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAAt                             <  1:2135265/150‑1 (MQ=255)
                                                                                                                          cgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAAt                             <  1:2470867/150‑1 (MQ=255)
                                                                                                                             gATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCAc                          <  1:2705544/150‑1 (MQ=255)
                                                                                                                             gATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCAc                          >  1:3342507/1‑150 (MQ=255)
                                                                                                                               tctcAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACct                        <  1:606637/150‑1 (MQ=255)
                                                                                                                               tctcAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACct                        <  1:3190850/150‑1 (MQ=255)
                                                                                                                                      gtcgtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGAtt                 >  2:5118550/1‑150 (MQ=255)
                                                                                                                                       tcgtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTa                <  1:1217888/150‑1 (MQ=255)
                                                                                                                                       tcgtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTa                >  2:1630548/1‑150 (MQ=255)
                                                                                                                                         gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATc              <  2:4622278/150‑1 (MQ=255)
                                                                                                                                         gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATc              >  2:2495878/1‑150 (MQ=255)
                                                                                                                                         gtcTGTGGGATGCTCAAAATCATGCCATGTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATc              <  2:2684381/150‑1 (MQ=255)
                                                                                                                                          tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCg             >  1:3801601/1‑150 (MQ=255)
                                                                                                                                              tgGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCct         <  1:1057559/150‑1 (MQ=255)
                                                                                                                                                 gATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCctgct      <  2:3080288/150‑1 (MQ=255)
                                                                                                                                                   tGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGt    <  1:473609/150‑1 (MQ=255)
                                                                                                                                                     cTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTa  >  2:268845/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTA  >  CP009273/2601255‑2601553

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: