Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,597,173 T→G 100% L666R (CTG→CGG)  hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,597,1730TG97.6% 1120.5 / 4.2 286L666R (CTG→CGG) hyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (147/132);  total (150/136)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 294 total aligned reads displayed.

GGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGT  >  CP009273/2597026‑2597319
                                                                                                                                                   |                                                                                                                                                  
gggCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGtgcggc                                                                                                                                                  >  2:1129549/1‑150 (MQ=255)
    taagagCACAGAACCATTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCACCGGCTGTCCAAAGAAATACAGCATCTGCAAAGAGGCGTCTTTCGTCTATATAGCCTGTATGTGGTCGCAGCACTGGTTGTGCGGCGAAt                                                                                                                                              >  2:719767/3‑150 (MQ=255)
    tCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAAt                                                                                                                                              <  1:949512/150‑1 (MQ=255)
     cAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATc                                                                                                                                             <  1:1809586/150‑1 (MQ=255)
      agagCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGAAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCg                                                                                                                                            >  1:1040746/1‑150 (MQ=255)
        agCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCt                                                                                                                                          <  1:2912342/150‑1 (MQ=255)
        agCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCt                                                                                                                                          <  1:3501504/150‑1 (MQ=255)
        agCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCt                                                                                                                                          <  1:2319911/150‑1 (MQ=255)
          cacaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTAt                                                                                                                                        >  1:308462/1‑150 (MQ=255)
          cacaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTAt                                                                                                                                        >  1:2799279/1‑150 (MQ=255)
          cacaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTAt                                                                                                                                        >  2:1669355/1‑150 (MQ=255)
           acaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTAtt                                                                                                                                       <  2:1092278/150‑1 (MQ=255)
           acaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTAtt                                                                                                                                       <  2:4393933/150‑1 (MQ=255)
           acaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTAtt                                                                                                                                       >  2:4864534/1‑150 (MQ=255)
           acaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTAtt                                                                                                                                       <  1:5154144/150‑1 (MQ=255)
            caGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTg                                                                                                                                      >  2:4261583/1‑150 (MQ=255)
             aGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGc                                                                                                                                     >  1:1671961/1‑150 (MQ=255)
                aCCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGt                                                                                                                                  <  1:3410198/150‑1 (MQ=255)
                  ccTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCt                                                                                                                                <  1:8939/150‑1 (MQ=255)
                  ccTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCt                                                                                                                                >  1:3476082/1‑150 (MQ=255)
                        gggATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACCGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGaa                                                                                                                          >  1:5325018/1‑150 (MQ=255)
                         ggAAGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGaaa                                                                                                                         <  1:1314155/150‑1 (MQ=255)
                          gATGAGCGGGTGATCCGCCCGATCGTGAGCGCCACCCTACGGCTGGCCAAAGACATACCGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAAt                                                                                                                        <  2:4437474/150‑1 (MQ=255)
                           aTGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATc                                                                                                                       <  2:1429113/150‑1 (MQ=255)
                              aGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCAcc                                                                                                                    <  2:3016496/150‑1 (MQ=255)
                              aGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCAcc                                                                                                                    <  2:509362/150‑1 (MQ=255)
                                cGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCAt                                                                                                                  >  1:1094133/1‑150 (MQ=255)
                                cGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCAt                                                                                                                  <  1:2988953/150‑1 (MQ=255)
                                cGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCAt                                                                                                                  <  1:287978/150‑1 (MQ=255)
                                     gATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGAc                                                                                                             >  2:2471435/1‑150 (MQ=255)
                                      aTCCGCCCCATCGTGAGCGCCACCCATAGGCTGGCCAAAGAAATTAAGCAGCTTAAAAGCGGCGACTTTCGTCTCTTTTGCCTGTATGTGTTCGACTCTATGGTTGTGCGGCTTATCGCTATTGCCGTCTAAGGAAATCTCCATGAGACa                                                                                                            <  2:494352/150‑1 (MQ=255)
                                      aTCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATg                                                                                                                 >  2:4512230/1‑145 (MQ=255)
                                       tCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGTGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACaa                                                                                                           <  1:235129/150‑1 (MQ=255)
                                          gCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTTCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTTAGGAAATCACCATGAGACAAAct                                                                                                        <  1:4136437/150‑1 (MQ=255)
                                             ccATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCtt                                                                                                     <  1:1575017/150‑1 (MQ=255)
                                              cATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCttt                                                                                                    >  2:1963322/1‑150 (MQ=255)
                                              cATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCttt                                                                                                    <  1:3719494/150‑1 (MQ=255)
                                               aTCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTg                                                                                                   >  2:3859496/1‑150 (MQ=255)
                                                tCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGc                                                                                                  <  2:3968234/150‑1 (MQ=255)
                                                 cGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCg                                                                                                 >  2:2328891/1‑150 (MQ=255)
                                                   tGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGAc                                                                                               <  1:1739118/150‑1 (MQ=255)
                                                     aGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACgg                                                                                             >  2:3324222/1‑150 (MQ=255)
                                                     aGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACgg                                                                                             <  2:838394/150‑1 (MQ=255)
                                                        gcCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGata                                                                                          <  2:2751751/150‑1 (MQ=255)
                                                            cccATCGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTg                                                                                      <  1:656124/150‑1 (MQ=255)
                                                            cccAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTg                                                                                      <  2:3119369/150‑1 (MQ=255)
                                                             ccAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTgg                                                                                     >  2:3541440/1‑150 (MQ=255)
                                                             ccAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTgg                                                                                     <  1:4753786/150‑1 (MQ=255)
                                                              cAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGt                                                                                    >  2:1724456/1‑150 (MQ=255)
                                                                 cGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCGATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAt                                                                                 >  1:795075/1‑150 (MQ=255)
                                                                  ggCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAtt                                                                                <  1:1892230/150‑1 (MQ=255)
                                                                  ggCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAtt                                                                                <  1:2528489/150‑1 (MQ=255)
                                                                    cTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAtttt                                                                              <  1:5235476/150‑1 (MQ=255)
                                                                      ggCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTg                                                                            <  2:3331472/150‑1 (MQ=255)
                                                                         cAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGt                                                                         >  2:3770931/1‑150 (MQ=255)
                                                                          aaaGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGtt                                                                        >  2:4139785/1‑150 (MQ=255)
                                                                           aaGAAATACAGCATCTGCAAAGCGGCGTCTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTa                                                                       <  1:4571347/150‑1 (MQ=255)
                                                                            aGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAg                                                                      >  1:3271006/1‑150 (MQ=255)
                                                                            aGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAg                                                                      >  2:2056429/1‑150 (MQ=255)
                                                                            aGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAg                                                                      >  1:450315/1‑150 (MQ=255)
                                                                             gAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGc                                                                     <  1:3111761/150‑1 (MQ=255)
                                                                             gAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGAGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGc                                                                     <  1:4154612/150‑1 (MQ=255)
                                                                               aaTACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGcac                                                                   >  2:1257096/1‑150 (MQ=255)
                                                                               aaTACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGcac                                                                   >  1:4667866/1‑150 (MQ=255)
                                                                                   cAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGc                                                               >  1:2667314/1‑150 (MQ=255)
                                                                                    aGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGcc                                                              <  2:4793247/150‑1 (MQ=255)
                                                                                       aTCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTg                                                           >  2:571486/1‑150 (MQ=255)
                                                                                         cTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGAt                                                         <  2:551147/150‑1 (MQ=255)
                                                                                         cTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGGCATTTTTGCGTTAGCACAGGCCGTGAt                                                         >  2:2163120/1‑150 (MQ=255)
                                                                                          tGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGAtt                                                        >  1:2427509/1‑150 (MQ=255)
                                                                                          tGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGAtt                                                        >  1:1391112/1‑150 (MQ=255)
                                                                                            cAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTTATCGCTTTTGCCGTCTAAGGACATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTct                                                      <  2:4601411/150‑1 (MQ=255)
                                                                                             aaaGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTctg                                                     >  2:2355759/1‑150 (MQ=255)
                                                                                                gcggcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACTGGCCGTGATTctgctg                                                  >  2:4670995/1‑150 (MQ=255)
                                                                                                  ggcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGAt                                                >  2:4082475/1‑150 (MQ=255)
                                                                                                  ggcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGAt                                                >  2:933979/1‑150 (MQ=255)
                                                                                                  ggcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGAt                                                <  1:5178690/150‑1 (MQ=255)
                                                                                                     gACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCt                                             <  2:654056/150‑1 (MQ=255)
                                                                                                       cTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTaa                                           <  2:1789184/150‑1 (MQ=255)
                                                                                                        tttCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGAACAGGCCGTGATTCTGCTGATGCTAAc                                          >  1:4502775/1‑150 (MQ=255)
                                                                                                          tCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAAccc                                        <  1:3192686/150‑1 (MQ=255)
                                                                                                           cGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAAcccc                                       >  1:3739776/1‑150 (MQ=255)
                                                                                                           cGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAAcccc                                       >  2:2720874/1‑150 (MQ=255)
                                                                                                            gTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCa                                      >  2:4724192/1‑150 (MQ=255)
                                                                                                              ctctATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGAAGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTTATTCTGCTGATGCTAACCCCACt                                    >  2:5271790/1‑150 (MQ=255)
                                                                                                               tctATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACtt                                   >  1:4307357/1‑150 (MQ=255)
                                                                                                               tctATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACtt                                   >  1:395214/1‑150 (MQ=255)
                                                                                                                  aTTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACttttt                                <  1:2257902/150‑1 (MQ=255)
                                                                                                                  aTTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACttttt                                <  1:4903382/150‑1 (MQ=255)
                                                                                                                    tGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTAc                              >  1:4508257/1‑150 (MQ=255)
                                                                                                                        tGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGt                          <  1:5032960/150‑1 (MQ=255)
                                                                                                                          tATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTAt                        <  1:359161/150‑1 (MQ=255)
                                                                                                                          tATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTAt                        >  1:4621230/1‑150 (MQ=255)
                                                                                                                            tgtgGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTAttt                      >  2:2588431/1‑150 (MQ=255)
                                                                                                                               ggTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTccc                   >  1:5285272/1‑150 (MQ=255)
                                                                                                                                 tCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCgg                 <  1:4864534/150‑1 (MQ=255)
                                                                                                                                   gccgcACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCa               <  2:105681/150‑1 (MQ=255)
                                                                                                                                   gccgcACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCa               >  2:1299559/1‑150 (MQ=255)
                                                                                                                                          tGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATAcgc        >  2:357609/1‑150 (MQ=255)
                                                                                                                                                tgcggcTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGt  >  2:3538477/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGT  >  CP009273/2597026‑2597319

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