Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,176,513 C→T 25.7% A206A (GCG→GCA yegX ← putative family 25 glycosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,176,5130CT25.7% 439.9 / 165.6 229A206A (GCG→GCAyegXputative family 25 glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base C (78/92);  new base T (35/24);  total (113/116)
Fisher's exact test for biased strand distribution p-value = 9.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 232 total aligned reads displayed.

CATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTACGTCCAGCACGGCAGG  >  CP009273/2176365‑2176662
                                                                                                                                                    |                                                                                                                                                     
   tAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCa                                                                                                                                                   <  1:4566586/150‑1 (MQ=255)
   tAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCa                                                                                                                                                   <  1:1536868/150‑1 (MQ=255)
    aaaTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAg                                                                                                                                                  >  1:1288290/1‑150 (MQ=255)
    aaaTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAg                                                                                                                                                  >  1:2415109/1‑150 (MQ=255)
     aaTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGa                                                                                                                                                 >  1:3931840/1‑150 (MQ=255)
       tACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGAtt                                                                                                                                               <  1:630290/150‑1 (MQ=255)
       tACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGAtt                                                                                                                                               <  1:481697/150‑1 (MQ=255)
       tACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGAtt                                                                                                                                               <  2:4727220/150‑1 (MQ=255)
        aCATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTg                                                                                                                                              >  1:4042267/1‑150 (MQ=255)
        aCATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTg                                                                                                                                              <  2:2787360/150‑1 (MQ=255)
          aTTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGgt                                                                                                                                            <  1:3651236/150‑1 (MQ=255)
          aTTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGgt                                                                                                                                            >  2:439321/1‑150 (MQ=255)
           ttAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGgtg                                                                                                                                           >  2:3660943/1‑150 (MQ=255)
           ttAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGgtg                                                                                                                                           <  2:3999089/150‑1 (MQ=255)
            tAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATCTTCATTGAAATAGCCCGCCAGAGTGgtgt                                                                                                                                          >  2:3878551/1‑150 (MQ=255)
            tAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGgtgt                                                                                                                                          >  2:2456099/1‑150 (MQ=255)
            tAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGgtgt                                                                                                                                          <  2:1682853/150‑1 (MQ=255)
             aaaaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGtgtg                                                                                                                                         >  2:3325579/1‑150 (MQ=255)
             aaaaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGtgtg                                                                                                                                         <  1:3829305/150‑1 (MQ=255)
               aaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGAt                                                                                                                                       <  1:4345167/150‑1 (MQ=255)
                agccACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACTTTGATAGTAGTGAGCCACCCACCACGGTTATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                                                                                                                      <  2:2686945/148‑1 (MQ=255)
                aTCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                                                                                                                      >  1:3330536/1‑150 (MQ=255)
                aTCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                                                                                                                      <  2:2382631/150‑1 (MQ=255)
                aTCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                                                                                                                      <  1:2901461/150‑1 (MQ=255)
                aTCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                                                                                                                      <  1:2413123/150‑1 (MQ=255)
                 tCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaa                                                                                                                                     <  2:3750396/150‑1 (MQ=255)
                  ccACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATaaa                                                                                                                                    >  2:476607/1‑150 (MQ=255)
                     ccGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaaaa                                                                                                                                 >  2:3205333/1‑150 (MQ=255)
                       ggACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACg                                                                                                                               >  1:1134835/1‑150 (MQ=255)
                        gACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACgg                                                                                                                              >  2:3519774/1‑150 (MQ=255)
                         aCCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGc                                                                                                                             >  1:2255752/1‑150 (MQ=255)
                             ttGATGCCATCTACCTGTCCACGTTCGTAATGCTGCCAGAAGCGCCAGGCAATTCCTTATTTGTCCGTACTACTTTTATAGTAGTGAGCCACCCACCACGTATATTCATTTAAATAGCCCGCATTATTGTTTTTATAAAAAACGGCTCCt                                                                                                                         <  2:2983839/150‑1 (MQ=255)
                               gATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGa                                                                                                                       >  1:1874660/1‑150 (MQ=255)
                                aTGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAg                                                                                                                      >  1:279775/1‑150 (MQ=255)
                                     aTCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaat                                                                                                                 >  1:1361109/1‑150 (MQ=255)
                                     aTCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaat                                                                                                                 >  1:887516/1‑150 (MQ=255)
                                     aTCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaat                                                                                                                 <  2:1456467/150‑1 (MQ=255)
                                     aTCTACCGGTCCACGGTCGGAATGCTGCCCGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTGGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaat                                                                                                                 <  2:1454321/150‑1 (MQ=255)
                                      tCTACCTGGCCACGGTCGGAATGCTGCCAGAAGCGCCCGGCCATTCCGTCTTTGTCCGGACGAGGTTGATAGTCGTGACCCACCCACCACGGATATTATTTTAAATTGCCCGCCATATGGGTGTGATAAAAAACGGCTCCTGAGTAaata                                                                                                                <  1:3150217/150‑1 (MQ=255)
                                       cTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTTTTTGTAGTGAGCCACCCACCAGGTATATTCTTTTAAATAGCCCGCCATTTTGGTGTGATAAAAAACGGCTCCTGAGTAaataa                                                                                                               <  1:4696196/150‑1 (MQ=255)
                                        tACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaataat                                                                                                              >  1:2034711/1‑150 (MQ=255)
                                             gTCCACGGTCGGAATGCTGCCAGACGCGCCCGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGTATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCt                                                                                                         <  2:3033190/150‑1 (MQ=255)
                                               ccACGGTCGGAATGCTGCCAGAAGCGCCCGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCttt                                                                                                       >  2:2667901/1‑150 (MQ=255)
                                                ctcGGTCGGATTGCTGCCAGAAGCGCCAGGCCTTTCCTTCTTTGTCCGGACTCCTTTGATTGTAGTGACCCACCCACCACGGATTTTAATTTAAATAGCCTGCCTGATTGGTTTTATAAAAAACGGCTCCTGAGTAAATAATCGGCtttt                                                                                                      <  1:3405243/148‑1 (MQ=255)
                                                cACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCtttt                                                                                                      <  2:3964845/150‑1 (MQ=255)
                                                cACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCtttt                                                                                                      >  2:2286589/1‑150 (MQ=255)
                                                cACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCtttt                                                                                                      >  2:4737236/1‑150 (MQ=255)
                                                 aCGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCttttt                                                                                                     <  2:4759892/150‑1 (MQ=255)
                                                  cGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTc                                                                                                    <  2:2325869/150‑1 (MQ=255)
                                                   ggTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTcc                                                                                                   >  2:2562898/1‑150 (MQ=255)
                                                    gtgggATTGCTGCCAGACGCGCCTGGCCATTCGGTATTTGTCGGGACGACGTTGATTGTCGTGAGCCACCCACAAGGTTTTTTCATTGAAATTGCCCGCCATTTGGGTGTTATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTccc                                                                                                  <  1:906883/147‑1 (MQ=255)
                                                    gTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTccc                                                                                                  >  1:257803/1‑150 (MQ=255)
                                                     tCGTAAGGCTGCCAGACGCGCCAGGCCATTCCGAATTTGTCCGGACGAGGTTGTTAGTAGTGAGCCACCCCCCAGGGATATTCTTTGAAATAGCCCGCCTTATTGGTTTTATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTcccc                                                                                                 <  1:3906078/150‑1 (MQ=255)
                                                        gAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTa                                                                                              >  2:3330889/1‑150 (MQ=255)
                                                        gAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTa                                                                                              >  2:1583645/1‑150 (MQ=255)
                                                         aaTGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTAc                                                                                             >  1:3745120/1‑150 (MQ=255)
                                                           tgctgcCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACtt                                                                                           >  2:4454598/1‑150 (MQ=255)
                                                           tgctgcCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACtt                                                                                           <  1:4016861/150‑1 (MQ=255)
                                                            gttgcCATAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACTTTTATAGTAGTGAGCCACCCACCACGTATATACATTGAAATAGCCTGCCAGATTGGTTTTATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACttt                                                                                          <  1:2785535/148‑1 (MQ=255)
                                                             ctgcCAGAAGCGCACGGCCATTCCGTAATGGTCCGGACGACTTTGATTGTCGTGAGCCACCCACCACGGATTTTAATTGAAATCGCCCGCCATTTTGGTGTTATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACtttt                                                                                         <  1:974287/150‑1 (MQ=255)
                                                              tgcCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACttttt                                                                                        <  2:659516/150‑1 (MQ=255)
                                                               gcCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACtttttt                                                                                       >  2:1956150/1‑150 (MQ=255)
                                                                    aaGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGAcc                                                                                  >  2:1247065/1‑150 (MQ=255)
                                                                     aGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCa                                                                                 <  2:3058974/150‑1 (MQ=255)
                                                                       cgcCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTTGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCAtt                                                                               <  2:584481/150‑1 (MQ=255)
                                                                          cAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCAttttt                                                                            >  2:4832520/1‑150 (MQ=255)
                                                                          cAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCAttttt                                                                            >  1:1090631/1‑150 (MQ=255)
                                                                          cAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCAttttt                                                                            >  2:2520024/1‑150 (MQ=255)
                                                                           aGGCCATGCCGTCTTTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTa                                                                           <  2:4168756/150‑1 (MQ=255)
                                                                            ggCCATTCCGTCATTGACCGGACGACGTTGATAGTAGTGAGCCACCAACCACGGAAATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAg                                                                          >  2:1655328/1‑150 (MQ=255)
                                                                              ccATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGcc                                                                        <  2:1512423/150‑1 (MQ=255)
                                                                                aTTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCAc                                                                      >  1:3292269/1‑150 (MQ=255)
                                                                                   ccGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGa                                                                   >  2:539222/1‑150 (MQ=255)
                                                                                             ccGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCt                                                         >  1:703894/1‑150 (MQ=255)
                                                                                              cGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCtt                                                        <  2:4667079/150‑1 (MQ=255)
                                                                                                   gacgTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTa                                                   >  2:3922990/1‑150 (MQ=255)
                                                                                                    acgTTGATAGTCGTGAGCCACCCACCACGGATATTCATTGAATTAGCCCGCCATTTTGGTGTTATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTaa                                                  <  1:4174218/150‑1 (MQ=255)
                                                                                                    acgTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTaa                                                  >  1:3568334/1‑150 (MQ=255)
                                                                                                      gTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAAtt                                                >  2:4390966/1‑150 (MQ=255)
                                                                                                      gTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAAtt                                                >  2:878863/1‑150 (MQ=255)
                                                                                                         tatatTAGTGAGCCACCCACCACTGATATTAATTTAAATAGCCCGCCAGTTTGGTGTTATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCtt                                             <  1:1250647/149‑1 (MQ=255)
                                                                                                         gATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCtt                                             >  2:4053163/1‑150 (MQ=255)
                                                                                                         gATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCAGTTTTAGCCACTGACTTACCCGCTTGCGTAATTCtt                                             >  2:4054980/1‑150 (MQ=255)
                                                                                                          aTAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCttt                                            >  1:4764263/1‑150 (MQ=255)
                                                                                                           tagtagTGAGCCACCCACCCCGGATATTCTTTGAAATAGCCCGCCAGATGGGTGTTATAAAACACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCtttt                                           <  2:836327/150‑1 (MQ=255)
                                                                                                               agtGACCCACCCACCACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCg                                       <  1:660445/150‑1 (MQ=255)
                                                                                                                   aGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATaa                                   >  1:4179816/1‑150 (MQ=255)
                                                                                                                      cacccaccACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAAttt                                >  2:701623/1‑150 (MQ=255)
                                                                                                                       acccaccACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAAtttt                               >  2:2232093/1‑150 (MQ=255)
                                                                                                                       acccaccACGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAAtttt                               >  2:2232081/1‑150 (MQ=255)
                                                                                                                         ccaccaCGGATATTCATTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTcc                             >  1:1939925/1‑150 (MQ=255)
                                                                                                                            ccacGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCg                          >  1:2939203/1‑150 (MQ=255)
                                                                                                                                ggATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGttct                      >  2:1358012/1‑150 (MQ=255)
                                                                                                                                    atTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTa                  <  1:3205333/150‑1 (MQ=255)
                                                                                                                                    atTAATTTAAATAGCCCGCCAGTTTGGTTTTATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTa                  <  1:1003957/150‑1 (MQ=255)
                                                                                                                                     tactttGAAATAGCCCGCCAGTTGGGTGTAATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTAc                 <  1:439321/146‑1 (MQ=255)
                                                                                                                                        aTTGAAATAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTACGTc              >  2:2177987/1‑150 (MQ=255)
                                                                                                                                               tAGCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTACGTCCAGCACg       <  2:1134835/150‑1 (MQ=255)
                                                                                                                                                ggCCTGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTACGTCCAGCACgg      >  1:1093586/2‑150 (MQ=255)
                                                                                                                                                    cGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTACGTCCAGCACGGCAgg  >  1:1666698/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTTCCCCGTACTTTTTTCGACCATTTTTAGCCACTGACTTACCCGCTTGCGTAATTCTTTTGCCGATAATTTTCCCCGTTCTTCTACGTCCAGCACGGCAGG  >  CP009273/2176365‑2176662

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