Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,039,409 T→G 19.4% intergenic (‑149/+118) ygfF ← / ← gcvP putative NAD(P)‑dependent oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/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,4090TG19.4% 192.8 / 19.6 139intergenic (‑149/+118)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (31/78);  new base G (27/3);  total (58/81)
Fisher's exact test for biased strand distribution p-value = 9.31e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

TGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGA  >  CP009273/3039262‑3039556
                                                                                                                                                   |                                                                                                                                                    
tttttttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGt                                                                                                                                                      <  1:3384274/146‑1 (MQ=255)
tGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGt                                                                                                                                                      >  2:53763/1‑148 (MQ=255)
       ccGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTc                                                                                                                                               <  2:3486995/148‑1 (MQ=255)
          agaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                                                                            >  2:1887765/1‑148 (MQ=255)
          agaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCAcc                                                                                                                                            >  2:89311/1‑148 (MQ=255)
           gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAccc                                                                                                                                           <  2:3543843/148‑1 (MQ=255)
           gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGGGCCCGCAccc                                                                                                                                           >  1:1592478/1‑148 (MQ=255)
              aaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTaa                                                                                                                                        <  2:2990397/148‑1 (MQ=255)
                attaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAcc                                                                                                                                      <  1:3207261/148‑1 (MQ=255)
                 ttaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAccc                                                                                                                                     <  2:1182616/148‑1 (MQ=255)
                 ttaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAccc                                                                                                                                     <  2:2157665/148‑1 (MQ=255)
                  taattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCAAAcccc                                                                                                                                    >  2:345764/1‑147 (MQ=255)
                    attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCGCACCCAAAcccgca                                                                                                                                  >  2:1642200/1‑145 (MQ=255)
                    attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                  <  2:533075/148‑1 (MQ=255)
                     ttCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAAACCtctc                                                                                                                                 >  1:3236255/1‑148 (MQ=255)
                      tCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                                                                                <  2:344169/148‑1 (MQ=255)
                      tCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                                                                                <  2:1498284/148‑1 (MQ=255)
                         ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCTCACCCAAACCCTCTCCCCa                                                                                                                             >  2:3313900/1‑148 (MQ=255)
                          tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCGCACCCAAACCCAATACCCaa                                                                                                                            >  2:947961/1‑148 (MQ=255)
                           gACGTATCTGTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                           <  2:2428393/148‑1 (MQ=255)
                            aCGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                          <  2:2389426/148‑1 (MQ=255)
                             cGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCCCACCCCAAaccccaccccccacc                                                                                                                         >  1:1615876/1‑143 (MQ=255)
                             cGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGGGCCCGCACCCCAACACACACCCCcaaag                                                                                                                         >  1:3360602/1‑146 (MQ=255)
                               tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCTAACCCGCACCCCAAAgggg                                                                                                                       >  2:1885505/1‑148 (MQ=255)
                               tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGGGCCGGCCTGGTCCCCTCACCCAACCCCTCACCCCCAAtggg                                                                                                                       >  2:649386/1‑148 (MQ=255)
                                 tCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCCGTGCCGGCCTGGTGCCCGCACCCGAACCCACTCCCCCAAGGGGCg                                                                                                                     >  2:170230/1‑148 (MQ=255)
                                      gagaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgggg                                                                                                                <  1:548213/148‑1 (MQ=255)
                                       agagTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGTGAGGGTATATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                               <  2:2924861/148‑1 (MQ=255)
                                       agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGAAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                               <  1:960411/148‑1 (MQ=255)
                                       agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                               <  1:960951/148‑1 (MQ=255)
                                        gaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCCCCCGAAACCACACCCCACACGGGGGCgggggc                                                                                                              >  1:1776254/1‑146 (MQ=255)
                                          ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTg                                                                                                            <  1:2322533/148‑1 (MQ=255)
                                          ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTg                                                                                                            <  1:2322519/148‑1 (MQ=255)
                                           tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGAGCCAGTGCCGGCCTGGTCCCCTCACCCAAACACCAAACCCAAAGCGGCCAGGGGACTGt                                                                                                           >  1:2445024/1‑148 (MQ=255)
                                            gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTc                                                                                                          <  1:1055935/148‑1 (MQ=255)
                                                aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCAAACCCCCACCCCAAAGGGGCGGGGGGGCTggggggg                                                                                                      >  2:1539899/1‑142 (MQ=255)
                                                         cTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGt                                                                                            <  2:2787863/148‑1 (MQ=255)
                                                         cTCAATCTTCAATGTTGGCACGCTTATCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTGGCCCCCACCCCAACCCACAAACCCAAACGCCGAGGGGACGGTCTGTGCAccatttt                                                                                             >  1:2766301/1‑142 (MQ=255)
                                                          tCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTACCCACACCCGAAACCACTACCCAAAGGGGGGAGGGGCCGGGCTGCGCACGTTTGTa                                                                                           >  2:1355774/1‑148 (MQ=255)
                                                              tCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                       <  2:2920743/148‑1 (MQ=255)
                                                               cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCCCCCCAACCCTCTCCCCAAACGGGCGGGGGGCCAGTGGGTGCACTTTTGTACTTTg                                                                                      >  2:1576788/1‑148 (MQ=255)
                                                               cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTg                                                                                      <  1:2927519/148‑1 (MQ=255)
                                                                     tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTGACCCACACCCAAACCCACACCCCCAACGGGGGAGGGGGCCGGCGGAGCCCATGTGTGCATTGTCAtaa                                                                                >  1:2921192/1‑146 (MQ=255)
                                                                         ggCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACt                                                                            <  1:3548213/148‑1 (MQ=255)
                                                                           cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCGGTTCCCCCCCCCCCAAACCCCACCCCCACCGCGCGCGGGGGCCGGCAGCGCCCGTTTGTACATTGTCATCTGACaaa                                                                          >  2:1077027/1‑148 (MQ=255)
                                                                           cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCGCCTCACCCAAAACCACTACCCAAAGGGGCGAGGGGCCGGTCTGAGCACTTTTGTACGTTGTCATCTGACTaa                                                                          >  2:2851081/1‑148 (MQ=255)
                                                                               cTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg                                                                      <  1:2019463/148‑1 (MQ=255)
                                                                                  ttCCCACTCGCCCTTTCTGGGAGAGGGCGGGGTTGAGGGTAAATATTCGCGCCAGTGCCGGACTGTTCCCCTCAACCTAACCCTCTCCCCAAAGGGGCGAGGGTACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcg                                                                   <  1:55212/148‑1 (MQ=37)
                                                                                   tCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                  >  2:3424415/1‑148 (MQ=255)
                                                                                   tCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGGGGGGACTGTCGGAGCACTTGTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                  >  2:542404/1‑148 (MQ=255)
                                                                                      cccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCACTCCCCCAAGGGGCGAGGGGACTGGCTGAGCACGTTTGTGCGTTGTCCTCTGACGAAAAAGGCGCCGa                                                               >  2:3058980/1‑148 (MQ=255)
                                                                                        cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGGGCCCACCCCCCAACCCACACCCCCCACCGGCGCGGGGGCGGTCGGTGCCCGTTTTTTATTTGTGATCTGTATAAAAAGTCGccccag                                                             >  1:429823/1‑144 (MQ=255)
                                                                                           gCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCACACCCCAACCCACACCCCCAAGGGGCGGGGGGGCTGTCTGGGGACGTTTGTACGTTGTGCTATGACGAAAAAATCGCCGacccgc                                                          >  2:2387249/1‑143 (MQ=255)
                                                                                            cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGcc                                                         <  2:277930/148‑1 (MQ=255)
                                                                                            cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGcc                                                         <  1:1945383/148‑1 (MQ=255)
                                                                                            cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCCAAACCCCTCCCCCCAGGGGGCAGGGGCCTGTCAGAGCACGGTTTTTCTTTTTCAAATGACTAAAAAAGCGGCGCAGCGcc                                                         >  2:2778684/1‑148 (MQ=255)
                                                                                             ccTTTCCGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCt                                                        <  1:1539899/148‑1 (MQ=255)
                                                                                             ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCt                                                        <  1:2958552/148‑1 (MQ=255)
                                                                                              cTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt                                                       <  1:346430/148‑1 (MQ=255)
                                                                                               tttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCttt                                                      >  2:2972615/1‑148 (MQ=255)
                                                                                               ttacAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCttt                                                      <  2:1183430/145‑1 (MQ=255)
                                                                                                  cAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                   >  2:2241956/1‑148 (MQ=255)
                                                                                                  cAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCGAACCCACTACCCAAAGGGGCGCGGGGACTGTCTGAGCACTTTTGTACGTTGTCATATGACTAAAAATTCGCCGAAGCGCCGTTAGa                                                   >  2:2241978/1‑148 (MQ=255)
                                                                                                     ggAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                <  1:305352/148‑1 (MQ=255)
                                                                                                      gagaGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAAt                                               <  2:3373576/148‑1 (MQ=255)
                                                                                                        gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAg                                             <  1:2866661/148‑1 (MQ=255)
                                                                                                        gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAg                                             <  1:3018368/148‑1 (MQ=255)
                                                                                                        gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAg                                             <  2:1753748/148‑1 (MQ=255)
                                                                                                        agtggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAg                                             <  2:3510763/145‑1 (MQ=255)
                                                                                                           ggCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCg                                          <  1:3465643/148‑1 (MQ=255)
                                                                                                            gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGa                                         <  1:1622538/148‑1 (MQ=255)
                                                                                                            gCCGGGGTGAGGGTAAATATGCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGa                                         <  1:2100581/148‑1 (MQ=255)
                                                                                                             ccGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTACCCCTCTCCCCAAAGGGTCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                        <  1:2778684/148‑1 (MQ=255)
                                                                                                             ccGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                        <  2:3477560/148‑1 (MQ=255)
                                                                                                             ccGGGGTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                        <  2:2097446/148‑1 (MQ=255)
                                                                                                                gggTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCa                                     <  1:170230/148‑1 (MQ=255)
                                                                                                                  gTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                   <  2:1592485/148‑1 (MQ=255)
                                                                                                                  gTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                   <  2:1592478/148‑1 (MQ=255)
                                                                                                                    gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGa                                 <  1:1375674/148‑1 (MQ=255)
                                                                                                                    gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGa                                 <  1:1375642/148‑1 (MQ=255)
                                                                                                                     aGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGaa                                <  2:2003052/148‑1 (MQ=255)
                                                                                                                      gggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCACTCCCCCAAGGGGGGAGGGGGCGGTCTGAGGACGTTTGTTCATTGTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTTGAAAAAAGTCGAATAACTGaaa                               >  2:2201167/1‑147 (MQ=255)
                                                                                                                          aaaTATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                           <  2:1978023/148‑1 (MQ=255)
                                                                                                                            atatTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCAAACCCACTACCCCAACGGGCGACGGGGCTGTGTGTGCCCGTTTTTTAGTTTTCATATGACTAAAAAAGCGCCCCACCGCCTTTTGAAAAAAGTAGAATCACTGAATTAATgg                         >  2:2817262/1‑148 (MQ=255)
                                                                                                                             tatTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCACTACCCAAAGGGGCGAGGGGACTGTCTGATCACTTTTGTACTTTTTCATATGAATAAAAAGGCGCCGAAGCGCCTTTATAAAAAAATCGAATCAGTTAATTAATGGt                        >  2:3422822/1‑148 (MQ=37)
                                                                                                                                tCGCGCCAGTGCCGGCCTGGTGCCCGCACCCCAAACCACACCCCCCACGGGCGCGGGGAGTGGCTGTTGACGGTTTTACTTTGTGATATTACTAAAAAAGCGCCGCAGCGCCCTTAGAAAAAAATCGTATAAGATAATTAATGGTaat                     >  2:654675/1‑146 (MQ=255)
                                                                                                                                  gcgcCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCTAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCg                   <  1:1029791/148‑1 (MQ=255)
                                                                                                                                      cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTaa               <  2:2680487/148‑1 (MQ=255)
                                                                                                                                      cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTaa               <  1:1642200/148‑1 (MQ=255)
                                                                                                                                      cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTaa               <  1:1642091/148‑1 (MQ=255)
                                                                                                                                       aGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAAt              <  2:3197033/148‑1 (MQ=255)
                                                                                                                                       aGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAAt              <  1:3161219/148‑1 (MQ=255)
                                                                                                                                           ccGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGt          <  2:1092769/148‑1 (MQ=255)
                                                                                                                                               ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGc      <  2:2605032/148‑1 (MQ=255)
                                                                                                                                               ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGc      <  2:1012604/148‑1 (MQ=255)
                                                                                                                                                cTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCa     <  2:2445024/148‑1 (MQ=255)
                                                                                                                                                cTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCa     >  2:2281847/1‑148 (MQ=255)
                                                                                                                                                cTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCa     <  1:417799/148‑1 (MQ=255)
                                                                                                                                                  gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAgg   <  1:89311/148‑1 (MQ=255)
                                                                                                                                                   ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATAGGTACGCAGGa  >  1:678039/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGA  >  CP009273/3039262‑3039556

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: