Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,041,026 C→A 26.5% D459Y (GAC→TAC)  gcvP ← glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,041,0260CA26.5% 681.4 / 219.7 302D459Y (GAC→TAC) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (118/104);  new base A (49/31);  total (167/135)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATTT  >  CP009273/3040878‑3041175
                                                                                                                                                    |                                                                                                                                                     
gggATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCg                                                                                                                                                      <  2:210789/150‑1 (MQ=255)
gggATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCg                                                                                                                                                      >  2:482733/1‑150 (MQ=255)
  gATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTg                                                                                                                                                    >  2:4716716/1‑150 (MQ=255)
  gATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTg                                                                                                                                                    >  2:2315242/1‑150 (MQ=255)
      atcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGcc                                                                                                                                                <  2:3591130/150‑1 (MQ=255)
         gCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACg                                                                                                                                             <  2:3131342/150‑1 (MQ=255)
         gCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACg                                                                                                                                             >  2:1661340/1‑150 (MQ=255)
          ccTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGt                                                                                                                                            <  1:1249723/150‑1 (MQ=255)
            tGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCt                                                                                                                                          >  1:4783524/1‑150 (MQ=255)
                tCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGt                                                                                                                                      >  1:3072105/1‑150 (MQ=255)
                 cAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTc                                                                                                                                     >  2:4220257/1‑150 (MQ=255)
                    cgcCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAg                                                                                                                                  >  1:984014/1‑150 (MQ=255)
                    cgcCAGATCTTTACGCTCCAGCGAGTGCAAATAGCGCATCATATCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAg                                                                                                                                  >  1:2389032/1‑150 (MQ=255)
                     gcCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGc                                                                                                                                 >  2:677372/1‑150 (MQ=255)
                      ccAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCg                                                                                                                                >  1:2964160/1‑150 (MQ=255)
                       cAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgt                                                                                                                               >  2:900314/1‑150 (MQ=255)
                       cAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgt                                                                                                                               >  1:1550471/1‑150 (MQ=255)
                       cAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgt                                                                                                                               >  1:437210/1‑150 (MQ=255)
                        aGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCgtg                                                                                                                              >  2:3492265/1‑150 (MQ=255)
                          aTCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTc                                                                                                                            >  1:3995265/1‑150 (MQ=255)
                          aTCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTc                                                                                                                            <  2:249642/150‑1 (MQ=255)
                           tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                           <  1:2449495/150‑1 (MQ=255)
                           tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                           <  1:325507/150‑1 (MQ=255)
                            cTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGa                                                                                                                          <  2:2465981/150‑1 (MQ=255)
                            cTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGa                                                                                                                          >  2:4361637/1‑150 (MQ=255)
                             tttACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGAt                                                                                                                         <  2:4661434/150‑1 (MQ=255)
                               tACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGt                                                                                                                       >  1:655535/1‑150 (MQ=255)
                                 cGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTcc                                                                                                                     >  1:3776447/1‑150 (MQ=255)
                                  gCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCa                                                                                                                    <  1:1379067/150‑1 (MQ=255)
                                  gCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCACAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCa                                                                                                                    >  2:1600669/1‑150 (MQ=255)
                                      cAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGcc                                                                                                                <  1:3901292/150‑1 (MQ=255)
                                          gAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTgg                                                                                                            >  1:856344/1‑150 (MQ=255)
                                             tGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTa                                                                                                         >  1:2099759/1‑150 (MQ=255)
                                             tGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTa                                                                                                         >  1:4405660/1‑150 (MQ=255)
                                                atatAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCg                                                                                                      >  1:2525962/1‑150 (MQ=255)
                                                  ataGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGcc                                                                                                    >  1:2791075/1‑150 (MQ=255)
                                                  ataGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGcc                                                                                                    >  1:4602282/1‑150 (MQ=255)
                                                   taGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGccc                                                                                                   >  2:1260718/1‑150 (MQ=255)
                                                    aGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCca                                                                                                  <  2:1805696/150‑1 (MQ=255)
                                                      cgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagc                                                                                                >  2:1066561/1‑150 (MQ=255)
                                                        catcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagcag                                                                                              <  2:1928296/150‑1 (MQ=255)
                                                           catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                           <  2:4166862/150‑1 (MQ=255)
                                                            atTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACg                                                                                          >  1:2340663/1‑150 (MQ=255)
                                                              ttCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGtt                                                                                        >  2:3528087/1‑150 (MQ=255)
                                                              ttCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGtt                                                                                        <  2:1348160/150‑1 (MQ=255)
                                                               tCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTg                                                                                       <  2:463226/150‑1 (MQ=255)
                                                                 ggTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaa                                                                                     >  2:181927/1‑150 (MQ=255)
                                                                    ttCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAg                                                                                  >  2:4326212/1‑150 (MQ=255)
                                                                    ttCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAg                                                                                  <  2:519689/150‑1 (MQ=255)
                                                                     tCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGc                                                                                 >  1:781377/1‑150 (MQ=255)
                                                                       gCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTg                                                                               <  2:1130682/150‑1 (MQ=255)
                                                                         tgtgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCa                                                                             >  1:4097404/1‑150 (MQ=255)
                                                                                 cGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGttt                                                                     <  2:4518156/150‑1 (MQ=255)
                                                                                 cGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGttt                                                                     <  2:3376652/150‑1 (MQ=255)
                                                                                    ttAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCa                                                                  >  1:2189370/1‑150 (MQ=255)
                                                                                     tAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAc                                                                 <  1:2268018/150‑1 (MQ=255)
                                                                                      aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                                >  1:224122/1‑150 (MQ=255)
                                                                                         cacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGt                                                             >  1:2716595/1‑150 (MQ=255)
                                                                                          acCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTg                                                            >  1:3415946/1‑150 (MQ=255)
                                                                                            cGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgt                                                          <  2:2693390/150‑1 (MQ=255)
                                                                                            cGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgt                                                          >  1:2133839/1‑150 (MQ=255)
                                                                                              gATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttg                                                        <  2:1638704/150‑1 (MQ=255)
                                                                                               aTGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgt                                                       <  2:706042/150‑1 (MQ=255)
                                                                                                tGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgtt                                                      <  1:2172695/150‑1 (MQ=255)
                                                                                                  ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                    <  1:1769755/150‑1 (MQ=255)
                                                                                                  ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                    <  1:1748210/150‑1 (MQ=255)
                                                                                                   gTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtca                                                   <  1:702991/150‑1 (MQ=255)
                                                                                                    tGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtcat                                                  >  1:2321111/1‑150 (MQ=255)
                                                                                                         aTTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAgg                                             >  2:1952816/1‑150 (MQ=255)
                                                                                                         aTTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAgg                                             >  1:834804/1‑150 (MQ=255)
                                                                                                              gtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGAt                                        <  2:2054719/150‑1 (MQ=255)
                                                                                                                  tcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCa                                    <  1:319694/150‑1 (MQ=255)
                                                                                                                   cgcgCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCaa                                   >  2:3225023/1‑150 (MQ=255)
                                                                                                                   cgcgCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCaa                                   >  2:4260182/1‑150 (MQ=255)
                                                                                                                   cgcgCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCaa                                   >  2:1760954/1‑150 (MQ=255)
                                                                                                                   cgcgCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCaa                                   >  2:1760848/1‑150 (MQ=255)
                                                                                                                     cgcAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAcc                                 <  1:1945260/150‑1 (MQ=255)
                                                                                                                     cgcAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAcc                                 >  2:145050/1‑150 (MQ=255)
                                                                                                                     cgcAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAcc                                 >  2:145347/1‑150 (MQ=255)
                                                                                                                     cgcAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAcc                                 <  1:3492265/150‑1 (MQ=255)
                                                                                                                           aTCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTc                           >  2:530708/1‑150 (MQ=255)
                                                                                                                            tCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCa                          >  2:4409944/1‑150 (MQ=255)
                                                                                                                              gCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGa                        <  2:4695769/150‑1 (MQ=255)
                                                                                                                                aGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAat                      >  2:2989718/1‑150 (MQ=255)
                                                                                                                                aGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAat                      <  1:482831/150‑1 (MQ=255)
                                                                                                                                  gCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAatat                    >  2:1278741/1‑150 (MQ=255)
                                                                                                                                      gATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCt                <  1:230156/150‑1 (MQ=255)
                                                                                                                                       aTAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCACCCGCGTTCAGAATATCGCTa               >  2:4078018/1‑150 (MQ=255)
                                                                                                                                       aTAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTa               >  2:2871236/1‑150 (MQ=255)
                                                                                                                                         agagCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg             >  1:1104414/1‑150 (MQ=255)
                                                                                                                                         agagCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg             >  1:4685777/1‑150 (MQ=255)
                                                                                                                                         agagCGGCTGTCGTGAGCCACCTCTTTGTACAGCGTGTCGATGTACAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAATCTGCAATACGTTTTCACGCGTGGTTGTTTAAAAAAGGGTGATCCCAACCGCGTTCAGAATATAGCTACg             >  1:3633614/1‑150 (MQ=255)
                                                                                                                                           agCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCa           >  2:1415357/1‑150 (MQ=255)
                                                                                                                                               gCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGtt       >  1:1007154/1‑150 (MQ=255)
                                                                                                                                               gCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGtt       <  2:3290352/150‑1 (MQ=255)
                                                                                                                                                cTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTg      <  2:2743417/150‑1 (MQ=255)
                                                                                                                                                cTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTg      >  1:41869/1‑150 (MQ=255)
                                                                                                                                                cTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTg      >  1:43604/1‑150 (MQ=255)
                                                                                                                                                cTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTg      <  2:2286333/150‑1 (MQ=255)
                                                                                                                                                    cGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGAttt  >  1:4428252/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATTT  >  CP009273/3040878‑3041175

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: