Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,601,404 C→T 21.8% 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,4040CT21.8% 566.2 / 144.0 244L377F (CTC→TTC) hyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (86/104);  new base T (25/28);  total (112/132)
Fisher's exact test for biased strand distribution p-value = 8.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

GGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTA  >  CP009273/2601256‑2601553
                                                                                                                                                    |                                                                                                                                                     
ggAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCt                                                                                                                                                      <  2:3264688/150‑1 (MQ=255)
  aGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCa                                                                                                                                                    <  1:877473/150‑1 (MQ=255)
     tAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCaaaa                                                                                                                                                 >  1:4529526/1‑150 (MQ=255)
     tAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCaaaa                                                                                                                                                 >  1:2158790/1‑150 (MQ=255)
      aGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAAt                                                                                                                                                <  1:2790049/150‑1 (MQ=255)
      aGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAAt                                                                                                                                                >  2:1130495/1‑150 (MQ=255)
      aGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAAt                                                                                                                                                >  2:2756247/1‑150 (MQ=255)
          aTTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATg                                                                                                                                            >  1:637145/1‑150 (MQ=255)
          aTTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATg                                                                                                                                            <  1:4584303/150‑1 (MQ=255)
           ttGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGc                                                                                                                                           >  2:3415566/1‑150 (MQ=255)
             ggCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCa                                                                                                                                         <  1:4194554/150‑1 (MQ=255)
               cGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCAtt                                                                                                                                       <  2:2205693/150‑1 (MQ=255)
                gggCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCAttt                                                                                                                                      >  1:3738997/1‑150 (MQ=255)
                gggCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCAttt                                                                                                                                      <  1:3390408/150‑1 (MQ=255)
                gggCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCAttt                                                                                                                                      >  1:2557780/1‑150 (MQ=255)
                 ggCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGCTGCTCAAAATCATGCCATTTa                                                                                                                                     >  1:3603284/1‑150 (MQ=255)
                 ggCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTa                                                                                                                                     >  1:2969734/1‑150 (MQ=255)
                    cGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAccg                                                                                                                                  <  1:1710428/150‑1 (MQ=255)
                       tGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTAccgccg                                                                                                                               >  2:4369271/1‑150 (MQ=255)
                       tGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAccgccg                                                                                                                               >  2:3125940/1‑150 (MQ=255)
                       tGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAccgccg                                                                                                                               >  2:1120178/1‑150 (MQ=255)
                       tGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAccgccg                                                                                                                               <  1:3969674/150‑1 (MQ=255)
                        gggAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGt                                                                                                                              <  1:1130197/150‑1 (MQ=255)
                          gAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGc                                                                                                                            <  1:2087108/150‑1 (MQ=255)
                             tttttGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGt                                                                                                                         >  1:2645903/1‑150 (MQ=255)
                               tttGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGttt                                                                                                                       >  1:1961666/1‑150 (MQ=255)
                                ttGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTg                                                                                                                      <  1:4145921/150‑1 (MQ=255)
                                 tGCCGCTCTGCTGCACATCTTTCACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCTAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATGTACCGCCGTGCTGTTTgg                                                                                                                     <  2:1397680/150‑1 (MQ=255)
                                  gCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTggc                                                                                                                    >  1:3505082/1‑150 (MQ=255)
                                    cgcgCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTggcgg                                                                                                                  >  2:3626449/1‑150 (MQ=255)
                                      cgcTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGgcggcg                                                                                                                >  2:3603730/1‑150 (MQ=255)
                                          gctgcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgc                                                                                                            >  2:4110905/1‑150 (MQ=255)
                                            tgcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgcgc                                                                                                          >  1:2854860/1‑150 (MQ=255)
                                             gcACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgcgct                                                                                                         <  1:733699/150‑1 (MQ=255)
                                                   cTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTgcgctggcgct                                                                                                   >  2:3429175/1‑150 (MQ=255)
                                                     tAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGcgctggcgctgg                                                                                                 >  1:809198/1‑150 (MQ=255)
                                                     tAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGAACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCAGTGGGATGTTCAAAATCATGCCATTTACCGCCGAGCTGTTAGGCGGCGGAGcgctggcgctgg                                                                                                 >  1:392871/1‑150 (MQ=255)
                                                        aCCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAg                                                                                              <  1:2228089/150‑1 (MQ=255)
                                                        aCCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAg                                                                                              <  2:1941796/150‑1 (MQ=255)
                                                         ccACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAgg                                                                                             <  1:1844471/150‑1 (MQ=255)
                                                         ccACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAgg                                                                                             <  1:2189901/150‑1 (MQ=255)
                                                           acaGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGa                                                                                           >  1:2897708/1‑150 (MQ=255)
                                                           acaGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGa                                                                                           <  2:3973104/150‑1 (MQ=255)
                                                            caGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGAt                                                                                          <  1:1303234/150‑1 (MQ=255)
                                                              gTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgc                                                                                        >  1:1795033/1‑150 (MQ=255)
                                                               tCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACTTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgcc                                                                                       <  2:1763353/150‑1 (MQ=255)
                                                               tCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgcc                                                                                       <  1:4387196/150‑1 (MQ=255)
                                                               tCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgcc                                                                                       <  1:4387016/150‑1 (MQ=255)
                                                                  ggCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATgccgcc                                                                                    >  1:3309955/1‑150 (MQ=255)
                                                                   gCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGTGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGccc                                                                                   >  1:3537590/1‑150 (MQ=255)
                                                                   gCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGccc                                                                                   >  2:1322120/1‑150 (MQ=255)
                                                                    cAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCt                                                                                  >  1:616944/1‑150 (MQ=255)
                                                                     aaaaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCtt                                                                                 >  1:1262309/1‑150 (MQ=255)
                                                                     aaaaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCtt                                                                                 <  1:2304890/150‑1 (MQ=255)
                                                                     aaaaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCtt                                                                                 <  1:3826137/150‑1 (MQ=255)
                                                                       aaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCa                                                                               <  2:1043343/150‑1 (MQ=255)
                                                                       aaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCa                                                                               <  2:2346568/150‑1 (MQ=255)
                                                                       aaaCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCa                                                                               <  1:2675192/150‑1 (MQ=255)
                                                                         aCGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAAc                                                                             <  2:2013134/150‑1 (MQ=255)
                                                                           gctgctGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAt                                                                           <  2:3736317/150‑1 (MQ=255)
                                                                           gctgctGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAt                                                                           >  2:1437389/1‑150 (MQ=255)
                                                                            ctgctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAtt                                                                          <  2:3109997/150‑1 (MQ=255)
                                                                             tgctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAttt                                                                         <  2:2001497/150‑1 (MQ=255)
                                                                             tgctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAttt                                                                         >  1:1615155/1‑150 (MQ=255)
                                                                              gctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAtttt                                                                        <  2:3841232/150‑1 (MQ=255)
                                                                              gctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAtttt                                                                        >  1:2010115/1‑150 (MQ=255)
                                                                               ctgTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAttttt                                                                       >  1:1272383/1‑150 (MQ=255)
                                                                                  ttCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCtt                                                                    <  2:1860281/150‑1 (MQ=255)
                                                                                    cTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAg                                                                  >  2:3477028/1‑150 (MQ=255)
                                                                                      gCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCg                                                                <  1:1854017/150‑1 (MQ=255)
                                                                                        ggTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGac                                                              >  1:4112567/1‑149 (MQ=255)
                                                                                        ggTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGaa                                                              >  1:2990247/1‑150 (MQ=255)
                                                                                          ttCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAAtt                                                            <  2:1049728/150‑1 (MQ=255)
                                                                                            ccGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTa                                                          <  1:701875/150‑1 (MQ=255)
                                                                                              ggCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATg                                                        <  1:2503668/150‑1 (MQ=255)
                                                                                              ggCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATg                                                        >  2:2954025/1‑150 (MQ=255)
                                                                                              ggCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATg                                                        >  1:3001142/1‑150 (MQ=255)
                                                                                               gCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGa                                                       <  2:3933276/150‑1 (MQ=255)
                                                                                                 aaTGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGAcc                                                     >  2:2869668/1‑150 (MQ=255)
                                                                                                   tGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCAt                                                   <  1:120196/150‑1 (MQ=255)
                                                                                                    gTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCAtt                                                  <  1:2927392/150‑1 (MQ=255)
                                                                                                          cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAccgcc                                            >  2:4046595/1‑150 (MQ=255)
                                                                                                          cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTAccgcc                                            >  2:4208335/1‑150 (MQ=255)
                                                                                                               gTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACt                                       <  2:1712974/150‑1 (MQ=255)
                                                                                                                 aCGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTgg                                     <  2:2705115/150‑1 (MQ=255)
                                                                                                                  cGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGc                                    <  2:1499240/150‑1 (MQ=255)
                                                                                                                   ggCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCa                                   <  2:4516815/150‑1 (MQ=255)
                                                                                                                   ggCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCa                                   <  1:1322556/150‑1 (MQ=255)
                                                                                                                          gcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATc                            <  1:867091/150‑1 (MQ=255)
                                                                                                                               ctcAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACctg                       <  2:2513264/150‑1 (MQ=255)
                                                                                                                               ctcAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACctg                       >  2:3980338/1‑150 (MQ=255)
                                                                                                                                     gtcgtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGACTTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGAtt                 <  2:24277/150‑1 (MQ=255)
                                                                                                                                     gtcgtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGAtt                 >  1:1187547/1‑150 (MQ=255)
                                                                                                                                     gtcgtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGAtt                 >  1:2725649/1‑150 (MQ=255)
                                                                                                                                        gtcTGTGGGATGTTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATc              >  2:6819/1‑150 (MQ=255)
                                                                                                                                          cTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTGCTTAGCGAAGTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGt            >  2:2098665/1‑150 (MQ=255)
                                                                                                                                           tgtgGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTc           <  1:3795871/150‑1 (MQ=255)
                                                                                                                                              gggATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACAATTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCctg        <  1:220322/150‑1 (MQ=255)
                                                                                                                                                 aTGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCctgctg     <  2:60961/150‑1 (MQ=255)
                                                                                                                                                   gCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGtt   <  2:2149225/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTA  >  CP009273/2601256‑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: