Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,595,418 C→A 39.7% S81* (TCG→TAG)  hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,595,4180CA39.7% 313.3 / 256.6 238S81* (TCG→TAG) hyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (77/66);  new base A (43/51);  total (120/118)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

TATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGA  >  CP009273/2595269‑2595567
                                                                                                                                                     |                                                                                                                                                     
tATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCtctc                                                                                                                                                       <  2:2456077/150‑1 (MQ=255)
   tAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTTATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGATTTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGc                                                                                                                                                    <  1:4169514/150‑1 (MQ=255)
       cTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCAttt                                                                                                                                                <  1:137221/150‑1 (MQ=255)
       cTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCAttt                                                                                                                                                >  1:444400/1‑150 (MQ=255)
        tttCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTGATGAATTTGCCGAGTTAACCCTGCGAATGGATATCCTCTCGGCATTTa                                                                                                                                               >  1:241065/1‑150 (MQ=255)
            cGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGt                                                                                                                                           <  2:3411140/150‑1 (MQ=255)
               cATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCt                                                                                                                                        <  2:2161050/150‑1 (MQ=255)
                  cacaTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGt                                                                                                                                     >  1:3884146/1‑150 (MQ=255)
                   acaTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGtt                                                                                                                                    <  1:2218890/150‑1 (MQ=255)
                    caTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa                                                                                                                                   >  2:2822970/1‑150 (MQ=255)
                     aTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTAt                                                                                                                                  >  1:1135728/1‑150 (MQ=255)
                         cTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTcc                                                                                                                              <  1:7669/150‑1 (MQ=255)
                          tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCt                                                                                                                             >  1:3862398/1‑150 (MQ=255)
                           ggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCtt                                                                                                                            <  1:831634/150‑1 (MQ=255)
                             tggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGGTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGc                                                                                                                          <  1:2445057/150‑1 (MQ=255)
                              gggCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCt                                                                                                                         >  1:377055/1‑150 (MQ=255)
                                gcggcgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCtgg                                                                                                                       <  1:1777961/150‑1 (MQ=255)
                                     gTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtggtg                                                                                                                  >  1:3989908/1‑150 (MQ=255)
                                       gATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                                                                                                                >  1:414117/1‑150 (MQ=255)
                                           ggCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGt                                                                                                            >  2:3427502/1‑150 (MQ=255)
                                            gCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGtt                                                                                                           <  1:284391/150‑1 (MQ=255)
                                            gCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGtt                                                                                                           <  1:1332707/150‑1 (MQ=255)
                                              atcatcAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtct                                                                                                         >  2:2276563/1‑150 (MQ=255)
                                              atcatcAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGGGGTTTGTtct                                                                                                         >  1:481445/1‑150 (MQ=255)
                                               tcatcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtctc                                                                                                        >  1:3822768/1‑150 (MQ=255)
                                                catcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtctct                                                                                                       >  2:2989759/1‑150 (MQ=255)
                                                  tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTctctct                                                                                                     >  1:979522/1‑150 (MQ=255)
                                                  tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTGTGTTctctct                                                                                                     >  1:1477521/1‑150 (MQ=255)
                                                     gCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTAtt                                                                                                  <  2:3339744/150‑1 (MQ=255)
                                                      cGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTc                                                                                                 <  2:2919753/150‑1 (MQ=255)
                                                      cGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTc                                                                                                 >  2:4022291/1‑150 (MQ=255)
                                                      cGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATGTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTc                                                                                                 >  1:4060907/1‑150 (MQ=255)
                                                       ggAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa                                                                                                <  2:3950352/150‑1 (MQ=255)
                                                        gAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt                                                                                               <  2:1470200/150‑1 (MQ=255)
                                                           ttACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGa                                                                                            <  2:2963642/150‑1 (MQ=255)
                                                           ttACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGa                                                                                            <  2:4216076/150‑1 (MQ=255)
                                                           ttACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGa                                                                                            <  2:2963628/150‑1 (MQ=255)
                                                            tACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAAACCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGaa                                                                                           >  2:1466277/1‑149 (MQ=255)
                                                              cGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACtt                                                                                         <  2:2415473/150‑1 (MQ=255)
                                                              cGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACtt                                                                                         <  2:1983259/150‑1 (MQ=255)
                                                               gCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTta                                                                                        <  2:1035006/150‑1 (MQ=255)
                                                                cAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTtat                                                                                       >  2:1319845/1‑150 (MQ=255)
                                                                 aaTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTTGGCATTTATGGTGCTGGTTATCTCCTTTCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTtata                                                                                      <  2:888844/150‑1 (MQ=255)
                                                                 aaTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTtata                                                                                      >  2:3693869/1‑150 (MQ=255)
                                                                 aaTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTtata                                                                                      <  2:601442/150‑1 (MQ=255)
                                                                  aTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTatat                                                                                     >  2:3772075/1‑150 (MQ=255)
                                                                   ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATg                                                                                    <  2:1632051/150‑1 (MQ=255)
                                                                    tACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATgc                                                                                   <  2:543105/150‑1 (MQ=255)
                                                                     aCATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATgcg                                                                                  >  1:652528/1‑150 (MQ=255)
                                                                      cATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATgcgc                                                                                 >  1:2729368/1‑150 (MQ=255)
                                                                           tGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATa                                                                            <  2:1583986/150‑1 (MQ=255)
                                                                              tgtAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGa                                                                         >  2:2100990/1‑150 (MQ=255)
                                                                               gtAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAg                                                                        <  1:1551943/150‑1 (MQ=255)
                                                                                 aaCTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAggg                                                                      >  1:148335/1‑150 (MQ=255)
                                                                                   cTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCa                                                                    <  2:1954808/150‑1 (MQ=255)
                                                                                          cGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTGTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGgcgc                                                             >  1:2837302/1‑150 (MQ=255)
                                                                                           gCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGcgcg                                                            >  1:3607027/1‑150 (MQ=255)
                                                                                              cGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCgcggcg                                                         <  2:3476133/150‑1 (MQ=255)
                                                                                                   tGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGAt                                                    >  2:3912278/1‑150 (MQ=255)
                                                                                                    gCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGTGAATACGAGGGCAAAGGCGCGGCGGCGATg                                                   <  1:4291084/150‑1 (MQ=255)
                                                                                                      ccccccTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATggg                                                 <  1:1267620/150‑1 (MQ=255)
                                                                                                             tttGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCttt                                          <  1:336771/150‑1 (MQ=255)
                                                                                                             tttGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCttt                                          >  1:4123624/1‑150 (MQ=255)
                                                                                                              ttGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTa                                         <  2:1078866/150‑1 (MQ=255)
                                                                                                               tGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTAt                                        <  1:3809061/150‑1 (MQ=255)
                                                                                                                 aaTTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGa                                      <  2:2629860/150‑1 (MQ=255)
                                                                                                                    ttGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAata                                   <  1:460269/150‑1 (MQ=255)
                                                                                                                    ttGCCGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAata                                   <  1:2778588/150‑1 (MQ=255)
                                                                                                                      gCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATAtt                                 <  2:3675722/150‑1 (MQ=255)
                                                                                                                       ccGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATAttt                                <  2:2336950/150‑1 (MQ=255)
                                                                                                                       ccGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATAttt                                >  2:3035131/1‑150 (MQ=255)
                                                                                                                        cGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATAtttt                               >  2:2165061/1‑150 (MQ=255)
                                                                                                                        cGATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATAtttt                               >  2:3084258/1‑150 (MQ=255)
                                                                                                                         gATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTc                              >  2:1040262/1‑150 (MQ=255)
                                                                                                                         gATTTAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTc                              <  2:4478763/150‑1 (MQ=255)
                                                                                                                          aTTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTca                             <  1:1587420/150‑1 (MQ=255)
                                                                                                                            ttAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatc                           <  2:3121656/150‑1 (MQ=255)
                                                                                                                             tAACCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatcg                          <  1:3465894/150‑1 (MQ=255)
                                                                                                                              aaCCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatcgc                         >  2:1385394/1‑150 (MQ=255)
                                                                                                                              aaCCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatcgc                         >  1:1042688/1‑150 (MQ=255)
                                                                                                                              aaCCCTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatcgc                         >  1:3759379/1‑150 (MQ=255)
                                                                                                                                 ccTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatcgcatc                      <  1:533416/150‑1 (MQ=255)
                                                                                                                                  cTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatcgcatcg                     >  1:2527675/1‑150 (MQ=255)
                                                                                                                                  cTGCGAATGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTcatcgcatcg                     >  1:4261417/1‑150 (MQ=255)
                                                                                                                                    gCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGAt                   <  2:560805/150‑1 (MQ=255)
                                                                                                                                     cGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATg                  >  2:688160/1‑150 (MQ=255)
                                                                                                                                      gAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATgg                 >  1:1659025/1‑150 (MQ=255)
                                                                                                                                      gAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGATTCTTTATGAATATTTTCATCGCATCGATgg                 >  2:4055779/1‑150 (MQ=255)
                                                                                                                                       aaTGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGt                >  2:3451926/1‑150 (MQ=255)
                                                                                                                                       aaTGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGt                <  2:2527675/150‑1 (MQ=255)
                                                                                                                                       aaTGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGAGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGt                <  2:966524/150‑1 (MQ=255)
                                                                                                                                        aTGGATAGCCTCTAGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGtt               >  2:4511905/1‑150 (MQ=255)
                                                                                                                                           gATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGcc            >  2:3685124/1‑150 (MQ=255)
                                                                                                                                           gATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGcc            >  2:2608354/1‑150 (MQ=255)
                                                                                                                                            aTAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGccc           >  2:4254812/1‑150 (MQ=255)
                                                                                                                                              aGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCctg         <  1:2656598/150‑1 (MQ=255)
                                                                                                                                              aGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCctg         >  2:38/1‑150 (MQ=255)
                                                                                                                                                     cGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGa  <  1:3137371/150‑1 (MQ=255)
                                                                                                                                                     cGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGa  <  1:619727/150‑1 (MQ=255)
                                                                                                                                                     cGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGa  >  2:3293184/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGA  >  CP009273/2595269‑2595567

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: