Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,680,642 C→A 91.4% E68D (GAG→GAT 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,680,6420CA91.4% 480.0 / 38.3 183E68D (GAG→GATfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base A (92/75);  total (99/84)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 187 total aligned reads displayed.

TAATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGC  >  CP009273/1680495‑1680786
                                                                                                                                                   |                                                                                                                                                
  agtCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTtatc                                                                                                                                                <  2:2012688/146‑1 (MQ=255)
       cTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAg                                                                                                                                           >  2:2262288/1‑148 (MQ=255)
          gCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGaca                                                                                                                                        >  1:2950039/1‑148 (MQ=255)
           cccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaa                                                                                                                                       >  1:3280147/1‑148 (MQ=255)
             cGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaaca                                                                                                                                     <  1:512669/148‑1 (MQ=255)
               gTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAg                                                                                                                                   <  2:2157691/148‑1 (MQ=255)
                   gCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTa                                                                                                                               >  1:3006377/1‑148 (MQ=255)
                      aGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAAt                                                                                                                            >  1:3152303/1‑148 (MQ=255)
                       gCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCAACCGCCGCCTGACGAATGGCGCTCGCTGTATCTTCAGACAACAAGCCTAAATc                                                                                                                           >  1:3307924/1‑148 (MQ=255)
                        cACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCt                                                                                                                          >  2:1263289/1‑148 (MQ=255)
                        cACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCt                                                                                                                          >  2:1263000/1‑148 (MQ=255)
                          cTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTc                                                                                                                        >  1:2138459/1‑148 (MQ=255)
                          cTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTc                                                                                                                        <  2:1071456/148‑1 (MQ=255)
                           ttCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCa                                                                                                                       <  1:2019124/148‑1 (MQ=255)
                              gttcatgttcatGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTa                                                                                                                    <  2:2759469/148‑1 (MQ=255)
                               ttcatgttcatgTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTaa                                                                                                                   <  2:2004017/148‑1 (MQ=255)
                               ttcatgttcatgTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTaa                                                                                                                   >  2:2176376/1‑148 (MQ=255)
                                tcatgttcatgtTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAAc                                                                                                                  >  2:111172/1‑148 (MQ=255)
                                  atgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACtt                                                                                                                >  1:238603/1‑148 (MQ=255)
                                    gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                              >  1:694240/1‑148 (MQ=255)
                                    gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                              <  1:3149570/148‑1 (MQ=255)
                                    gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                              >  1:1245431/1‑148 (MQ=255)
                                     ttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTg                                                                                                             <  2:1301154/148‑1 (MQ=255)
                                      tcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGc                                                                                                            <  2:2636956/148‑1 (MQ=255)
                                        atgttACTGTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCg                                                                                                          <  1:879629/148‑1 (MQ=255)
                                         tgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGc                                                                                                         <  1:1992803/148‑1 (MQ=255)
                                          gttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCt                                                                                                        >  1:1288500/1‑148 (MQ=255)
                                           ttACTTTGCGTGCAGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCAACCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTACCTTAACTTTTGCCTCTg                                                                                                       >  2:1265355/1‑148 (MQ=255)
                                               tttGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACg                                                                                                   >  1:1789091/1‑148 (MQ=255)
                                                ttGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGc                                                                                                  >  2:2347041/1‑148 (MQ=255)
                                                ttGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGc                                                                                                  <  1:2567352/148‑1 (MQ=255)
                                                ttGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGc                                                                                                  >  1:2736996/1‑148 (MQ=255)
                                                  gCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCtt                                                                                                <  2:514840/148‑1 (MQ=255)
                                                   cGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTg                                                                                               <  2:798492/148‑1 (MQ=255)
                                                     tGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCCCCACCCCAAAGCCCCTCGCTTAATCTTCAGACAACAACCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGt                                                                                             <  2:3157854/148‑1 (MQ=255)
                                                      gCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGtt                                                                                            <  1:2543402/148‑1 (MQ=255)
                                                        cGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAg                                                                                          >  2:2437231/1‑148 (MQ=255)
                                                        cGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAg                                                                                          <  2:2514673/148‑1 (MQ=255)
                                                        cGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAg                                                                                          >  2:779235/1‑148 (MQ=255)
                                                           aGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgcgc                                                                                       <  2:1405619/148‑1 (MQ=255)
                                                             ccGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCa                                                                                     >  2:1906217/1‑148 (MQ=255)
                                                                  cTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCa                                                                                >  1:626371/1‑148 (MQ=255)
                                                                  cTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCa                                                                                <  2:42721/148‑1 (MQ=255)
                                                                   tGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCAt                                                                               <  1:3197710/148‑1 (MQ=255)
                                                                     ccAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGa                                                                             <  1:1657176/148‑1 (MQ=255)
                                                                       aGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAAt                                                                           <  1:457299/148‑1 (MQ=255)
                                                                       aGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAAt                                                                           >  1:433923/1‑148 (MQ=255)
                                                                       aGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAAt                                                                           <  1:954411/148‑1 (MQ=255)
                                                                        gATAGCCAGCGGGAATTAGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATc                                                                          <  2:1973021/148‑1 (MQ=255)
                                                                         aTAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCa                                                                         >  2:1251989/1‑148 (MQ=255)
                                                                          tAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAg                                                                        >  2:2053545/1‑148 (MQ=255)
                                                                            gCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTg                                                                      <  2:2754538/148‑1 (MQ=255)
                                                                             ccAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGa                                                                     <  1:716673/148‑1 (MQ=255)
                                                                              cAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAg                                                                    <  2:384924/148‑1 (MQ=255)
                                                                               aGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAgg                                                                   >  1:2216351/1‑148 (MQ=255)
                                                                                gCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGt                                                                  <  2:149768/148‑1 (MQ=255)
                                                                                 cGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTg                                                                 <  2:1504357/148‑1 (MQ=255)
                                                                                    gAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGc                                                              >  2:1500694/1‑148 (MQ=255)
                                                                                      aTTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAt                                                            >  2:859981/1‑148 (MQ=255)
                                                                                      aTTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAt                                                            >  2:1090100/1‑148 (MQ=255)
                                                                                      aTTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAt                                                            <  1:2437231/148‑1 (MQ=255)
                                                                                      aTTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAt                                                            <  1:779200/148‑1 (MQ=255)
                                                                                       ttCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAtt                                                           >  2:2206220/1‑148 (MQ=255)
                                                                                        tcgtcgTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAttt                                                          >  2:1845445/1‑148 (MQ=255)
                                                                                         cgtcgtCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAtttt                                                         <  1:938348/148‑1 (MQ=255)
                                                                                          gtcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtc                                                        <  2:1929378/148‑1 (MQ=255)
                                                                                           tcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtct                                                       <  1:1065638/148‑1 (MQ=255)
                                                                                           tcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtct                                                       >  1:1086600/1‑148 (MQ=255)
                                                                                           tcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtct                                                       >  1:1503039/1‑148 (MQ=255)
                                                                                              tcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCg                                                    >  2:2049175/1‑148 (MQ=255)
                                                                                               cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                   >  2:763825/1‑148 (MQ=255)
                                                                                               cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                   >  1:846400/1‑148 (MQ=255)
                                                                                               cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                   >  1:1394477/1‑148 (MQ=255)
                                                                                               cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                   >  1:3085920/1‑148 (MQ=255)
                                                                                               cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                   >  1:1348412/1‑148 (MQ=255)
                                                                                               cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCAGTAACTTTTGCCGCTGCACGCGTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                   >  1:846397/1‑148 (MQ=255)
                                                                                               cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAAGCTTCATTATCTTTTGCCGCTGCACGCTGGGTTAGCGCCCGCGCATGAATCAGTGAGGTGGGCAGTGTCTCGGt                                                   >  1:2077867/1‑148 (MQ=255)
                                                                                                   cTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaa                                               <  1:1142841/148‑1 (MQ=255)
                                                                                                   cTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaa                                               <  2:635527/148‑1 (MQ=255)
                                                                                                   cTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaa                                               <  1:49349/148‑1 (MQ=255)
                                                                                                        cTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCg                                          <  1:2892736/148‑1 (MQ=255)
                                                                                                              gTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATg                                    <  1:3063610/148‑1 (MQ=255)
                                                                                                               tACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGc                                   >  1:2556677/1‑148 (MQ=255)
                                                                                                                  ttCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTcc                                >  2:1566749/1‑148 (MQ=255)
                                                                                                                  ttCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTcc                                <  1:2480130/148‑1 (MQ=255)
                                                                                                                  ttCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTcc                                >  1:2288499/1‑148 (MQ=255)
                                                                                                                    cATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCag                              >  2:1786831/1‑148 (MQ=255)
                                                                                                                           cgccTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGt                       >  1:2382438/1‑148 (MQ=255)
                                                                                                                            gccTGACGAATGGCGCTCGCTTTCTCTTCAGACAACACGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGACTCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGtt                      <  1:1421130/148‑1 (MQ=255)
                                                                                                                            gccTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGtt                      <  1:3199008/148‑1 (MQ=255)
                                                                                                                                gACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGt                  <  2:2948805/148‑1 (MQ=255)
                                                                                                                                gACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGt                  >  1:2393616/1‑148 (MQ=255)
                                                                                                                                    aaTGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtgt              >  2:1535997/1‑148 (MQ=255)
                                                                                                                                      tGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGAGGAAAAGCTCAAGCGAGCGTTGAGTTTTTgc            >  2:1694577/1‑148 (MQ=255)
                                                                                                                                             cgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCcac     >  2:816073/1‑148 (MQ=255)
                                                                                                                                              gctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCcaca    >  2:2720149/1‑148 (MQ=255)
                                                                                                                                              gctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCcaca    >  2:2720263/1‑148 (MQ=255)
                                                                                                                                               ctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAg   >  2:3284414/1‑148 (MQ=255)
                                                                                                                                               ctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAg   >  2:1323889/1‑148 (MQ=255)
                                                                                                                                                tttCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGc  >  1:1225313/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TAATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGC  >  CP009273/1680495‑1680786

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: