Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,679,816 C→A 6.8% D344Y (GAC→TAC)  fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,679,8160CA 6.8% 1106.3 / 39.3 292D344Y (GAC→TAC) fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (155/117);  new base A (13/7);  total (168/124)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 313 total aligned reads displayed.

CGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACAT  >  CP009273/1679667‑1679965
                                                                                                                                                     |                                                                                                                                                     
cGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCAcgtc                                                                                                                                                       <  1:4144348/150‑1 (MQ=255)
 ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                                                                      <  1:4464386/150‑1 (MQ=255)
 ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                                                                      >  1:140702/1‑150 (MQ=255)
 ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                                                                      <  1:2082060/150‑1 (MQ=255)
 ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCAAGTTGATCGCCAcgtcg                                                                                                                                                      >  2:4228727/1‑150 (MQ=255)
  gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCAcgtcgt                                                                                                                                                     >  2:1422396/1‑150 (MQ=255)
   ttCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGtt                                                                                                                                                    <  1:2697982/150‑1 (MQ=255)
    tCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTc                                                                                                                                                   <  2:861508/150‑1 (MQ=255)
       aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGAGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                                                                >  1:3249633/1‑150 (MQ=255)
         aCCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                                                              >  1:4067837/1‑150 (MQ=255)
           ccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                                                            >  1:918458/1‑150 (MQ=255)
           ccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                                                            <  1:4734961/150‑1 (MQ=255)
            cACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCAc                                                                                                                                           >  2:3672481/1‑150 (MQ=255)
               tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTg                                                                                                                                        >  2:2617096/1‑150 (MQ=255)
                gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                                                                       <  1:859852/150‑1 (MQ=255)
                 cgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                                                                                                      >  1:1284682/1‑150 (MQ=255)
                    aGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                                                                                                   >  1:1915059/1‑150 (MQ=255)
                     gtgtTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                                                                                                                  <  2:2433035/150‑1 (MQ=255)
                      tgtttgttAAAACTTTCCATGCCTTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagcag                                                                                                                                 <  1:2757073/150‑1 (MQ=255)
                      tgtttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagcag                                                                                                                                 <  2:3599140/150‑1 (MQ=255)
                        tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCagag                                                                                                                               >  2:751530/1‑150 (MQ=255)
                            ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                                                                                                           <  1:1517017/150‑1 (MQ=255)
                              aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTAGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                                                                         >  1:2794480/1‑150 (MQ=255)
                               aaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGtt                                                                                                                        <  1:507777/150‑1 (MQ=255)
                               aaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGtt                                                                                                                        >  1:3775538/1‑150 (MQ=255)
                                 aCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCAGTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTAGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                                                                      >  1:2141768/1‑150 (MQ=255)
                                   tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATg                                                                                                                    <  2:2718717/150‑1 (MQ=255)
                                    ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                                                                   >  2:604637/1‑150 (MQ=255)
                                    ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                                                                   >  2:3433454/1‑150 (MQ=255)
                                     tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                                                                  <  2:3179257/150‑1 (MQ=255)
                                     tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTAGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                                                                  <  2:2311695/150‑1 (MQ=255)
                                      ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                                                                 >  2:1218356/1‑150 (MQ=255)
                                          gCCATCTGCCATCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                                                             <  2:3774230/150‑1 (MQ=255)
                                           ccATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCacac                                                                                                            >  1:2368808/1‑150 (MQ=255)
                                            cATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                                                                                                           <  2:3362810/150‑1 (MQ=255)
                                            cATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAGGGTTAATGCCTCACACt                                                                                                           >  2:4338263/1‑150 (MQ=255)
                                             aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtg                                                                                                          >  2:3062939/1‑150 (MQ=255)
                                              tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTAGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtgt                                                                                                         >  2:1588358/1‑150 (MQ=255)
                                                tGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTgtgt                                                                                                       >  2:4465889/1‑150 (MQ=255)
                                                 gCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                                                                      <  1:389992/150‑1 (MQ=255)
                                                   cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                                                                    >  1:765341/1‑150 (MQ=255)
                                                    aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTggg                                                                                                   >  2:2127143/1‑150 (MQ=255)
                                                      cAAGCGCACCGCTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                                                                                                 >  1:3983337/1‑150 (MQ=255)
                                                      cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                                                                                                 >  1:364962/1‑150 (MQ=255)
                                                        aGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCa                                                                                               <  1:3634818/150‑1 (MQ=255)
                                                          cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACt                                                                                             >  2:198642/1‑150 (MQ=255)
                                                                gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCACAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcccc                                                                                       <  2:2708896/150‑1 (MQ=255)
                                                                gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcccc                                                                                       >  2:526938/1‑150 (MQ=255)
                                                                gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTAGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcccc                                                                                       <  2:4717139/150‑1 (MQ=255)
                                                                 aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCg                                                                                      >  2:4003950/1‑150 (MQ=255)
                                                                 aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCg                                                                                      >  2:4181752/1‑150 (MQ=255)
                                                                     cAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCAt                                                                                  >  1:4623483/1‑150 (MQ=255)
                                                                     cAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCAt                                                                                  >  1:3130179/1‑150 (MQ=255)
                                                                     cAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCAt                                                                                  >  2:4579069/1‑150 (MQ=255)
                                                                       gtaaaTTGTTGATCACCATTGGACGGAAGACTTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACATCGTTCCCCATCACCTGACAGCAGATCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                                                                <  1:1001685/148‑1 (MQ=255)
                                                                        gAAATTGTGGATCACCATTGGACGGAAGACTTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGAt                                                                               <  1:2589579/150‑1 (MQ=255)
                                                                         aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGAtt                                                                              >  2:1282163/1‑150 (MQ=255)
                                                                            ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAg                                                                           >  1:4314026/1‑150 (MQ=255)
                                                                            ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAg                                                                           >  1:4313353/1‑150 (MQ=255)
                                                                             tgtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGc                                                                          <  1:4643500/150‑1 (MQ=255)
                                                                               tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTg                                                                        >  2:9216/1‑150 (MQ=255)
                                                                               tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTg                                                                        >  2:1677778/1‑150 (MQ=255)
                                                                                ggATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGc                                                                       <  2:769269/150‑1 (MQ=255)
                                                                                  aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGccc                                                                     <  1:383427/150‑1 (MQ=255)
                                                                                       cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                                                                <  2:891379/150‑1 (MQ=255)
                                                                                       cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                                                                >  1:1224995/1‑150 (MQ=255)
                                                                                       cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                                                                >  2:3340113/1‑150 (MQ=255)
                                                                                              cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTcc                                                         <  2:135718/150‑1 (MQ=255)
                                                                                                 aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCggg                                                      <  1:1809613/150‑1 (MQ=255)
                                                                                                  aGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGa                                                     <  1:142867/150‑1 (MQ=255)
                                                                                                      gttcagttcaAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga                                                 <  1:3039804/150‑1 (MQ=255)
                                                                                                       ttcagttcaAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTgag                                                <  2:3249633/150‑1 (MQ=255)
                                                                                                         cagttcaAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt                                              <  1:3179066/150‑1 (MQ=255)
                                                                                                          agttcaAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAtt                                             >  1:1146358/1‑150 (MQ=255)
                                                                                                          agttcaAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAtt                                             >  1:750699/1‑150 (MQ=255)
                                                                                                            ttcaAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTc                                           >  1:2269138/1‑150 (MQ=255)
                                                                                                                aaGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTGTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCa                                       <  1:887047/150‑1 (MQ=255)
                                                                                                                   ttACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAAtt                                    >  2:1756689/1‑150 (MQ=255)
                                                                                                                   ttACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTAGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAAtt                                    <  2:4605410/150‑1 (MQ=255)
                                                                                                                     aCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTcc                                  <  1:198642/150‑1 (MQ=255)
                                                                                                                     aCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTcc                                  >  1:2009585/1‑150 (MQ=255)
                                                                                                                      ccGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCg                                 <  2:4202095/150‑1 (MQ=255)
                                                                                                                       cGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGc                                >  1:4221679/1‑150 (MQ=255)
                                                                                                                       cGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGc                                <  1:4085838/150‑1 (MQ=255)
                                                                                                                        ggAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCa                               <  2:3886387/150‑1 (MQ=255)
                                                                                                                             cgcCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGcgcg                          <  1:2878466/150‑1 (MQ=255)
                                                                                                                              gcCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCgg                         >  1:4579751/1‑150 (MQ=255)
                                                                                                                               ccccccATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCggg                        <  1:3433454/150‑1 (MQ=255)
                                                                                                                                 ccccATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGcc                      <  1:2477175/150‑1 (MQ=255)
                                                                                                                                   ccATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCag                    >  2:2897810/1‑150 (MQ=255)
                                                                                                                                    cATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCaga                   >  1:1644495/1‑150 (MQ=255)
                                                                                                                                    cATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCaga                   >  1:3728980/1‑150 (MQ=255)
                                                                                                                                     aTGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCagag                  >  2:2687313/1‑150 (MQ=255)
                                                                                                                                     aTGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAGTGAGATTTCACCAATTCCGCAGCGCGGGCCagag                  >  1:206493/1‑150 (MQ=255)
                                                                                                                                         tGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGgcca              <  2:4097806/150‑1 (MQ=255)
                                                                                                                                            tCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGgccagcc           >  2:693500/1‑150 (MQ=255)
                                                                                                                                             cGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGgccagcca          <  1:4568426/150‑1 (MQ=255)
                                                                                                                                                   gtcgtTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCCGCGGAc    >  2:3820227/1‑150 (MQ=255)
                                                                                                                                                     cgtTCCCCATCACCTGACAGCAGAGCATGGTTCATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGAcat  <  2:2504406/150‑1 (MQ=255)
                                                                                                                                                     cgtTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGAcat  <  1:1177136/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACAT  >  CP009273/1679667‑1679965

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: