Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,039,703 C→A 86.0% 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,7030CA86.0% 595.0 / 105.1 250A900S (GCG→TCG) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (20/15);  new base A (104/111);  total (124/126)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 251 total aligned reads displayed.

GAGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGC  >  CP009273/3039555‑3039847
                                                                                                                                                    |                                                                                                                                                
gAGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGgcgc                                                                                                                                                  <  2:3474465/149‑1 (MQ=255)
   cAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGtt                                                                                                                                               <  1:1395397/149‑1 (MQ=255)
      aaCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAc                                                                                                                                            <  1:3229377/149‑1 (MQ=255)
       aCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCAcc                                                                                                                                           <  1:3333946/149‑1 (MQ=255)
       aCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCAcc                                                                                                                                           <  1:945129/149‑1 (MQ=255)
        cAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCa                                                                                                                                          >  1:2988909/1‑149 (MQ=255)
        cAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCa                                                                                                                                          >  1:196050/1‑149 (MQ=255)
          ggTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGc                                                                                                                                        >  2:348870/1‑149 (MQ=255)
          ggTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGc                                                                                                                                        <  1:3688580/149‑1 (MQ=255)
            ttACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCgg                                                                                                                                      >  2:1287529/1‑149 (MQ=255)
             tACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCggg                                                                                                                                     >  1:3632579/1‑149 (MQ=255)
             tACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCggg                                                                                                                                     >  1:995304/1‑149 (MQ=255)
             tACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCggg                                                                                                                                     >  1:3179249/1‑149 (MQ=255)
              aCGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGt                                                                                                                                    <  1:2189242/149‑1 (MQ=255)
                 gTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTAt                                                                                                                                 <  1:1394353/149‑1 (MQ=255)
                  tCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATc                                                                                                                                >  2:1141700/1‑149 (MQ=255)
                   cGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCt                                                                                                                               >  1:2099508/1‑149 (MQ=255)
                     ccGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTACGGCGAGTTCACCAGCGGGTTATCTTc                                                                                                                             <  2:2228237/149‑1 (MQ=255)
                         aaaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGc                                                                                                                         <  1:2255008/149‑1 (MQ=255)
                         aaaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGc                                                                                                                         >  1:82066/1‑149 (MQ=255)
                          aaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCg                                                                                                                        >  1:1931934/1‑149 (MQ=255)
                             atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGcc                                                                                                                     >  2:1519207/1‑149 (MQ=255)
                              tcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                                    >  1:1760526/1‑149 (MQ=255)
                               catcCAGAGGTTTCAGTGTCGGCCAGTATTTTTCTGCCACACCTGCCTGGATTACCGCAACTGCACGGCTGTACGGATGCGCCCTCTCGGCGACCAGTTCGCTCTGAATGTGGGGCGATTTCTCCAGCGGGTTATCTTCCAGCGGcctg                                                                                                                   <  1:4579448/149‑3 (MQ=255)
                                atcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGa                                                                                                                  <  2:1741876/149‑1 (MQ=255)
                                atcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGa                                                                                                                  <  2:1741575/149‑1 (MQ=255)
                                   cAGACGTTTCACTGTCGTCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGacac                                                                                                               <  1:117738/149‑1 (MQ=255)
                                   cAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGacac                                                                                                               <  1:474596/149‑1 (MQ=255)
                                   cAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGacac                                                                                                               <  1:1115995/149‑1 (MQ=255)
                                   cAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGacac                                                                                                               <  2:3040898/149‑1 (MQ=255)
                                   cAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGacac                                                                                                               <  2:3157579/149‑1 (MQ=255)
                                       cGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGc                                                                                                           <  2:3750245/149‑1 (MQ=255)
                                       cGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGc                                                                                                           <  2:532770/149‑1 (MQ=255)
                                         tttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCtt                                                                                                         >  2:3763525/1‑149 (MQ=255)
                                             aCTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAc                                                                                                     <  2:3017488/149‑1 (MQ=255)
                                              cTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                                    <  2:1355112/149‑1 (MQ=255)
                                                    gCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                              <  2:1664006/149‑1 (MQ=255)
                                                      cAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAt                                                                                            <  1:1732969/149‑1 (MQ=255)
                                                        gTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                          <  1:184203/149‑1 (MQ=255)
                                                        gTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                          >  2:451044/1‑149 (MQ=255)
                                                         tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTc                                                                                         >  1:4625045/1‑149 (MQ=255)
                                                         tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTc                                                                                         >  1:896501/1‑149 (MQ=255)
                                                           tttGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTg                                                                                       <  1:2528444/149‑1 (MQ=255)
                                                           tttGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTg                                                                                       >  1:839281/1‑149 (MQ=255)
                                                                  gCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAt                                                                                >  2:3179831/1‑149 (MQ=255)
                                                                     acacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                             >  2:1965365/1‑149 (MQ=255)
                                                                      cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGcc                                                                            >  2:441074/1‑149 (MQ=255)
                                                                       acCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCa                                                                           <  2:4727154/149‑1 (MQ=255)
                                                                          tGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCa                                                                        >  1:2006196/1‑149 (MQ=255)
                                                                            ccGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATc                                                                      <  1:3783995/149‑1 (MQ=255)
                                                                             cGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcg                                                                     >  1:255631/1‑149 (MQ=255)
                                                                              gggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgc                                                                    <  1:4587047/149‑1 (MQ=255)
                                                                              gggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgc                                                                    <  1:2355794/149‑1 (MQ=255)
                                                                                gAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGt                                                                  <  2:490802/149‑1 (MQ=255)
                                                                                gAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGt                                                                  <  2:491090/149‑1 (MQ=255)
                                                                                 aaTACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTc                                                                 <  1:1659181/149‑1 (MQ=255)
                                                                                 aaTACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTc                                                                 <  1:276078/149‑1 (MQ=255)
                                                                                 aaTACCGCAACTTCACGGCTGTACGGAGGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTc                                                                 <  1:3020217/149‑1 (MQ=255)
                                                                                      cGCAACGTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATaa                                                            <  2:2353342/149‑1 (MQ=255)
                                                                                          aCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                        <  2:301787/149‑1 (MQ=255)
                                                                                          aCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                        <  1:2411425/149‑1 (MQ=255)
                                                                                          aCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                        >  2:3207927/1‑149 (MQ=255)
                                                                                          aCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                        >  1:3530512/1‑149 (MQ=255)
                                                                                           cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGATTGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGa                                                       <  1:4395068/149‑1 (MQ=255)
                                                                                           cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGa                                                       >  2:1804359/1‑149 (MQ=255)
                                                                                           cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGa                                                       >  2:226584/1‑149 (MQ=255)
                                                                                           cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGa                                                       >  2:4193617/1‑149 (MQ=255)
                                                                                             tCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATc                                                     >  1:265617/1‑149 (MQ=255)
                                                                                              cACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                                    <  2:2588064/149‑1 (MQ=255)
                                                                                                 ggCTGTACGGATGCGCCCCCTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGt                                                 <  1:1622724/149‑1 (MQ=255)
                                                                                                 ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGt                                                 <  1:2296109/149‑1 (MQ=255)
                                                                                                 ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGt                                                 <  2:4610/149‑1 (MQ=255)
                                                                                                  gCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGtt                                                <  1:4078906/149‑1 (MQ=255)
                                                                                                  gCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGtt                                                <  1:1688227/149‑1 (MQ=255)
                                                                                                     gTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCa                                             <  2:112215/149‑1 (MQ=255)
                                                                                                        cGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACtt                                          <  1:623602/149‑1 (MQ=255)
                                                                                                         ggATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACttt                                         <  1:4503768/149‑1 (MQ=255)
                                                                                                          gATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                        >  1:3834021/1‑149 (MQ=255)
                                                                                                          gATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                        >  1:127571/1‑149 (MQ=255)
                                                                                                          gATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                        >  2:1753120/1‑149 (MQ=255)
                                                                                                           aTGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGc                                       >  1:148244/1‑149 (MQ=255)
                                                                                                               gcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTc                                   >  2:3120118/1‑149 (MQ=255)
                                                                                                               gcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTc                                   >  2:4137276/1‑149 (MQ=255)
                                                                                                               gcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTc                                   <  2:902030/149‑1 (MQ=255)
                                                                                                                cccACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCa                                  <  1:3296047/149‑1 (MQ=255)
                                                                                                                   aCTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAt                               <  2:3341010/149‑1 (MQ=255)
                                                                                                                   aCTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAt                               >  2:275597/1‑149 (MQ=255)
                                                                                                                     tCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTc                             <  2:1746931/149‑1 (MQ=255)
                                                                                                                       ggCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAg                           >  1:3922166/1‑149 (MQ=255)
                                                                                                                        gCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGt                          <  1:74749/149‑1 (MQ=255)
                                                                                                                            ccAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGt                      >  1:4170169/1‑149 (MQ=255)
                                                                                                                              aGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTc                    >  2:3179073/1‑149 (MQ=255)
                                                                                                                                ttCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCaa                  >  2:1401625/1‑149 (MQ=255)
                                                                                                                                    ctctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCAt              <  2:699520/149‑1 (MQ=255)
                                                                                                                                     tctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATc             >  2:2764268/1‑149 (MQ=255)
                                                                                                                                     tctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATc             >  1:1049049/1‑149 (MQ=255)
                                                                                                                                     tctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCGTTCAGATTCAGTCGGTTCAACCATc             <  2:2526647/149‑1 (MQ=255)
                                                                                                                                       tGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAg           >  2:1988488/1‑149 (MQ=255)
                                                                                                                                            gtgCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTAc      <  2:3287360/149‑1 (MQ=255)
                                                                                                                                              gcggcgAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTAccc    <  2:240667/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGC  >  CP009273/3039555‑3039847

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: