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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,041,0260CA 6.8% 1033.3 / 32.8 266D459Y (GAC→TAC) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (142/106);  new base A (10/8);  total (152/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 270 total aligned reads displayed.

CGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATT  >  CP009273/3040877‑3041174
                                                                                                                                                     |                                                                                                                                                    
cGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTc                                                                                                                                                      <  2:2627983/150‑1 (MQ=255)
 gggATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCg                                                                                                                                                     >  2:1635567/1‑150 (MQ=255)
   gATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTg                                                                                                                                                   >  2:1019611/1‑150 (MQ=255)
    atcatcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGa                                                                                                                                                  >  2:3092589/1‑150 (MQ=255)
    atcatcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGa                                                                                                                                                  <  2:2756646/150‑1 (MQ=255)
    atcatcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGa                                                                                                                                                  <  2:3258330/150‑1 (MQ=255)
      catcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGc                                                                                                                                                >  1:329918/1‑150 (MQ=255)
        tcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCa                                                                                                                                              <  2:1085490/150‑1 (MQ=255)
          gCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACg                                                                                                                                            >  2:3804659/1‑150 (MQ=255)
           ccTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGt                                                                                                                                           >  2:1893542/1‑150 (MQ=255)
            cTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTc                                                                                                                                          <  1:3748872/150‑1 (MQ=255)
             tGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCt                                                                                                                                         >  1:2197571/1‑150 (MQ=255)
             tGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCt                                                                                                                                         >  1:3127905/1‑150 (MQ=255)
              gATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCtt                                                                                                                                        >  2:3832928/1‑150 (MQ=255)
              gATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCtt                                                                                                                                        <  2:2085910/150‑1 (MQ=255)
                ttCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTg                                                                                                                                      >  2:2331294/1‑150 (MQ=255)
                ttCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTg                                                                                                                                      >  2:1731148/1‑150 (MQ=255)
                 tCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGt                                                                                                                                     <  2:44900/150‑1 (MQ=255)
                    gcgcCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCa                                                                                                                                  <  1:1268918/150‑1 (MQ=255)
                           aTCTTTACGCTCCAGAGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTc                                                                                                                           >  2:539244/1‑150 (MQ=255)
                            tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                          >  1:1099512/1‑150 (MQ=255)
                            tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                          <  2:2196924/150‑1 (MQ=255)
                            tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                          <  2:2123016/150‑1 (MQ=255)
                            tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                          <  1:515698/150‑1 (MQ=255)
                               ttACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATg                                                                                                                       >  2:766535/1‑150 (MQ=255)
                                 aCGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATTGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTc                                                                                                                     >  2:2892870/1‑150 (MQ=255)
                                 aCGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTc                                                                                                                     >  2:4066009/1‑150 (MQ=255)
                                  cGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTcc                                                                                                                    <  1:2484072/150‑1 (MQ=255)
                                   gCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCa                                                                                                                   >  2:143080/1‑150 (MQ=255)
                                   gCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCa                                                                                                                   <  1:333052/150‑1 (MQ=255)
                                    cTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAg                                                                                                                  >  1:343083/1‑150 (MQ=255)
                                        aGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCg                                                                                                              >  2:3750613/1‑150 (MQ=255)
                                        aGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCg                                                                                                              <  1:2820428/150‑1 (MQ=255)
                                        aGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCg                                                                                                              <  2:1507963/150‑1 (MQ=255)
                                           gAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTgg                                                                                                           <  2:2966673/150‑1 (MQ=255)
                                            aGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGt                                                                                                          >  1:1544206/1‑150 (MQ=255)
                                              tGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTa                                                                                                        >  1:2768767/1‑150 (MQ=255)
                                               gCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTAt                                                                                                       >  1:3554795/1‑150 (MQ=255)
                                                 atatAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCg                                                                                                     >  1:2953575/1‑150 (MQ=255)
                                                  tataGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGc                                                                                                    >  1:1536670/1‑150 (MQ=255)
                                                  tataGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGc                                                                                                    >  2:263553/1‑150 (MQ=255)
                                                       cgcATCATTTCGGTTTCGCTTTGGTAGCGATTAAAAACCGGATGGGGGAGGATATCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGATGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagc                                                                                               >  1:854628/1‑150 (MQ=255)
                                                        gcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagca                                                                                              <  1:2610985/150‑1 (MQ=255)
                                                         catcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagcag                                                                                             >  1:1130493/1‑150 (MQ=255)
                                                          atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                            >  2:3685939/1‑150 (MQ=255)
                                                            catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                          <  2:3990076/150‑1 (MQ=255)
                                                            catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                          >  1:1402541/1‑150 (MQ=255)
                                                            catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                          >  2:3237517/1‑150 (MQ=255)
                                                              tttCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGt                                                                                        >  1:1622757/1‑150 (MQ=255)
                                                               ttCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGtt                                                                                       <  2:3445667/150‑1 (MQ=255)
                                                                  ggTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaa                                                                                    <  2:2235001/150‑1 (MQ=255)
                                                                  ggTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaa                                                                                    >  2:3309975/1‑150 (MQ=255)
                                                                    tttCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaaaa                                                                                  >  2:2991394/1‑150 (MQ=255)
                                                                     ttCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAg                                                                                 >  1:624937/1‑150 (MQ=255)
                                                                       cGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCt                                                                               <  2:3525202/150‑1 (MQ=255)
                                                                            tgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCAtt                                                                          <  1:1965085/150‑1 (MQ=255)
                                                                              gTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTAc                                                                        <  1:2560404/150‑1 (MQ=255)
                                                                                 gctaTTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGtt                                                                     >  2:2074960/4‑150 (MQ=255)
                                                                                  cGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGttt                                                                    >  1:2861098/1‑150 (MQ=255)
                                                                                   gATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGGTGAAGAGCTGCATTACGttgt                                                                   >  2:2802056/1‑148 (MQ=255)
                                                                                    aTTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTc                                                                  >  1:2641027/1‑150 (MQ=255)
                                                                                     ttAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCa                                                                 <  1:3240499/150‑1 (MQ=255)
                                                                                      tAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAc                                                                >  1:3023952/1‑150 (MQ=255)
                                                                                       aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                               >  1:3064741/1‑150 (MQ=255)
                                                                                         acacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcgcg                                                             >  1:2413785/1‑150 (MQ=255)
                                                                                         acacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcgcg                                                             <  2:2540634/150‑1 (MQ=255)
                                                                                         acacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcgcg                                                             >  1:3722174/1‑150 (MQ=255)
                                                                                          cacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGt                                                            >  2:2315336/1‑150 (MQ=255)
                                                                                            ccGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTgg                                                          <  1:2273858/150‑1 (MQ=255)
                                                                                             cGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgt                                                         >  2:3571852/1‑150 (MQ=255)
                                                                                                aTGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgt                                                      >  2:1991420/1‑150 (MQ=255)
                                                                                                   ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                   <  2:1369641/150‑1 (MQ=255)
                                                                                                   ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                   <  2:2984367/150‑1 (MQ=255)
                                                                                                   ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                   >  2:3680405/1‑150 (MQ=255)
                                                                                                    gTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtca                                                  <  1:4168002/150‑1 (MQ=255)
                                                                                                        ggATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCaa                                              >  1:4144090/1‑150 (MQ=255)
                                                                                                         gATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAg                                             >  2:705359/1‑150 (MQ=255)
                                                                                                         gATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAg                                             <  2:809541/150‑1 (MQ=255)
                                                                                                          aTTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAgg                                            <  1:4051213/150‑1 (MQ=255)
                                                                                                             tcgtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTg                                         >  2:1021960/1‑150 (MQ=255)
                                                                                                             tcgtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTg                                         <  1:3605203/150‑1 (MQ=255)
                                                                                                                tcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATc                                      <  1:2428513/150‑1 (MQ=255)
                                                                                                                 cgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATcc                                     >  1:111018/1‑150 (MQ=255)
                                                                                                                 cgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATcc                                     <  2:2377652/150‑1 (MQ=255)
                                                                                                                 cgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATcc                                     <  1:202739/150‑1 (MQ=255)
                                                                                                                  gtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATccc                                    <  2:2441063/150‑1 (MQ=255)
                                                                                                                   tcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCa                                   >  1:1948943/1‑150 (MQ=255)
                                                                                                                    cgcgCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCaa                                  >  1:3993966/1‑150 (MQ=255)
                                                                                                                     gcgcAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAc                                 <  2:2750186/150‑1 (MQ=255)
                                                                                                                      cgcAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAcc                                >  1:632985/1‑150 (MQ=255)
                                                                                                                            aTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTc                          <  1:93518/150‑1 (MQ=255)
                                                                                                                                cAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGaa                      >  2:1569243/1‑150 (MQ=255)
                                                                                                                                 aGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAat                     >  1:2182056/1‑150 (MQ=255)
                                                                                                                                 aGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAat                     >  1:1384804/1‑150 (MQ=255)
                                                                                                                                  ggCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAata                    >  1:1712083/1‑150 (MQ=255)
                                                                                                                                    cTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATc                  <  1:849876/150‑1 (MQ=255)
                                                                                                                                      ggATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGc                >  2:732611/1‑150 (MQ=255)
                                                                                                                                        aTAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTa              >  2:478075/1‑150 (MQ=255)
                                                                                                                                          agagCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg            <  1:3504810/150‑1 (MQ=255)
                                                                                                                                                   gTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGAt   <  1:4089699/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATT  >  CP009273/3040877‑3041174

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: