Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,041,321 G→T 5.2% I360I (ATC→ATA 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,3210GT 5.2% 1251.6 / 26.3 324I360I (ATC→ATAgcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base G (148/158);  new base T (8/10);  total (156/168)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

ATTTCAGCCGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATT  >  CP009273/3041172‑3041470
                                                                                                                                                     |                                                                                                                                                     
ttttCAGCCGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCg                                                                                                                                                       <  2:736897/149‑1 (MQ=255)
 tttCAGCCGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGa                                                                                                                                                      >  1:3392834/1‑150 (MQ=255)
 tttCAGCCGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGa                                                                                                                                                      >  1:2045155/1‑150 (MQ=255)
    cAGCCGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATAc                                                                                                                                                   >  2:3390310/1‑150 (MQ=255)
      gCCGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGt                                                                                                                                                 >  1:603878/1‑150 (MQ=255)
        cGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCTATACGttt                                                                                                                                               <  1:233359/150‑1 (MQ=255)
        cGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGttt                                                                                                                                               >  2:2652018/1‑150 (MQ=255)
        cGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGttt                                                                                                                                               >  1:2746208/1‑150 (MQ=255)
           ttCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAg                                                                                                                                            <  2:5075/150‑1 (MQ=255)
            tCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAgg                                                                                                                                           <  2:1605201/150‑1 (MQ=255)
             cGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGc                                                                                                                                          >  1:3820663/1‑150 (MQ=255)
              ggCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATTGTGCGCATGGCGCAGTTTCAGACCTGTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGcc                                                                                                                                         >  2:3521746/1‑150 (MQ=255)
              ggCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGTGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGcc                                                                                                                                         <  2:3241380/150‑1 (MQ=255)
              ggCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGcc                                                                                                                                         >  1:415260/1‑150 (MQ=255)
                cACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCaa                                                                                                                                       <  1:3905591/150‑1 (MQ=255)
                cACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCATTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCaa                                                                                                                                       >  1:875750/1‑150 (MQ=255)
                 aCGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGGTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAAc                                                                                                                                      <  1:1577564/150‑1 (MQ=255)
                   gcgTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCg                                                                                                                                    <  1:954628/150‑1 (MQ=255)
                     gTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCggg                                                                                                                                  >  1:1475852/1‑150 (MQ=255)
                      tCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGc                                                                                                                                 >  2:2740026/1‑150 (MQ=255)
                        aGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCg                                                                                                                               >  2:1684719/1‑150 (MQ=255)
                          tACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTg                                                                                                                             >  1:1968022/1‑150 (MQ=255)
                              cccGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGACATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATaa                                                                                                                         >  1:1231975/1‑150 (MQ=255)
                                cGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAAc                                                                                                                       >  1:4068008/1‑150 (MQ=255)
                                  cTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACgg                                                                                                                     >  2:392288/1‑150 (MQ=255)
                                       tCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCATTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATAc                                                                                                                >  1:1560939/1‑150 (MQ=255)
                                        cGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACa                                                                                                               <  1:2888379/150‑1 (MQ=255)
                                         ggCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAg                                                                                                              <  2:1900459/150‑1 (MQ=255)
                                          gCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAgg                                                                                                             >  2:2238508/1‑150 (MQ=255)
                                           ccACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGc                                                                                                            >  2:1010511/1‑150 (MQ=255)
                                           ccACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGc                                                                                                            >  2:688765/1‑150 (MQ=255)
                                            cACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCTATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCt                                                                                                           >  1:4308210/1‑150 (MQ=255)
                                              cTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTgg                                                                                                         <  2:2302879/150‑1 (MQ=255)
                                              cTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTgg                                                                                                         <  1:3516524/150‑1 (MQ=255)
                                               ttCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGc                                                                                                        >  1:2096631/1‑150 (MQ=255)
                                               ttCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGc                                                                                                        <  2:2301844/150‑1 (MQ=255)
                                                tCCACACACAAGTTGTCGAAATAGTGCTCATGGCGCAGTTTCAGACATTTTTTTTGCAGTCCCGACTCCAGGAAATCGGTCATACGGTTAATGAGGTTAGAGATACGTTTCAGTCCAACCGAGCCGTGATAAACGGCATACAGGCTGGCg                                                                                                       >  2:836428/1‑150 (MQ=255)
                                                tCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCg                                                                                                       <  2:2076955/150‑1 (MQ=255)
                                                     acacaAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGtt                                                                                                  <  1:3266022/150‑1 (MQ=255)
                                                      cacaAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGttt                                                                                                 <  2:2868212/150‑1 (MQ=255)
                                                         aaGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGNTAAACGGCATACAGGCTGGCGATGTTTGcc                                                                                              <  1:4277824/150‑1 (MQ=255)
                                                          aGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCca                                                                                             <  2:3141801/150‑1 (MQ=255)
                                                           ggTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCcag                                                                                            >  1:2538526/1‑150 (MQ=255)
                                                            gtgtCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCcagc                                                                                           <  2:2599780/150‑1 (MQ=255)
                                                               tCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGt                                                                                        <  2:1580687/150‑1 (MQ=255)
                                                                cGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTa                                                                                       <  2:34341/150‑1 (MQ=255)
                                                                cGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTa                                                                                       <  2:3608823/150‑1 (MQ=255)
                                                                 gAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTAc                                                                                      >  2:4189002/1‑150 (MQ=255)
                                                                      aGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTggg                                                                                 >  2:2804074/1‑150 (MQ=255)
                                                                      aGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTggg                                                                                 <  2:4089963/150‑1 (MQ=255)
                                                                          cgcATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGt                                                                             <  1:2331284/150‑1 (MQ=255)
                                                                          cgcATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGt                                                                             <  2:1333719/150‑1 (MQ=255)
                                                                          cgcATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGt                                                                             >  1:226742/1‑150 (MQ=255)
                                                                               ggCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACaaa                                                                        <  2:1832936/150‑1 (MQ=255)
                                                                                 cgcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATg                                                                      <  2:1117958/150‑1 (MQ=255)
                                                                                 cgcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATg                                                                      >  1:3419597/1‑150 (MQ=255)
                                                                                     gTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTgg                                                                  <  2:2913579/150‑1 (MQ=255)
                                                                                     gTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTgg                                                                  >  1:493074/1‑150 (MQ=255)
                                                                                      tttCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGa                                                                 >  1:666632/1‑150 (MQ=255)
                                                                                       ttCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAg                                                                <  1:866265/150‑1 (MQ=255)
                                                                                        tCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGt                                                               <  1:2516509/150‑1 (MQ=255)
                                                                                         cAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGtt                                                              <  2:1393560/150‑1 (MQ=255)
                                                                                         cAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGtt                                                              <  2:1787238/150‑1 (MQ=255)
                                                                                          aGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTc                                                             >  2:3674382/1‑150 (MQ=255)
                                                                                          aGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTc                                                             <  1:4626644/150‑1 (MQ=255)
                                                                                          aGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTc                                                             >  2:687552/1‑150 (MQ=255)
                                                                                             ccTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCTATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                          >  1:1829836/1‑150 (MQ=255)
                                                                                              cTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                         >  1:3988812/1‑150 (MQ=255)
                                                                                              cTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                         >  2:822444/1‑150 (MQ=255)
                                                                                               tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCttt                                                        <  2:2536926/150‑1 (MQ=255)
                                                                                               tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCttt                                                        <  1:1114715/150‑1 (MQ=255)
                                                                                                 tttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTtct                                                      >  1:3810349/1‑150 (MQ=255)
                                                                                                    gttgCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAc                                                   <  2:439328/150‑1 (MQ=255)
                                                                                                     ttgCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcg                                                  >  2:1549355/1‑150 (MQ=255)
                                                                                                      tgCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcgg                                                 <  1:2839847/150‑1 (MQ=255)
                                                                                                          ggCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGa                                             <  1:2869352/150‑1 (MQ=255)
                                                                                                              cgccgcCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGACATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATg                                         >  2:4285603/1‑150 (MQ=255)
                                                                                                               gccgccAGGATATCGGTCAGACGGTGAATGCGGTTAGCTATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGt                                        >  2:1451290/1‑150 (MQ=255)
                                                                                                               gccgccAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGt                                        <  1:1011769/150‑1 (MQ=255)
                                                                                                                 cgccAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTg                                      <  1:2283942/150‑1 (MQ=255)
                                                                                                                  gccAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGc                                     >  1:4509704/1‑150 (MQ=255)
                                                                                                                   ccAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCt                                    <  1:1390850/150‑1 (MQ=255)
                                                                                                                   ccAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCt                                    >  1:3686756/1‑150 (MQ=255)
                                                                                                                     aGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTcg                                  >  2:4576390/1‑150 (MQ=255)
                                                                                                                     aGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTcg                                  <  1:1246484/150‑1 (MQ=255)
                                                                                                                     aGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTcg                                  >  2:2109970/1‑150 (MQ=255)
                                                                                                                        atatCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGa                               <  2:2036614/150‑1 (MQ=255)
                                                                                                                          atCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGt                             <  1:2211792/150‑1 (MQ=255)
                                                                                                                                cAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCAt                       >  1:3048879/1‑150 (MQ=255)
                                                                                                                                 aGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATc                      <  1:2494137/150‑1 (MQ=255)
                                                                                                                                 aGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATc                      >  1:2272558/1‑150 (MQ=255)
                                                                                                                                  gACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCg                     <  1:2138171/150‑1 (MQ=255)
                                                                                                                                       tGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATg                <  1:1980982/150‑1 (MQ=255)
                                                                                                                                       tGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATg                >  1:2592718/1‑150 (MQ=255)
                                                                                                                                            gCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAg           <  1:3727754/150‑1 (MQ=255)
                                                                                                                                            gCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAg           <  1:4324741/150‑1 (MQ=255)
                                                                                                                                                ttAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGcgcg       >  2:1325767/1‑150 (MQ=255)
                                                                                                                                                    cGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTAt   >  2:3839633/1‑150 (MQ=255)
                                                                                                                                                    cGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTAt   <  1:2659760/150‑1 (MQ=255)
                                                                                                                                                     gATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTAtt  >  2:881061/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTTCAGCCGCTTCGGCACGCGTCAGTACGCCCGCTTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATT  >  CP009273/3041172‑3041470

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: