Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,039,703 C→A 100% 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,7030CA99.0% 626.5 / ‑1.9 193A900S (GCG→TCG) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (94/97);  total (94/99)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 193 total aligned reads displayed.

CAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTA  >  CP009273/3039558‑3039842
                                                                                                                                                 |                                                                                                                                           
    aCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAc                                                                                                                                       >  2:499996/1‑148 (MQ=255)
     cAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAcc                                                                                                                                      >  2:2773934/1‑148 (MQ=255)
     cAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAcc                                                                                                                                      <  2:2533244/148‑1 (MQ=255)
     cAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAcc                                                                                                                                      >  1:71953/1‑148 (MQ=255)
       ggTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAg                                                                                                                                    >  1:3365428/1‑148 (MQ=255)
          tACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCgg                                                                                                                                 >  2:2469500/1‑148 (MQ=255)
           aCGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCggg                                                                                                                                <  2:2023540/148‑1 (MQ=255)
            cGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGtt                                                                                                                                          >  1:961437/1‑137 (MQ=255)
             ggTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGtt                                                                                                                              <  2:3130227/148‑1 (MQ=255)
              gTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTa                                                                                                                             <  1:677020/148‑1 (MQ=255)
                cGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATc                                                                                                                           <  2:46281/148‑1 (MQ=255)
                cGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATc                                                                                                                           >  1:1689343/1‑148 (MQ=255)
                  ccGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGACTACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAAAGGGCGGCGAGTTCACCAGCGGGTGATCtt                                                                                                                         >  2:1168472/1‑148 (MQ=255)
                  ccGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCtt                                                                                                                         >  2:3481856/1‑148 (MQ=255)
                  ccGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCtt                                                                                                                         <  2:170926/148‑1 (MQ=255)
                  ccGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCtt                                                                                                                         <  1:335060/148‑1 (MQ=255)
                   cGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTc                                                                                                                        <  1:1499349/148‑1 (MQ=255)
                    gTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTcc                                                                                                                       <  2:1875665/148‑1 (MQ=255)
                    gTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTcc                                                                                                                       <  2:1667552/148‑1 (MQ=255)
                     tAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCa                                                                                                                      >  1:596255/1‑148 (MQ=255)
                      aaaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAg                                                                                                                     >  1:3191143/1‑148 (MQ=255)
                      aaaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAg                                                                                                                     >  1:249308/1‑148 (MQ=255)
                         catcatCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCgg                                                                                                                  >  2:2193647/1‑148 (MQ=255)
                         catcatCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCgg                                                                                                                  >  2:1454586/1‑148 (MQ=255)
                            catcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                               >  2:3127785/1‑148 (MQ=255)
                            catcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCa                                                                                                               >  2:120310/1‑148 (MQ=255)
                             atcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAg                                                                                                              <  1:1342861/148‑1 (MQ=255)
                              tcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGa                                                                                                             <  1:3059550/148‑1 (MQ=255)
                               ccAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGac                                                                                                            >  1:83728/1‑148 (MQ=255)
                                  gACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAcc                                                                                                         <  1:1624036/148‑1 (MQ=255)
                                    cGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCgg                                                                                                       >  2:3264370/1‑148 (MQ=255)
                                     gTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGc                                                                                                      <  2:2808416/148‑1 (MQ=255)
                                     gTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGc                                                                                                      <  2:1236095/148‑1 (MQ=255)
                                      tttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCt                                                                                                     <  1:644665/148‑1 (MQ=255)
                                      tttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCt                                                                                                     <  1:2958729/148‑1 (MQ=255)
                                      tttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCt                                                                                                     <  2:2366217/148‑1 (MQ=255)
                                      tttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCt                                                                                                     <  1:2410100/148‑1 (MQ=255)
                                        tCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCttt                                                                                                   <  1:657909/148‑1 (MQ=255)
                                         cACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTc                                                                                                  >  1:2063813/1‑148 (MQ=255)
                                            tGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                               >  1:2940257/1‑148 (MQ=255)
                                             gTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCt                                                                                              >  2:888186/1‑148 (MQ=255)
                                             gTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCt                                                                                              <  1:2539765/148‑1 (MQ=255)
                                               cGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTgg                                                                                            >  2:787713/1‑148 (MQ=255)
                                                  ccAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                         <  2:2093988/148‑1 (MQ=255)
                                                   cAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCaa                                                                                        <  1:3119938/148‑1 (MQ=255)
                                                    aGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAt                                                                                       >  2:2450905/1‑148 (MQ=255)
                                                      tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                     <  2:2849888/148‑1 (MQ=255)
                                                      tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                     <  1:1666719/148‑1 (MQ=255)
                                                      tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                     >  2:1278292/1‑148 (MQ=255)
                                                       aTTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTc                                                                                    >  2:585311/1‑148 (MQ=255)
                                                            tCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTgcgc                                                                               >  1:1069232/1‑148 (MQ=255)
                                                                ccACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAt                                                                           <  2:836731/148‑1 (MQ=255)
                                                                 cacaCCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATa                                                                          >  1:1850172/1‑148 (MQ=255)
                                                                  acacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAg                                                                         >  2:1076697/1‑148 (MQ=255)
                                                                   cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                        >  2:1158291/1‑148 (MQ=255)
                                                                    acCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGcc                                                                       >  1:865805/1‑148 (MQ=255)
                                                                    acCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGcc                                                                       >  1:3255508/1‑148 (MQ=255)
                                                                      cTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAg                                                                     <  2:2572724/148‑1 (MQ=255)
                                                                            ggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgc                                                               >  2:3510909/1‑148 (MQ=255)
                                                                             gAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgcg                                                              <  2:345096/148‑1 (MQ=255)
                                                                               aTACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTc                                                            >  1:2890554/1‑148 (MQ=255)
                                                                                 aCCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGa                                                          <  1:2469500/148‑1 (MQ=255)
                                                                                     cAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATaaa                                                      <  2:2903706/148‑1 (MQ=255)
                                                                                     cAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATaaa                                                      <  1:1860234/148‑1 (MQ=255)
                                                                                      aaCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAg                                                     <  2:1069220/148‑1 (MQ=255)
                                                                                       aCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGc                                                    >  1:1193125/1‑148 (MQ=255)
                                                                                       aCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGc                                                    >  2:2059321/1‑148 (MQ=255)
                                                                                        cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                   <  1:186348/148‑1 (MQ=255)
                                                                                        cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                   <  1:186945/148‑1 (MQ=255)
                                                                                        cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                   <  1:802017/148‑1 (MQ=255)
                                                                                        cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                   >  2:2453606/1‑148 (MQ=255)
                                                                                         ttCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGa                                                  >  2:1426428/1‑148 (MQ=255)
                                                                                            aCGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                               <  1:499996/148‑1 (MQ=255)
                                                                                            aCGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                               <  1:572968/148‑1 (MQ=255)
                                                                                              ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAg                                             <  1:2916465/148‑1 (MQ=255)
                                                                                                   tACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCTTCGATAAAGCGATCCAGTTCCa                                        >  1:1027552/1‑148 (MQ=255)
                                                                                                   tACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCa                                        <  1:2022526/148‑1 (MQ=255)
                                                                                                    aCGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCAc                                       >  1:2702138/1‑148 (MQ=255)
                                                                                                          gcgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCt                                 <  1:3507688/148‑1 (MQ=255)
                                                                                                          gcgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCt                                 <  1:2123203/148‑1 (MQ=255)
                                                                                                            gcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCttt                               <  1:3169828/148‑1 (MQ=255)
                                                                                                            gcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCttt                               <  1:1493662/148‑1 (MQ=255)
                                                                                                                   cGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTc                        >  2:2331655/1‑148 (MQ=255)
                                                                                                                   cGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTc                        >  2:1732969/1‑148 (MQ=255)
                                                                                                                    ggCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCa                       <  1:2153506/148‑1 (MQ=255)
                                                                                                                        aCCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCg                   >  2:2568776/1‑148 (MQ=255)
                                                                                                                         ccAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCgg                  <  2:71953/148‑1 (MQ=255)
                                                                                                                            gTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTc               >  1:438429/1‑148 (MQ=255)
                                                                                                                             ttCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCa              <  2:2839378/148‑1 (MQ=255)
                                                                                                                             ttCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCa              >  1:3403831/1‑148 (MQ=255)
                                                                                                                             ttCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCa              <  1:1063986/148‑1 (MQ=255)
                                                                                                                               cGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAAc            <  1:41832/148‑1 (MQ=255)
                                                                                                                                gCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCACTTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAAcc           <  1:1816097/148‑1 (MQ=255)
                                                                                                                                 ctctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCa          >  2:3469974/1‑148 (MQ=255)
                                                                                                                                  tctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCAt         <  2:1946494/148‑1 (MQ=255)
                                                                                                                                  tctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCAt         <  1:1454586/148‑1 (MQ=255)
                                                                                                                                   ctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATc        <  2:2064460/148‑1 (MQ=255)
                                                                                                                                      aaTGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCTACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGc     <  2:2043865/148‑1 (MQ=255)
                                                                                                                                       aTGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGACAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCg    <  1:3341787/148‑1 (MQ=255)
                                                                                                                                         gtgCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTa  <  2:3365428/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTA  >  CP009273/3039558‑3039842

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: