Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 826,737 +GT coding (177/285 nt) ihfB → integration host factor (IHF), DNA‑binding protein, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306826,7351.G94.2% 1190.8 / 59.4 311coding (175/285 nt)ihfBintegration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base . (10/8);  new base G (113/180);  total (123/188)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01
*AP012306826,7352.T94.2% 1201.2 / 55.9 311coding (175/285 nt)ihfBintegration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base . (10/8);  new base T (113/180);  total (123/188)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 323 total aligned reads displayed.

ACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑GTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAA  >  AP012306/826587‑826875
                                                                                                                                                     ||                                                                                                                                            
aCTTGCCACCCAGCAATCGCACATCCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑g                                                                                                                                             <  2:2186201/150‑1 (MQ=255)
 tttGCCACCGAGAAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gt                                                                                                                                            <  2:403147/149‑1 (MQ=255)
 cTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gt                                                                                                                                            <  2:745555/150‑1 (MQ=255)
    gCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgta                                                                                                                                         <  2:2219771/150‑4 (MQ=255)
     ccACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgtac                                                                                                                                        <  1:3685254/150‑5 (MQ=255)
      cACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCCGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCAAGCGTATTGACATCCGCGGCTTCGGCAGCTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgtacc                                                                                                                                       >  2:3986105/1‑145 (MQ=255)
      cACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgtacc                                                                                                                                       >  1:2534775/1‑145 (MQ=255)
        cccAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgtaccgg                                                                                                                                     >  2:516359/1‑143 (MQ=255)
        cccAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgtaccgg                                                                                                                                     >  2:1162869/1‑143 (MQ=255)
           aGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACg                                                                                                                                    >  1:418856/1‑150 (MQ=255)
             cAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTa                                                                                                                                  <  1:4002126/150‑1 (MQ=255)
             cAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTa                                                                                                                                  >  1:2734574/1‑150 (MQ=255)
              aaTCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTaa                                                                                                                                 <  2:2728813/150‑1 (MQ=255)
              aaTCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTaa                                                                                                                                 <  1:3964200/150‑1 (MQ=255)
               aTCGCACATCCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAAt                                                                                                                                <  1:2270043/150‑1 (MQ=255)
                 cGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATcc                                                                                                                              <  1:116194/150‑1 (MQ=255)
                 cGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATcc                                                                                                                              <  2:3938146/150‑1 (MQ=255)
                  gCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCg                                                                                                                             >  2:1266713/1‑150 (MQ=255)
                       ttCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGAc                                                                                                                        <  2:903500/150‑1 (MQ=255)
                       tccccGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGAc                                                                                                                        <  1:1546879/148‑1 (MQ=255)
                        tCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACt                                                                                                                       <  1:2895906/150‑1 (MQ=255)
                          ccGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACATAATCCGAAGACTgg                                                                                                                     <  1:243966/150‑1 (MQ=255)
                           cGCCAAGACGGTTGAAGATGCAGTAGAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGc                                                                                                                    >  1:1469877/1‑150 (MQ=255)
                           cGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGc                                                                                                                    >  1:3741735/1‑150 (MQ=255)
                                 gACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATaaa                                                                                                              <  1:224592/150‑1 (MQ=255)
                                  aCGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAg                                                                                                             <  1:3121867/150‑1 (MQ=255)
                                   cGGTTGAAGATGCAGTAAAAGAGATGCGGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCAAAGACTGGCGATAAAGt                                                                                                            >  2:1771426/1‑150 (MQ=255)
                                      ttGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGa                                                                                                         <  2:617067/150‑1 (MQ=255)
                                       tGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGaa                                                                                                        >  1:3535505/1‑150 (MQ=255)
                                        gAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAAc                                                                                                       <  2:2384747/150‑1 (MQ=255)
                                          aGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTg                                                                                                     <  2:3546564/150‑1 (MQ=255)
                                          aGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTg                                                                                                     <  2:2285367/150‑1 (MQ=255)
                                           gATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTgg                                                                                                    <  2:2261009/150‑1 (MQ=255)
                                            aTGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTgga                                                                                                   <  2:3048653/150‑1 (MQ=255)
                                              gCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTggaag                                                                                                 <  1:11983/150‑1 (MQ=255)
                                              gCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTggaag                                                                                                 <  2:1291280/150‑1 (MQ=255)
                                                 gTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTggaaggaa                                                                                              >  1:2631205/1‑150 (MQ=255)
                                                   aaaaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGaaaa                                                                                            <  2:1933834/150‑1 (MQ=255)
                                                   aaaaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGaaaa                                                                                            >  2:1483552/1‑150 (MQ=255)
                                                    aaaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAAt                                                                                           <  1:3923234/150‑1 (MQ=255)
                                                     aaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATa                                                                                          <  1:4149986/150‑1 (MQ=255)
                                                       gagaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACg                                                                                        <  1:232364/150‑1 (MQ=255)
                                                       gagaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACg                                                                                        <  2:784583/150‑1 (MQ=255)
                                                       gagaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACg                                                                                        <  1:3494260/150‑1 (MQ=255)
                                                        agaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGt                                                                                       >  2:1349596/1‑150 (MQ=255)
                                                         gaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGtt                                                                                      <  2:29699/150‑1 (MQ=255)
                                                         gaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGtt                                                                                      <  2:204635/150‑1 (MQ=255)
                                                          aTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTc                                                                                     <  1:2599492/150‑1 (MQ=255)
                                                          aTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTc                                                                                     >  2:1023619/1‑150 (MQ=255)
                                                           tGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTcc                                                                                    <  2:2431065/150‑1 (MQ=255)
                                                           tGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTcc                                                                                    <  2:55150/150‑1 (MQ=255)
                                                            gCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCt                                                                                   <  2:3417274/150‑1 (MQ=255)
                                                               ggAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCAc                                                                                <  1:2841521/150‑1 (MQ=255)
                                                               ggAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCAc                                                                                >  2:1962303/1‑150 (MQ=255)
                                                                 aGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACtt                                                                              >  2:3775652/1‑150 (MQ=255)
                                                                  gCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACttt                                                                             <  2:1869267/150‑1 (MQ=255)
                                                                   cATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTctttaaa                                                                            <  1:1797900/150‑7 (MQ=255)
                                                                    atatGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTaa                                                                           <  1:3795848/150‑1 (MQ=255)
                                                                     tatGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTaaa                                                                          <  2:4153602/150‑1 (MQ=255)
                                                                      atGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAAc                                                                         >  2:1592525/1‑150 (MQ=255)
                                                                         gCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTg                                                                      <  1:807287/150‑1 (MQ=255)
                                                                         gCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTg                                                                      >  1:35977/1‑150 (MQ=255)
                                                                          ccTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTgg                                                                     >  2:3799522/1‑150 (MQ=255)
                                                                            tCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTa                                                                   <  1:1112275/150‑1 (MQ=255)
                                                                             cGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTaa                                                                  >  1:2234571/1‑150 (MQ=255)
                                                                              gACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTaaa                                                                 <  1:3311215/150‑1 (MQ=255)
                                                                              gACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGCaaa                                                                 >  1:531251/1‑150 (MQ=255)
                                                                               aCTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAg                                                                >  2:3387043/1‑150 (MQ=255)
                                                                                     gcgcAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTgc                                                          <  2:1666750/150‑1 (MQ=255)
                                                                                          gggCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGcgcg                                                       >  2:722305/1‑148 (MQ=255)
                                                                                           ggCGTGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATc                                                    <  1:3520576/150‑1 (MQ=255)
                                                                                            gCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATcg                                                   <  2:1715601/150‑1 (MQ=255)
                                                                                            gCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATcg                                                   <  2:2346119/150‑1 (MQ=255)
                                                                                              gAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATcgcg                                                 >  1:3282527/1‑150 (MQ=255)
                                                                                              gAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATcgcg                                                 >  2:2127057/1‑150 (MQ=255)
                                                                                              gAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATcgcg                                                 >  1:1620586/1‑150 (MQ=255)
                                                                                               aGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCCCCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGTGCTATAgcgc                                                >  2:3188881/1‑150 (MQ=255)
                                                                                               aGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCgcgc                                                >  1:1104052/1‑150 (MQ=255)
                                                                                                   tATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAat                                            <  1:1502589/150‑1 (MQ=255)
                                                                                                   tATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAat                                            >  1:1705184/1‑150 (MQ=255)
                                                                                                    aTTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAata                                           <  2:4241773/150‑1 (MQ=255)
                                                                                                        aaaTCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTa                                       >  2:1231191/1‑150 (MQ=255)
                                                                                                                ggTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAg                               >  1:2007533/1‑150 (MQ=255)
                                                                                                                  tttCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGGAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtt                             >  1:1380959/1‑150 (MQ=255)
                                                                                                                   ttCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttt                            <  1:932495/150‑1 (MQ=255)
                                                                                                                     cGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttttt                          >  1:3207631/1‑150 (MQ=255)
                                                                                                                       gCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTa                        >  2:3626522/1‑150 (MQ=255)
                                                                                                                       gCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTa                        >  2:2531650/1‑150 (MQ=255)
                                                                                                                              ctctTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAAc                 <  1:577711/150‑1 (MQ=255)
                                                                                                                                ctTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACtt               <  1:3053250/150‑1 (MQ=255)
                                                                                                                                ctTTGCACTACCGCGCACCACGTGTACCGGACGGAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACtt               >  2:3320890/1‑150 (MQ=255)
                                                                                                                                  ttGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGa             <  2:149689/150‑1 (MQ=255)
                                                                                                                                  ttGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGa             <  1:4078178/150‑1 (MQ=255)
                                                                                                                                   tGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGaa            >  1:707670/1‑150 (MQ=255)
                                                                                                                                   tGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGaa            <  1:1319396/150‑1 (MQ=255)
                                                                                                                                      aCTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga         >  1:2925528/1‑150 (MQ=255)
                                                                                                                                       cTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgag        <  2:2678572/150‑1 (MQ=255)
                                                                                                                                        tACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgaga       >  1:1008550/1‑150 (MQ=255)
                                                                                                                                         aCCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagag      >  2:443299/1‑150 (MQ=255)
                                                                                                                                            gcgcACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGaaa   <  2:2821673/150‑1 (MQ=255)
                                                                                                                                                     ||                                                                                                                                            
ACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑GTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAA  >  AP012306/826587‑826875

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑