Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,050 T→G 10.1% intergenic (‑149/+118) ygfF ← / ← gcvP putative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,046,0500TG10.1% 445.2 / 11.4 204intergenic (‑149/+118)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (74/105);  new base G (18/7);  total (92/112)
Fisher's exact test for biased strand distribution p-value = 4.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

TTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAG  >  NC_000913/3045902‑3046198
                                                                                                                                                    |                                                                                                                                                    
ttGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGtt                                                                                                                                                     >  2:3019884/1‑150 (MQ=255)
   tttttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTccc                                                                                                                                                  <  2:1147132/150‑1 (MQ=255)
     tttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTcccca                                                                                                                                                >  2:2858573/1‑149 (MQ=255)
     tttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                                                                                                >  2:2881842/1‑150 (MQ=255)
     tttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTAGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTccccg                                                                                                                                                >  2:2473378/1‑149 (MQ=255)
      ttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTc                                                                                                                                               <  2:1755566/150‑1 (MQ=255)
            gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTa                                                                                                                                         <  1:1369653/150‑1 (MQ=255)
            gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTa                                                                                                                                         <  1:3975619/150‑1 (MQ=255)
            gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCccccaa                                                                                                                                         >  1:3857065/1‑148 (MQ=255)
                   taattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                  >  2:2211343/1‑150 (MQ=255)
                        cGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACACCCCa                                                                                                                             >  2:1691509/1‑150 (MQ=255)
                         gTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCCAACCCAAAACCACACCCCac                                                                                                                            >  2:697668/1‑149 (MQ=255)
                              cGTAACGTTGAGAGTGCTCACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGACTGTTCCCCTCACCCTAACCCTCTCCCCAAAgggg                                                                                                                       <  2:2240544/150‑1 (MQ=37)
                                 aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCCCCCAAACCCACTCCCCAACGGGGCGa                                                                                                                    >  1:1928939/1‑150 (MQ=255)
                                   cTTTTAGATTGCTAACACACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCCGGGTGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                  <  1:1786108/150‑1 (MQ=255)
                                   cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                  <  2:2042543/150‑1 (MQ=255)
                                    tttGAGATTGCTAACAAACTGGCTCAATCTGCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTCAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCGCACCCGCTCCCCACACGGGCGCggg                                                                                                                 >  2:2794181/1‑150 (MQ=255)
                                      tGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                               >  2:2086768/1‑150 (MQ=255)
                                               tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                                                                      >  2:939215/1‑150 (MQ=255)
                                               tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAAACCACTCCCCAAAGGGGCGAGGGGGCGGTCTGAg                                                                                                      >  1:3527035/1‑150 (MQ=37)
                                               tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCACACCCCAAAGGGGCGCGGGGGCTGTCTGAg                                                                                                      >  1:1910173/1‑150 (MQ=255)
                                               tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCCCCCAAACCCTCACCCCACACGGGCGGGGGGACAGTCTGAg                                                                                                      >  1:878804/1‑150 (MQ=255)
                                                   aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACt                                                                                                  <  1:2104365/150‑1 (MQ=255)
                                                      cTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACtttt                                                                                               <  1:2328196/150‑1 (MQ=255)
                                                      cTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACtttt                                                                                               <  1:1003150/150‑1 (MQ=255)
                                                         gCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCGAAACCACACCCCAAAGGGGCGAGGGGACGGGCTTGTCACTTTTGTa                                                                                            >  2:1539448/1‑150 (MQ=25)
                                                           tCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACt                                                                                          <  1:338602/150‑1 (MQ=255)
                                                           tCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCTCACCCCAAACGGGCGAGGGGGCGGTCTGAGCACTTTTGTAcg                                                                                          >  2:865868/1‑149 (MQ=255)
                                                            cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTCCtt                                                                                         >  2:3679243/1‑150 (MQ=255)
                                                            cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGCGGGGACTGTCTGGGCACGTTTGGAcgt                                                                                         >  2:2506114/1‑148 (MQ=255)
                                                            cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACACCCCCAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTAcat                                                                                         >  1:3619192/1‑148 (MQ=255)
                                                            cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCAAACCCACTCCCCACAGGGGCGCGGGGGCAGTCTGAGCACATTTGTACtt                                                                                         >  2:424210/1‑150 (MQ=255)
                                                             aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAACATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                        >  2:1041744/1‑150 (MQ=255)
                                                             aaTCTTCAAAGTTGCCACGCTCTTCCCCCACCCCCTTACAGGGAGAGGGCCGCGGTGAGGGTAAATATTCGCGCCAGTCCCGGCCTGTTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                        <  2:1446/150‑1 (MQ=25)
                                                              aTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCCCACCCAAACCCACACCCCCCAGGGGCGCGGGGCCTGTCTGAGCACAGTTGTACTTTg                                                                                       >  1:817718/1‑150 (MQ=255)
                                                               tCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGt                                                                                      >  2:3330274/1‑150 (MQ=255)
                                                               tCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCGAACCCACACCCCCAAGGGGCGAGGGGCCGGTCTGAGCACATTTGTACttttt                                                                                      >  2:1896230/1‑148 (MQ=255)
                                                                 ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCa                                                                                    <  2:2801358/150‑1 (MQ=255)
                                                                 ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCCCTTTTGTACTTTGTCa                                                                                    >  2:3001648/1‑150 (MQ=255)
                                                                    actGTTGGCACGCTTTTCCCCCTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                                                                                 <  2:1785699/148‑1 (MQ=255)
                                                                    aaTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                                                                                 <  1:207217/150‑1 (MQ=255)
                                                                    aaTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                                                                                 <  1:2550171/150‑1 (MQ=255)
                                                                       gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCTCACCCAAACCCACTACCCCAAGGGGCGGGGGGGCTGTCTGAGCACATTTGTAAGTTGTGATCTGAc                                                                              >  1:2237704/1‑150 (MQ=255)
                                                                       gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                                                                              <  2:281799/150‑1 (MQ=255)
                                                                       gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                                                                              <  2:2923603/150‑1 (MQ=255)
                                                                        ttGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACt                                                                             <  2:2719709/150‑1 (MQ=255)
                                                                            cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCGAACCCACACCCCAAAGGGGCGCGGGGGCTGTCTGAGCCCGTTTGTACTTTGTCATATGACTaaaa                                                                         >  1:2167443/1‑150 (MQ=255)
                                                                             aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaaa                                                                        <  2:1478454/150‑1 (MQ=255)
                                                                             aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaaa                                                                        >  2:2691369/1‑150 (MQ=255)
                                                                             aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCAAACCCGCTACCCCAACGGGGGGGGGGCCGGTCTGAGCACGTTTGTACTTTGTGATATGAATaaaaa                                                                        >  1:3239927/1‑150 (MQ=255)
                                                                             aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCAAACCCACTACCCAAAGGGGCGCGGGGGCTGTCTGAGGACGTTTGTTCTTTGTAATCTGTCTaaaaa                                                                        >  2:879470/1‑150 (MQ=255)
                                                                              cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg                                                                       <  1:1642940/150‑1 (MQ=255)
                                                                              cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCCGACTAAAAAg                                                                       <  2:3887709/150‑1 (MQ=255)
                                                                                cTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgc                                                                     <  2:2463529/150‑1 (MQ=255)
                                                                                        ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc                                                             <  2:166753/150‑1 (MQ=255)
                                                                                        ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc                                                             <  2:2969130/150‑1 (MQ=255)
                                                                                          tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCGCACCCCAAACCACAACCCAAACGGGGGCGGGGGCGGTGTGAGCACTTTTGTACTTTTTAAAATGACAAAAAATTCGCCGAAgcgc                                                           >  1:746550/1‑150 (MQ=255)
                                                                                          tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCACACCCAAAACCCATCCCCCAAGGGGCGGGGGGGCTGTGTGGGCACTGTTGTTATTTGTCATATTAATAAAAAGGCGCCGAAgcgc                                                           >  1:3660616/1‑150 (MQ=255)
                                                                                          tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                           >  2:2699450/1‑150 (MQ=255)
                                                                                          tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCGAAACCCCACCCCCAAGGGGGGCGGGGGCAGGGTGTGCCCTGTTGTTCATTGTAATCGGACAAAAAAGGCGCCCAAgcgc                                                           >  1:1611323/1‑150 (MQ=255)
                                                                                          tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCCAACCCACTACCCCAAGGGGCGGGGGGGCGGTGTGAGCCCATTTGTACTTTTTGATATTAAGAAAAAGGCGcccccgcgc                                                           >  1:3243680/1‑143 (MQ=255)
                                                                                           cGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGcc                                                          <  1:1774189/150‑1 (MQ=255)
                                                                                           cGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGcc                                                          <  1:1188417/150‑1 (MQ=255)
                                                                                               cTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTa                                                      <  1:1382313/150‑1 (MQ=255)
                                                                                                tttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAg                                                     >  2:321591/1‑150 (MQ=255)
                                                                                                tttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAg                                                     <  2:529115/150‑1 (MQ=255)
                                                                                                tttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAg                                                     <  2:3922669/150‑1 (MQ=255)
                                                                                                tttAAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCCCTTTTGTACTGTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAg                                                     >  2:564250/1‑150 (MQ=255)
                                                                                                 ttCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCttttta                                                    <  1:1016050/150‑4 (MQ=25)
                                                                                                 ttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                    <  2:1973897/150‑1 (MQ=255)
                                                                                                 ttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                    <  2:2615043/150‑1 (MQ=255)
                                                                                                    cgggAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTGCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                 <  1:1189806/149‑1 (MQ=255)
                                                                                                    aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                 <  1:3865390/150‑1 (MQ=255)
                                                                                                      ggAGAGGGCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                               <  2:3505990/150‑1 (MQ=255)
                                                                                                      ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                               <  2:3747311/150‑1 (MQ=255)
                                                                                                         gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                            <  1:2103541/150‑1 (MQ=255)
                                                                                                         gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                            <  1:435421/150‑1 (MQ=255)
                                                                                                           gggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCAAACCCACACCCCCAAGGGGCGAGGGGGCAGTCTGAGCACTTTTTTACGTTTTGAAATGACTAAAAATGCGCCCAAGCCCCATTAGAAAATAGTcca                                          >  1:2613627/1‑148 (MQ=255)
                                                                                                             gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCTCACCCAAACCCACACCCCACAGGGGCGCGGGGGCGGTCTGGGCCCGTTTTTTCTTTGTCCTCAGTCGAAAAATGCCCCGCAGCGCCATTAGGAAAAAGTCGAAt                                        >  2:3110420/1‑150 (MQ=255)
                                                                                                                        ggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTTAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATTGTCGAATCAGGGAAtttc                             >  2:2431257/1‑148 (MQ=255)
                                                                                                                         gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                            >  2:646585/1‑150 (MQ=255)
                                                                                                                         gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                            <  2:544709/150‑1 (MQ=255)
                                                                                                                         gTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                            <  2:2185497/150‑1 (MQ=255)
                                                                                                                          tAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTg                           <  2:822715/150‑1 (MQ=255)
                                                                                                                          tAAACATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTg                           <  1:1041744/150‑1 (MQ=255)
                                                                                                                              tatTCGCGCCAGTGCCGGCCTGTTGCCCACACCCAAAACCACACCCCAAAGCGGCGAGGGGCCTGTCTGTGCACATTTGTACTTTTTGAACTGTCTAAAAAGGCGCCGAAGCGCCGATAGAAAAAAGTAGAAAAAATGAATTAAAGGTAt                       >  1:900000/1‑150 (MQ=255)
                                                                                                                               cttCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTAtt                      <  1:3177840/149‑1 (MQ=255)
                                                                                                                                ttCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTc                     <  1:3012805/150‑1 (MQ=255)
                                                                                                                                  cgcgCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACAACCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTTTACTTTGTCATCTTACTAAAAAGGCGCCGAAGCGCCTTTAGAAAAAAGTCGAAACACTGAATTACTTGTATTCGc                   >  2:3241376/1‑150 (MQ=255)
                                                                                                                                   gcgcCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCt                  >  2:3050091/1‑150 (MQ=255)
                                                                                                                                       cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc              <  1:1668614/150‑1 (MQ=255)
                                                                                                                                          tGCCGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACATTGTCAACTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAAATACTGGTATTCGCTAATCGGt           >  2:4098988/1‑150 (MQ=255)
                                                                                                                                               gCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCa      <  2:3619192/150‑1 (MQ=255)
                                                                                                                                               gCCTGTTCCCCACACCCTAACCCACTACCCAAAGGGGCCAGGGGACTGTCTGAGCACTTTTGTACGGTGTCAACTTACTAAAAAGGCGCCGAAGCGCCATTAGAAAATAGTAGCATCAGTGAATTAAATGTATTAGCGAATCGGTACGCa      >  2:3219117/1‑150 (MQ=255)
                                                                                                                                               gCCTGTTCCCCACACCCAAACCCACACCCCAAAGGGGCGAGGGGACAGTCTGAGCACTTTTGTACATTGTCAACTGAATAAAAAGGCGCCGAACCGCCTTTAGAAAATAGTCGAATCAGTGAAATACAGGTATTCGCTAAACGGTAAGCa      >  2:3229860/1‑150 (MQ=255)
                                                                                                                                               gCCTGGACCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCa      <  1:2506114/150‑1 (MQ=255)
                                                                                                                                                ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg     <  2:3958285/150‑1 (MQ=255)
                                                                                                                                                ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg     <  1:3724893/150‑1 (MQ=255)
                                                                                                                                                  tGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGACTCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa   >  2:2359715/1‑150 (MQ=255)
                                                                                                                                                  tGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa   <  2:545428/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAG  >  NC_000913/3045902‑3046198

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: