Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,680,642 C→A 84.9% 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,6420CA84.9% 638.3 / 124.5 279E68D (GAG→GATfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (32/10);  new base A (124/113);  total (156/123)
Fisher's exact test for biased strand distribution p-value = 4.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 284 total aligned reads displayed.

AATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCT  >  CP009273/1680496‑1680787
                                                                                                                                                  |                                                                                                                                                 
aaTTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTtatc                                                                                                                                                 <  2:4644195/149‑1 (MQ=255)
aaTTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTtatc                                                                                                                                                 >  2:3459288/1‑149 (MQ=255)
 aTTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTatct                                                                                                                                                >  1:4360922/1‑149 (MQ=255)
  ttCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCtt                                                                                                                                               >  1:1043714/1‑149 (MQ=255)
  ttCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCtt                                                                                                                                               >  1:1680803/1‑149 (MQ=255)
   tCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTc                                                                                                                                              <  2:2476925/149‑1 (MQ=255)
      cTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGa                                                                                                                                           >  2:891803/1‑149 (MQ=255)
      cTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGa                                                                                                                                           <  1:3673649/149‑1 (MQ=255)
          cccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaac                                                                                                                                       <  2:3162391/149‑1 (MQ=255)
           ccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaaca                                                                                                                                      >  1:3315045/1‑149 (MQ=255)
           ccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaaca                                                                                                                                      >  2:1892671/1‑149 (MQ=255)
           ccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaaca                                                                                                                                      >  1:4011961/1‑149 (MQ=255)
            cGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGAcaacaa                                                                                                                                     <  2:2236926/149‑1 (MQ=255)
              gTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGc                                                                                                                                   >  1:171799/1‑149 (MQ=255)
              gTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGc                                                                                                                                   >  2:1616845/1‑149 (MQ=255)
              gTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGc                                                                                                                                   >  1:3413749/1‑149 (MQ=255)
              gTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGc                                                                                                                                   >  2:4161271/1‑149 (MQ=255)
               ttAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGcc                                                                                                                                  >  2:4307650/1‑149 (MQ=255)
               ttAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGcc                                                                                                                                  <  2:4680851/149‑1 (MQ=255)
               ttAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGcc                                                                                                                                  <  1:449532/149‑1 (MQ=255)
                   ccAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTaaa                                                                                                                              <  1:4534473/149‑1 (MQ=255)
                     aGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATc                                                                                                                            >  1:405153/1‑149 (MQ=255)
                     aGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATc                                                                                                                            >  1:1198160/1‑149 (MQ=255)
                       cACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCtt                                                                                                                          <  2:3315045/149‑1 (MQ=255)
                        aCTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTc                                                                                                                         <  2:4554403/149‑1 (MQ=255)
                          ttCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCAt                                                                                                                       >  1:3497897/1‑149 (MQ=255)
                           tCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCAtt                                                                                                                      >  1:1570787/1‑149 (MQ=255)
                             gttcatgttcatGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTaa                                                                                                                    >  2:687954/1‑149 (MQ=255)
                             gttcatgttcatGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTaa                                                                                                                    >  2:229232/1‑149 (MQ=255)
                             gttcatgttcatGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTaa                                                                                                                    >  1:670840/1‑149 (MQ=255)
                                 atgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACttt                                                                                                                >  1:386964/1‑149 (MQ=255)
                                  tgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                               >  2:3042978/1‑149 (MQ=255)
                                   gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTg                                                                                                              <  2:4460012/149‑1 (MQ=255)
                                        tgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCt                                                                                                         >  1:547522/1‑149 (MQ=255)
                                         gttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTg                                                                                                        <  2:1711091/149‑1 (MQ=255)
                                          ttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGc                                                                                                       >  2:515644/1‑149 (MQ=255)
                                              tttGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGc                                                                                                   >  2:2746255/1‑149 (MQ=255)
                                               ttGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCt                                                                                                  <  1:2764233/149‑1 (MQ=255)
                                                 gCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTg                                                                                                >  1:2966174/1‑149 (MQ=255)
                                                 gCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTg                                                                                                <  1:599149/149‑1 (MQ=255)
                                                  cGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTgg                                                                                               <  2:2368782/149‑1 (MQ=255)
                                                  cGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTgg                                                                                               <  2:4764869/149‑1 (MQ=255)
                                                        ggAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgcg                                                                                         >  2:158265/1‑149 (MQ=255)
                                                        ggAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgcg                                                                                         >  1:3917950/1‑149 (MQ=255)
                                                        ggAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgcg                                                                                         >  1:2390340/1‑149 (MQ=255)
                                                         gAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgcgc                                                                                        >  2:132442/1‑149 (MQ=255)
                                                           gCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCa                                                                                      >  1:2891880/1‑149 (MQ=255)
                                                            ccGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAg                                                                                     >  1:1179519/1‑149 (MQ=255)
                                                               gTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAgcgc                                                                                  >  1:3083493/1‑149 (MQ=255)
                                                                  tGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATg                                                                               >  2:2915214/1‑149 (MQ=255)
                                                                    ccAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGac                                                                             <  1:2781648/149‑2 (MQ=255)
                                                                     cAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAAt                                                                            <  2:1469242/149‑1 (MQ=255)
                                                                     cAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAAt                                                                            >  1:1504250/1‑149 (MQ=255)
                                                                         tAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGt                                                                        <  2:1834583/149‑1 (MQ=255)
                                                                          aGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTg                                                                       >  2:4587804/1‑149 (MQ=255)
                                                                           gCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGa                                                                      >  1:2618976/1‑149 (MQ=255)
                                                                            ccAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAg                                                                     <  2:1220889/149‑1 (MQ=255)
                                                                               gCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTg                                                                  >  1:3438353/1‑149 (MQ=255)
                                                                                  ggAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGc                                                               <  1:225751/149‑1 (MQ=255)
                                                                                   gAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCa                                                              >  2:1758919/1‑149 (MQ=255)
                                                                                   gAATTCGTCGTCATGATGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGATTTATCTTCAGAAAACAAGCCTAAATCTTCAATAACTTTTGACGCTGAACGCTTGTTTATAGCCAGCGCATGAATCAGAGAGTTGGGCa                                                              >  1:3014174/1‑149 (MQ=255)
                                                                                          tcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                       <  2:3086610/149‑1 (MQ=255)
                                                                                          tcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                       >  2:1401843/1‑149 (MQ=255)
                                                                                          tcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                       <  2:810077/149‑1 (MQ=255)
                                                                                           cgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTcc                                                      <  1:1361373/149‑1 (MQ=255)
                                                                                             tcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                    >  1:3761924/1‑149 (MQ=255)
                                                                                              cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTc                                                   <  2:1432279/149‑1 (MQ=255)
                                                                                               aTGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCg                                                  >  1:1040363/1‑149 (MQ=255)
                                                                                                tGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGa                                                 >  1:3544069/1‑149 (MQ=255)
                                                                                                tGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGa                                                 >  2:1026663/1‑149 (MQ=255)
                                                                                                 gCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaa                                                <  2:386964/149‑1 (MQ=255)
                                                                                                  cTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaaa                                               <  2:2966174/149‑1 (MQ=255)
                                                                                                  cTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaaa                                               >  2:3389063/1‑149 (MQ=255)
                                                                                                   tGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAAt                                              <  2:3545301/149‑1 (MQ=255)
                                                                                                   tGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAAt                                              >  2:2012734/1‑149 (MQ=255)
                                                                                                    gTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATg                                             >  1:4772326/1‑149 (MQ=255)
                                                                                                      ccTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCg                                           <  1:2153872/149‑1 (MQ=255)
                                                                                                        tGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGa                                         >  2:2777911/1‑149 (MQ=255)
                                                                                                         gCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGaa                                        <  1:1650329/149‑1 (MQ=255)
                                                                                                          ccAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGaaa                                       >  1:4713463/1‑149 (MQ=255)
                                                                                                             gTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGc                                    <  1:2048885/149‑1 (MQ=255)
                                                                                                             gTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGc                                    <  2:4179088/149‑1 (MQ=255)
                                                                                                                cTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTcc                                 <  1:635428/149‑1 (MQ=255)
                                                                                                                   cATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagc                              >  2:4097876/1‑149 (MQ=255)
                                                                                                                     tCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagcga                            <  1:4141746/149‑1 (MQ=255)
                                                                                                                      ccgccgCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagcgag                           <  2:2869308/149‑1 (MQ=255)
                                                                                                                      ccgccgCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagcgag                           >  2:1038288/1‑149 (MQ=255)
                                                                                                                           gccTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGACATGCGGAAATGCTCCAGCGAGCGTTg                      <  2:254673/149‑1 (MQ=255)
                                                                                                                             cTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAg                    >  2:2196383/1‑149 (MQ=255)
                                                                                                                               gACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGtt                  >  2:3175632/1‑149 (MQ=255)
                                                                                                                                    aTGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTgc             <  2:3307596/149‑1 (MQ=255)
                                                                                                                                    aTGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTgc             >  2:2452014/1‑149 (MQ=255)
                                                                                                                                    aTGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTgc             <  2:4050765/149‑1 (MQ=255)
                                                                                                                                         gctcgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGcccc        <  1:1750305/149‑1 (MQ=255)
                                                                                                                                          ctcgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCca       >  1:2728461/1‑149 (MQ=255)
                                                                                                                                            cgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCcaca     >  2:3955472/1‑149 (MQ=255)
                                                                                                                                            cgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCcaca     <  2:1678178/149‑1 (MQ=255)
                                                                                                                                               tttCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCt  >  1:1985911/1‑149 (MQ=255)
                                                                                                                                               tttATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCt  <  1:4634070/149‑1 (MQ=255)
                                                                                                                                               tttATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCt  <  2:3710955/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCT  >  CP009273/1680496‑1680787

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: