Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,605,773 A→T 8.2% D172V (GAT→GTT)  hyfR → hydrogenase‑4 transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,605,7730AT 8.2% 1026.3 / 64.1 292D172V (GAT→GTT) hyfRhydrogenase‑4 transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (135/133);  new base T (8/16);  total (143/149)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 307 total aligned reads displayed.

CTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGT  >  CP009273/2605625‑2605921
                                                                                                                                                    |                                                                                                                                                    
 tggcGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGAt                                                                                                                                                     <  1:2044603/149‑1 (MQ=255)
 tggcGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGAt                                                                                                                                                     >  2:606511/1‑149 (MQ=255)
 tggcGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGAt                                                                                                                                                     >  1:4104330/1‑149 (MQ=255)
   gcgcGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCt                                                                                                                                                   >  1:5031146/1‑149 (MQ=255)
       ggATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATc                                                                                                                                               >  2:1535922/1‑149 (MQ=255)
         aTAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGc                                                                                                                                             >  2:3304994/1‑149 (MQ=255)
          tAATCACCCCGATTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCt                                                                                                                                            <  1:1743508/149‑1 (MQ=255)
            aTCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGa                                                                                                                                          >  1:338627/1‑149 (MQ=255)
            aTCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGa                                                                                                                                          <  1:4381044/149‑1 (MQ=255)
            aTCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGa                                                                                                                                          <  2:1338597/149‑1 (MQ=255)
            aTCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATATTGATCTCGACGATCTGATCGCTGa                                                                                                                                          >  1:526509/1‑149 (MQ=255)
                  ccGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGc                                                                                                                                    >  1:1246613/1‑149 (MQ=255)
                  ccGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGc                                                                                                                                    >  2:953830/1‑149 (MQ=255)
                  ccGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGc                                                                                                                                    <  2:2508905/149‑1 (MQ=255)
                     gTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCg                                                                                                                                 <  2:1760780/149‑1 (MQ=255)
                      ttCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGt                                                                                                                                >  1:975912/1‑149 (MQ=255)
                       tCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTg                                                                                                                               <  1:2075082/149‑1 (MQ=255)
                         tGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTgag                                                                                                                             >  2:883835/1‑149 (MQ=255)
                         tGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTgag                                                                                                                             <  2:4890521/149‑1 (MQ=255)
                             aTCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGTGAGATcc                                                                                                                         >  1:4476431/1‑149 (MQ=255)
                               cATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCaat                                                                                                                       >  1:4079500/1‑149 (MQ=255)
                                 tGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGTGAGATCcatca                                                                                                                     <  2:4582958/149‑1 (MQ=255)
                                 tGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCcatca                                                                                                                     <  1:2349406/149‑1 (MQ=255)
                                  ggCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCcatcat                                                                                                                    >  2:5064024/1‑149 (MQ=255)
                                     tCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCAtttt                                                                                                                 <  1:685264/149‑1 (MQ=255)
                                       aaTGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCAtttttt                                                                                                               <  2:1907047/149‑1 (MQ=255)
                                        atgatgAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTc                                                                                                              >  2:4776137/1‑149 (MQ=255)
                                         tgatgaAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCg                                                                                                             >  1:3522178/1‑149 (MQ=255)
                                           atgaAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGt                                                                                                           >  2:3473617/1‑149 (MQ=255)
                                             gaagaaCAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCt                                                                                                         <  1:2029865/149‑1 (MQ=255)
                                                gaaCAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGc                                                                                                      <  2:792352/149‑1 (MQ=255)
                                                 aaCAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCt                                                                                                     <  1:427238/149‑1 (MQ=255)
                                                  aCAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCtt                                                                                                    >  1:3285173/1‑149 (MQ=255)
                                                     gcagGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTtcag                                                                                                 <  1:606511/149‑1 (MQ=255)
                                                        ggTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTtcagtca                                                                                              >  2:2459608/1‑149 (MQ=255)
                                                          tGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTcagtcagt                                                                                            <  2:3156456/149‑1 (MQ=255)
                                                          tGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTcagtcagt                                                                                            >  1:1828634/1‑149 (MQ=255)
                                                           ggCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTa                                                                                           <  2:4916093/149‑1 (MQ=255)
                                                             cgcTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATg                                                                                         <  2:2102771/149‑1 (MQ=255)
                                                              gcTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATgg                                                                                        >  1:2962255/1‑149 (MQ=255)
                                                              gcTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATgg                                                                                        >  2:5243276/1‑149 (MQ=255)
                                                                 gCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTAc                                                                                     >  1:3699167/1‑149 (MQ=255)
                                                                   gTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTg                                                                                   <  1:1414948/149‑1 (MQ=255)
                                                                       gCTTGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCg                                                                               <  1:2824932/149‑1 (MQ=255)
                                                                           gagaAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGatca                                                                           <  1:1505312/149‑1 (MQ=255)
                                                                           gagaAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGatca                                                                           <  1:728344/149‑1 (MQ=255)
                                                                             gaAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGatcatc                                                                         <  1:2225746/149‑1 (MQ=255)
                                                                                aGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGaa                                                                      <  2:3818072/149‑1 (MQ=255)
                                                                                aGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGaa                                                                      >  1:4567052/1‑149 (MQ=255)
                                                                                aGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGaa                                                                      >  1:3461575/1‑149 (MQ=255)
                                                                                aGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGaa                                                                      >  2:4676127/1‑149 (MQ=255)
                                                                                aGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGaa                                                                      >  2:1065227/1‑149 (MQ=255)
                                                                                 gATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGaaa                                                                     <  1:3590880/149‑1 (MQ=255)
                                                                                 gATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGaaa                                                                     >  1:3976226/1‑149 (MQ=255)
                                                                                    cgcgATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGaa                                                                  <  2:1577147/149‑1 (MQ=255)
                                                                                    cgcgATCATCATCGGGGACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGaa                                                                  <  1:3269668/149‑1 (MQ=255)
                                                                                     gcgATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAAc                                                                 <  2:4484825/149‑1 (MQ=255)
                                                                                        atcatcatcGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACgag                                                              <  2:1449155/149‑1 (MQ=255)
                                                                                        atcatcatcGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACgag                                                              >  2:2348322/1‑149 (MQ=255)
                                                                                          catcatcGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGaa                                                            <  2:4959440/149‑1 (MQ=255)
                                                                                          catcatcGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGaa                                                            <  1:3882537/149‑1 (MQ=255)
                                                                                               tcGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCa                                                       <  1:4754119/149‑1 (MQ=255)
                                                                                               tcGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCa                                                       <  1:2854075/149‑1 (MQ=255)
                                                                                               tcGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCa                                                       >  1:1010470/1‑149 (MQ=255)
                                                                                                cGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAg                                                      <  1:2617158/149‑1 (MQ=255)
                                                                                                  ggTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGcc                                                    >  2:1948672/1‑149 (MQ=255)
                                                                                                  ggTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGACCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGcc                                                    <  2:75767/149‑1 (MQ=255)
                                                                                                     aCTGGTTGATATTACTACTGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCtgt                                                 <  1:4948314/149‑1 (MQ=255)
                                                                                                     aCTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGTTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCtgt                                                 <  1:3683031/149‑1 (MQ=255)
                                                                                                     aCTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCtgt                                                 >  2:1334405/1‑149 (MQ=255)
                                                                                                     aCTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCtgt                                                 <  1:3121234/149‑1 (MQ=255)
                                                                                                      cTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCtgtg                                                >  2:4345452/1‑149 (MQ=255)
                                                                                                      cTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATAGAAAGAACGAGAAGTTCAGCCtgtg                                                >  2:5108809/1‑149 (MQ=255)
                                                                                                        ggTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTgtggt                                              >  1:1644054/1‑149 (MQ=255)
                                                                                                        ggTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATAGAAAGAACGAGAAGTTCAGCCTgtggt                                              <  2:4104330/149‑1 (MQ=255)
                                                                                                          ttGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGc                                            <  1:2419610/149‑1 (MQ=255)
                                                                                                             atatTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAga                                         >  2:1856808/1‑148 (MQ=255)
                                                                                                             atatTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGc                                         >  1:258111/1‑149 (MQ=255)
                                                                                                             atatTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGc                                         >  1:3767455/1‑149 (MQ=255)
                                                                                                                ttACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGAGAt                                      <  2:3080940/149‑1 (MQ=255)
                                                                                                                 tactacTGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATc                                     <  2:2226072/149‑1 (MQ=255)
                                                                                                                  aCTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCt                                    <  1:1867083/149‑1 (MQ=255)
                                                                                                                   cTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCtt                                   <  1:4676127/149‑1 (MQ=255)
                                                                                                                   cTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATAGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCtt                                   >  2:4863392/1‑149 (MQ=255)
                                                                                                                      aTGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCt                                >  2:3601779/1‑149 (MQ=255)
                                                                                                                        gCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGc                              <  1:2115817/149‑1 (MQ=255)
                                                                                                                               tGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACAt                       >  1:4833596/1‑149 (MQ=255)
                                                                                                                                     aTCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTGTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCg                 <  2:4444465/149‑1 (MQ=255)
                                                                                                                                      tCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGACGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCgt                >  1:3005145/1‑149 (MQ=255)
                                                                                                                                      tCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCgt                >  1:1023974/1‑149 (MQ=255)
                                                                                                                                      tCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCgt                >  2:2066968/1‑149 (MQ=255)
                                                                                                                                       cTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCgtg               <  1:2242495/149‑1 (MQ=255)
                                                                                                                                       cTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCgtg               >  1:1937631/1‑149 (MQ=255)
                                                                                                                                         tGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTc             <  1:1085003/149‑1 (MQ=255)
                                                                                                                                         tGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATAGAAAGAACGAGAAATTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTc             <  2:2962255/149‑1 (MQ=255)
                                                                                                                                          gATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCt            >  2:967384/1‑149 (MQ=255)
                                                                                                                                                 acgaTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGc     <  2:3283699/149‑1 (MQ=255)
                                                                                                                                                 acgaTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGc     >  2:2680308/1‑149 (MQ=255)
                                                                                                                                                  cgaTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCt    <  2:4737725/149‑1 (MQ=255)
                                                                                                                                                    aTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGAGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGt  <  1:1659741/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGCGCGGATAATCACCCCGGTTCTGCAATCATGGCTCAATGATGAAGAACAGCAGGTGGCGCTGCGTTTGCTGGAGAAAGATCGCGATCATCATCGGGTACTGGTTGATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGT  >  CP009273/2605625‑2605921

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: