Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,039,703 C→A 6.6% A900S (GCG→TCG)  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,039,7030CA 6.6% 984.6 / 29.6 257A900S (GCG→TCG) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (112/127);  new base A (10/7);  total (122/135)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 262 total aligned reads displayed.

AGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACC  >  CP009273/3039556‑3039851
                                                                                                                                                   |                                                                                                                                                    
   aGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCa                                                                                                                                                 >  2:110454/1‑150 (MQ=255)
    gAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCAc                                                                                                                                                <  2:2776925/150‑1 (MQ=255)
          gTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCgg                                                                                                                                          <  1:1570001/150‑1 (MQ=255)
           ttACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCggg                                                                                                                                         <  1:3765337/150‑1 (MQ=255)
            tACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGt                                                                                                                                        >  1:301215/1‑150 (MQ=255)
            tACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGATCTGAATGTGCGGCGCGTTCACCAGCGGGt                                                                                                                                        <  1:2028713/150‑1 (MQ=255)
             aCGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGtt                                                                                                                                       <  2:3423327/150‑1 (MQ=255)
             aCGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGtt                                                                                                                                       <  2:4190907/150‑1 (MQ=255)
              cGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTa                                                                                                                                      >  2:1998775/1‑150 (MQ=255)
                 tCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCt                                                                                                                                   <  1:3354773/150‑1 (MQ=255)
                  cGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCtt                                                                                                                                  >  1:4213262/1‑150 (MQ=255)
                    ccGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTcc                                                                                                                                <  1:2892321/150‑1 (MQ=255)
                     cGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCa                                                                                                                               <  1:4363721/150‑1 (MQ=255)
                      gTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAg                                                                                                                              >  2:2769505/1‑150 (MQ=255)
                      gTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAg                                                                                                                              <  1:673956/150‑1 (MQ=255)
                      gTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAg                                                                                                                              >  1:3842023/1‑150 (MQ=255)
                        aaaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCg                                                                                                                            <  1:816174/150‑1 (MQ=255)
                        aaaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCg                                                                                                                            <  1:3260438/150‑1 (MQ=255)
                           catcatCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGcc                                                                                                                         <  2:3957954/150‑1 (MQ=255)
                            atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                                        <  2:1698090/150‑1 (MQ=255)
                            atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                                        >  2:1975540/1‑150 (MQ=255)
                            atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                                        <  1:3783085/150‑1 (MQ=255)
                            atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                                        <  2:812863/150‑1 (MQ=255)
                            atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                                        >  2:4634447/1‑150 (MQ=255)
                              catcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGa                                                                                                                      >  2:3522859/1‑150 (MQ=255)
                               atcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGac                                                                                                                     <  1:744248/150‑1 (MQ=255)
                                tcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGaca                                                                                                                    <  2:4574547/150‑1 (MQ=255)
                                tcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGaca                                                                                                                    <  1:3352669/150‑1 (MQ=255)
                                  cAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAcc                                                                                                                  >  1:3660039/1‑150 (MQ=255)
                                  cAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGacaca                                                                                                                  <  1:2628679/150‑2 (MQ=255)
                                       gTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCtt                                                                                                             <  2:2506631/150‑1 (MQ=255)
                                          tCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCa                                                                                                          >  2:3379086/1‑150 (MQ=255)
                                           cACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAc                                                                                                         <  1:4401165/150‑1 (MQ=255)
                                            aCTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                                        <  1:1092162/150‑1 (MQ=255)
                                              tgtggtccAGTATTTGTCTGGCACTCGTGCCGGGAATACCGCAACTTCACGGCTGTTCGGATGCGCCGACTCGGCGACCTGTTCGCTCTGAATGTGCGGCGGGTTCACCAGGTGGTTATCTTCCAGCGGGCAGACTCCGGCTTTGAGCTg                                                                                                      <  2:3218389/144‑1 (MQ=255)
                                              tGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTg                                                                                                      <  1:3649633/150‑1 (MQ=255)
                                              tGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTg                                                                                                      >  1:1846185/1‑150 (MQ=255)
                                               gTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTgg                                                                                                     <  1:1539576/150‑1 (MQ=255)
                                                  ggCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                                  >  2:439049/1‑150 (MQ=255)
                                                   gCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCaa                                                                                                 >  2:3843531/1‑150 (MQ=255)
                                                     cAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAtt                                                                                               >  2:2285586/1‑150 (MQ=255)
                                                     cAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAtt                                                                                               >  2:1613522/1‑150 (MQ=255)
                                                     cAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAtt                                                                                               >  2:4099587/1‑150 (MQ=255)
                                                      aGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                              >  1:2812113/1‑150 (MQ=255)
                                                      aGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                              >  2:63496/1‑150 (MQ=255)
                                                        tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCt                                                                                            <  2:1804606/150‑1 (MQ=255)
                                                        tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCt                                                                                            <  2:1804432/150‑1 (MQ=255)
                                                        tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCt                                                                                            <  2:1905470/150‑1 (MQ=255)
                                                               cTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGa                                                                                     <  1:1352217/150‑1 (MQ=255)
                                                                tGCCACACCTGCCGGGAAGACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAt                                                                                    <  1:3474667/150‑1 (MQ=255)
                                                                    acacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGcc                                                                                >  2:2622949/1‑150 (MQ=255)
                                                                     cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCa                                                                               <  1:1178105/150‑1 (MQ=255)
                                                                       ccTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGc                                                                             <  1:4249600/150‑1 (MQ=255)
                                                                         tGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGgcgc                                                                                                                                                      <  1:1437361/75‑1 (MQ=255)
                                                                         tGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCAt                                                                           >  2:2981850/1‑150 (MQ=255)
                                                                         tGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCAt                                                                           >  2:2126907/1‑150 (MQ=255)
                                                                         tGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCAt                                                                           >  1:3134247/1‑150 (MQ=255)
                                                                            cGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgc                                                                        >  2:1303060/1‑150 (MQ=255)
                                                                             gggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgcg                                                                       >  2:4633640/1‑150 (MQ=255)
                                                                              ggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCGGCATCGCGt                                                                      >  1:968680/1‑150 (MQ=255)
                                                                               gAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTc                                                                     >  1:4562332/1‑150 (MQ=255)
                                                                                aaTACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCg                                                                    <  1:3612989/150‑1 (MQ=255)
                                                                                aaTACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCg                                                                    <  2:1436809/150‑1 (MQ=255)
                                                                                aaTACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCg                                                                    <  2:3167936/150‑1 (MQ=255)
                                                                                   aCCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATa                                                                 <  2:2124553/150‑1 (MQ=255)
                                                                                       cAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGc                                                             <  1:4087627/150‑1 (MQ=255)
                                                                                       cAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGc                                                             >  2:3385454/1‑150 (MQ=255)
                                                                                        aaCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCAATAAAGCg                                                            >  2:4482221/1‑150 (MQ=255)
                                                                                          cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGAt                                                          >  1:2061127/1‑150 (MQ=255)
                                                                                            tCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                                        <  2:1202160/150‑1 (MQ=255)
                                                                                            tCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                                        <  2:4253352/150‑1 (MQ=255)
                                                                                             cACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCa                                                       >  2:1538258/1‑150 (MQ=255)
                                                                                               cGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGt                                                     <  1:2793344/150‑1 (MQ=255)
                                                                                                ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGtt                                                    <  1:2444288/150‑1 (MQ=255)
                                                                                                    gTACGGCTGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCAc                                                >  1:1477051/1‑150 (MQ=255)
                                                                                                      aCGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACtt                                              <  1:99568/150‑1 (MQ=255)
                                                                                                      aCGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACtt                                              <  1:525004/150‑1 (MQ=255)
                                                                                                       cGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACttt                                             >  1:2359782/1‑150 (MQ=255)
                                                                                                       cGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACttt                                             <  1:560762/150‑1 (MQ=255)
                                                                                                       cGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACttt                                             >  2:4358807/1‑150 (MQ=255)
                                                                                                       cGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACttt                                             >  2:3655983/1‑150 (MQ=255)
                                                                                                        ggATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                            >  2:2483668/1‑150 (MQ=255)
                                                                                                        ggATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                            <  1:1034659/150‑1 (MQ=255)
                                                                                                        ggATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                            <  1:1034640/150‑1 (MQ=255)
                                                                                                            gcgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCttt                                        <  2:886164/150‑1 (MQ=255)
                                                                                                            gcgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCttt                                        <  1:3712615/150‑1 (MQ=255)
                                                                                                             cgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTc                                       >  2:3112694/1‑150 (MQ=255)
                                                                                                                       gCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTc                             <  1:3319452/150‑1 (MQ=255)
                                                                                                                       gCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTc                             <  1:2591641/150‑1 (MQ=255)
                                                                                                                       gCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTc                             <  1:2091648/150‑1 (MQ=255)
                                                                                                                          aCCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGt                          <  1:3297232/150‑1 (MQ=255)
                                                                                                                            cAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTc                        >  1:3889175/1‑150 (MQ=255)
                                                                                                                              gTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCaa                      >  2:1370539/1‑150 (MQ=255)
                                                                                                                                     ctGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAg               <  2:302993/150‑1 (MQ=255)
                                                                                                                                     ctGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAg               <  2:3561715/150‑1 (MQ=255)
                                                                                                                                      tGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGc              >  2:4070636/1‑150 (MQ=255)
                                                                                                                                       gAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCg             <  1:3158759/150‑1 (MQ=255)
                                                                                                                                              cggcgCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGc      >  1:3553586/1‑150 (MQ=255)
                                                                                                                                               ggcgCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGcc     >  1:934322/1‑150 (MQ=255)
                                                                                                                                                  gcgTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAcc  <  2:730272/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACC  >  CP009273/3039556‑3039851

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: