Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,605,882 C→A 5.5% S208R (AGC→AGA hyfR → hydrogenase‑4 transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,605,8820CA 5.5% 1147.2 / 30.5 298S208R (AGC→AGAhyfRhydrogenase‑4 transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (159/122);  new base A (9/8);  total (168/130)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGTTATTACTCTT  >  CP009273/2605733‑2606031
                                                                                                                                                     |                                                                                                                                                     
gATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGc                                                                                                                                                       <  1:3651656/150‑1 (MQ=255)
  tatTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGTTACTGGGCGATCATAGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGa                                                                                                                                                     <  2:4878909/150‑1 (MQ=255)
  tatTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGa                                                                                                                                                     >  1:1239788/1‑150 (MQ=255)
   atTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGAt                                                                                                                                                    <  2:3405158/150‑1 (MQ=255)
     tACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCt                                                                                                                                                  <  1:1777476/150‑1 (MQ=255)
      aCTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCtt                                                                                                                                                 >  2:4669139/1‑150 (MQ=255)
        tAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGAGATCTTTc                                                                                                                                               >  1:3792754/1‑150 (MQ=255)
         aaTGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCt                                                                                                                                              <  2:1071150/150‑1 (MQ=255)
         aaTGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCt                                                                                                                                              >  1:2618234/1‑150 (MQ=255)
           tGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGc                                                                                                                                            >  1:1102426/1‑150 (MQ=255)
            gCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGcc                                                                                                                                           <  2:72402/150‑1 (MQ=255)
             cAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCt                                                                                                                                          <  1:1915869/150‑1 (MQ=255)
               gTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTca                                                                                                                                        <  1:112739/150‑1 (MQ=255)
                tGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTcac                                                                                                                                       <  1:4531608/150‑1 (MQ=255)
                  cTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACAt                                                                                                                                     >  2:909053/1‑150 (MQ=255)
                  cTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACAt                                                                                                                                     >  2:780209/1‑150 (MQ=255)
                   tGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACAtt                                                                                                                                    >  1:2789191/1‑150 (MQ=255)
                    gTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATAGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATtg                                                                                                                                   >  1:4930810/1‑150 (MQ=255)
                      cacaTCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCTCATtgtg                                                                                                                                 <  1:3141846/150‑1 (MQ=255)
                      cacaTCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATtgtg                                                                                                                                 <  1:1012443/150‑1 (MQ=255)
                       acaTCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGAGATCTTTCTGCCTCACATTGTGc                                                                                                                                <  2:1623579/150‑1 (MQ=255)
                       aaaTCTTGATCTCGAAGATCTTATCGCTGAAGTCTCTCGTGATATCCATGATTTTTTAGGTCTGGCTTCAGTCAGTATGGTACTGTGCGATCATCGAAATAACTATAAGTTCAGCCTGTTGTGCAGCGATCTTTCTGCCTCAATTTGTGc                                                                                                                                <  2:3571801/148‑1 (MQ=255)
                         aTCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCgt                                                                                                                              <  1:3221974/150‑1 (MQ=255)
                                tctcGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGcc                                                                                                                       <  2:4405454/150‑1 (MQ=255)
                                    gacgaTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGc                                                                                                                   >  1:726893/1‑150 (MQ=255)
                                       gaTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGt                                                                                                                <  1:61411/150‑1 (MQ=255)
                                       gaTCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGt                                                                                                                <  1:927002/150‑1 (MQ=255)
                                          cTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATg                                                                                                             <  1:3872797/150‑1 (MQ=255)
                                              tCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTg                                                                                                         >  1:2896630/1‑150 (MQ=255)
                                                   gACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGaa                                                                                                    >  1:4312996/1‑150 (MQ=255)
                                                     cGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAg                                                                                                  <  1:1321849/150‑1 (MQ=255)
                                                      gTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAgt                                                                                                 <  1:652866/150‑1 (MQ=255)
                                                        cGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAgtgt                                                                                               >  2:1230688/1‑150 (MQ=255)
                                                         gCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTa                                                                                              <  2:4311460/150‑1 (MQ=255)
                                                             gTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGc                                                                                          >  1:1894834/1‑150 (MQ=255)
                                                             gTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGc                                                                                          <  2:2550268/150‑1 (MQ=255)
                                                               gagaTCCATCATTTTTTCGGTCTGGCTTGAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTg                                                                                        <  1:2372245/150‑1 (MQ=255)
                                                               gagaTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTg                                                                                        >  1:5003114/1‑150 (MQ=255)
                                                               gagaTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTg                                                                                        >  2:195989/1‑150 (MQ=255)
                                                               gagaTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTg                                                                                        >  2:1652918/1‑150 (MQ=255)
                                                               gagaTCCATCATGTTTTCGGTCGGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTg                                                                                        <  1:4965854/150‑1 (MQ=255)
                                                                agaTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGa                                                                                       >  1:3427330/1‑150 (MQ=255)
                                                                 gaTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGac                                                                                      >  2:792112/1‑150 (MQ=255)
                                                                 gaTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGac                                                                                      <  1:8158/150‑1 (MQ=255)
                                                                 gaTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGac                                                                                      <  1:4267254/150‑1 (MQ=255)
                                                                  aTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGaca                                                                                     >  1:1425939/1‑150 (MQ=255)
                                                                  aTCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGaca                                                                                     >  2:3639140/1‑150 (MQ=255)
                                                                     catcatTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACaa                                                                                  >  2:3312379/1‑150 (MQ=255)
                                                                        tatTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACg                                                                               <  1:152194/149‑1 (MQ=255)
                                                                          ttttttCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCt                                                                             <  2:378965/150‑1 (MQ=255)
                                                                            ttttCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTAc                                                                           <  1:4258117/150‑1 (MQ=255)
                                                                               tCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACaaa                                                                        <  1:1807968/150‑1 (MQ=255)
                                                                               tCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACaaa                                                                        >  1:273967/1‑150 (MQ=255)
                                                                                cGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAAc                                                                       <  1:3968829/150‑1 (MQ=255)
                                                                                  gTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTTTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAAccc                                                                     >  2:2797221/1‑150 (MQ=255)
                                                                                  gTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAAccc                                                                     <  2:2042213/150‑1 (MQ=255)
                                                                                    cTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGa                                                                   >  2:3129949/1‑150 (MQ=255)
                                                                                     tGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGAc                                                                  <  2:5142966/150‑1 (MQ=255)
                                                                                      ggCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACa                                                                 >  1:4493202/1‑150 (MQ=255)
                                                                                      ggCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACa                                                                 >  2:2656646/1‑150 (MQ=255)
                                                                                        cTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAAc                                                               >  1:3005939/1‑150 (MQ=255)
                                                                                         ttCAGTCAGTATGGTACTGGGCGATCATCGAAAGANNGNGNANTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAAcc                                                              >  1:484312/1‑150 (MQ=255)
                                                                                         ttCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAAcc                                                              >  2:4989998/1‑150 (MQ=255)
                                                                                           cagtcagtATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGa                                                            <  2:1326861/150‑1 (MQ=255)
                                                                                           cagtcagtATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGa                                                            <  1:4582931/150‑1 (MQ=255)
                                                                                             gtcagtATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGAcc                                                          <  1:4215391/150‑1 (MQ=255)
                                                                                                       tACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCAcc                                                >  1:3331459/1‑150 (MQ=255)
                                                                                                       tACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCAcc                                                <  1:1259302/150‑1 (MQ=255)
                                                                                                           gggCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGc                                            >  2:1915577/1‑150 (MQ=255)
                                                                                                            tgCGATCATCGATAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGACTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCa                                           <  2:2747233/149‑1 (MQ=255)
                                                                                                             gCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAg                                          >  1:4715865/1‑150 (MQ=255)
                                                                                                                   atcGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATc                                    >  2:1545744/1‑150 (MQ=255)
                                                                                                                      gAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGt                                 >  1:460210/1‑150 (MQ=255)
                                                                                                                       aaaGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGtt                                >  1:5139386/1‑150 (MQ=255)
                                                                                                                       aaaGAACGAGAAGTTCAGCCTGTGGTGCAGAGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGtt                                <  2:376695/150‑1 (MQ=255)
                                                                                                                        aaGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGttt                               >  2:4016516/1‑150 (MQ=255)
                                                                                                                             cGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTg                          >  2:2232271/1‑150 (MQ=255)
                                                                                                                             cGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGACACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTg                          <  1:4362291/150‑1 (MQ=255)
                                                                                                                               agaAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGc                        <  1:3082643/150‑1 (MQ=255)
                                                                                                                               agaAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGc                        >  2:2906015/1‑150 (MQ=255)
                                                                                                                                 aaGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCaa                      >  2:4808889/1‑150 (MQ=255)
                                                                                                                                   gTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAAcg                    <  2:1775918/150‑1 (MQ=255)
                                                                                                                                   gTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAAcg                    >  1:792665/1‑150 (MQ=255)
                                                                                                                                    ttCAGCCTGTGGTGCAGAGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAAcgc                   >  2:4059660/1‑150 (MQ=255)
                                                                                                                                     tCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAAcgcg                  <  2:185778/150‑1 (MQ=255)
                                                                                                                                      cAGCCTGTGGTGCAGCGATCTTTCTGCCTCACAGTGTGCGTGTCTGCCACGCAGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGa                 <  2:3264719/150‑1 (MQ=255)
                                                                                                                                       aGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGAc                >  1:3264075/1‑150 (MQ=255)
                                                                                                                                       aGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGAc                >  1:2105904/1‑150 (MQ=255)
                                                                                                                                        gCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCTAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACAGACCTTGACGCACCGTGCAGAAGATCTGATTCTCTGGCAACGCGAcc               >  2:3876977/1‑150 (MQ=255)
                                                                                                                                         ccTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGAccc              >  2:4219830/1‑150 (MQ=255)
                                                                                                                                         acTGTGGTGCAGCGATCTATCTGCCTCACATAGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGAAGATGATCTGTTTCTCAGGCTACGCGAccc              >  2:2511967/2‑150 (MQ=255)
                                                                                                                                          cTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCg             <  1:5216647/150‑1 (MQ=255)
                                                                                                                                            gtggtgCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGtt           <  1:1128952/150‑1 (MQ=255)
                                                                                                                                             tggtgCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGtta          >  1:2161204/1‑150 (MQ=255)
                                                                                                                                               gtgCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTGTCTCTGGCAACGCGACCCGttatt        >  1:1879379/1‑150 (MQ=255)
                                                                                                                                                 gCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGTTATTAc      >  2:4851678/1‑150 (MQ=255)
                                                                                                                                                 gCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGTTATTAc      >  1:839863/1‑150 (MQ=255)
                                                                                                                                                  cAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGTTATTAct     >  2:4076155/1‑150 (MQ=255)
                                                                                                                                                  cAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGTTATTAct     >  2:820047/1‑150 (MQ=255)
                                                                                                                                                     cGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGTTATTACTCtt  <  2:3005954/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATATTACTAATGCAGTGCTGTCACATCTTGATCTCGACGATCTGATCGCTGACGTCGCTCGTGAGATCCATCATTTTTTCGGTCTGGCTTCAGTCAGTATGGTACTGGGCGATCATCGAAAGAACGAGAAGTTCAGCCTGTGGTGCAGCGATCTTTCTGCCTCACATTGTGCGTGTCTGCCACGCTGTATGCCTGGCGAAAGTGTATTGCTGACACAAACGCTACAAACCCGACAACCGACCTTGACGCACCGTGCAGATGATCTGTTTCTCTGGCAACGCGACCCGTTATTACTCTT  >  CP009273/2605733‑2606031

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: