Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,047,996 G→A 5.8% A349V (GCC→GTC)  gcvP ← 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,047,9960GA 5.8% 2523.5 / 92.0 770A349V (GCC→GTC) gcvPglycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base G (355/369);  new base A (21/24);  total (376/394)
Fisher's exact test for biased strand distribution p-value = 8.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

Only 100 of 781 total aligned reads displayed.

TTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCTTTCGATACACCGATAATACG  >  NC_000913/3047848‑3048144
                                                                                                                                                    |                                                                                                                                                    
tttGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGGCCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACgg                                                                                                                                                      <  1:1674385/149‑1 (MQ=255)
  tGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCa                                                                                                                                                    <  4:5753289/149‑1 (MQ=255)
         cACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGc                                                                                                                                             >  3:852963/1‑149 (MQ=255)
          aCTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGTTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCt                                                                                                                                            >  2:562895/1‑149 (MQ=255)
          aCTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCt                                                                                                                                            >  3:600379/1‑149 (MQ=255)
           cTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTg                                                                                                                                           <  2:5186173/149‑1 (MQ=255)
            ttCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTgg                                                                                                                                          >  2:4040082/1‑149 (MQ=255)
             tCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGc                                                                                                                                         <  4:3994924/149‑1 (MQ=255)
                    tcaAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGttt                                                                                                                                  >  3:288650/2‑149 (MQ=255)
                    acaAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGACATACAGGCTGGCGATGttt                                                                                                                                  <  2:1164279/149‑1 (MQ=255)
                      aaGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGc                                                                                                                                <  3:3962909/149‑1 (MQ=255)
                      aaGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGc                                                                                                                                <  1:4307991/149‑1 (MQ=255)
                        ggTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCca                                                                                                                              <  3:703259/149‑1 (MQ=255)
                        ggTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGGTTGCca                                                                                                                              >  3:2552184/1‑149 (MQ=255)
                           gtCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCcagca                                                                                                                           >  2:7157331/1‑149 (MQ=255)
                            tCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCcagca                                                                                                                           >  1:1925006/1‑148 (MQ=255)
                              gAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTa                                                                                                                        >  4:520552/1‑149 (MQ=255)
                                 aTAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCt                                                                                                                     <  1:4089836/149‑1 (MQ=255)
                                 aTAGTGCGCATGGCACAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCt                                                                                                                     <  1:544213/149‑1 (MQ=255)
                                    gTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTggg                                                                                                                  <  1:6719475/149‑1 (MQ=255)
                                      gcgcATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGaa                                                                                                                <  2:6692877/149‑1 (MQ=255)
                                            ggCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACaa                                                                                                          <  2:2183031/149‑1 (MQ=255)
                                              cgcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAAt                                                                                                        <  2:4556206/149‑1 (MQ=255)
                                              cgcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAAt                                                                                                        <  4:6296210/149‑1 (MQ=255)
                                              cgcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAAt                                                                                                        <  2:6925224/149‑1 (MQ=255)
                                               gcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATg                                                                                                       >  1:3822614/1‑149 (MQ=255)
                                               gcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAACAGTACCTGGGAAGTACAAATg                                                                                                       >  3:4378089/1‑149 (MQ=255)
                                                 aGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGtt                                                                                                     <  3:4839021/149‑1 (MQ=255)
                                                  gTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTg                                                                                                    <  4:4249059/149‑1 (MQ=255)
                                                  gTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTg                                                                                                    <  3:1519611/149‑1 (MQ=255)
                                                    ttCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGa                                                                                                  <  4:1409088/149‑1 (MQ=255)
                                                    ttCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGa                                                                                                  <  1:3508359/149‑1 (MQ=255)
                                                    ttCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGa                                                                                                  <  4:4610580/149‑1 (MQ=255)
                                                    ttCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGa                                                                                                  <  2:5984783/149‑1 (MQ=255)
                                                    ttCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGa                                                                                                  <  1:6312184/149‑1 (MQ=255)
                                                     tCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAAt                                                                                                        <  1:1032833/142‑1 (MQ=255)
                                                        gACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCcagca                                                                                                                           <  2:427464/120‑1 (MQ=255)
                                                        gACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTg                                                                                                    <  3:3528206/143‑1 (MQ=255)
                                                        gACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTc                                                                                              <  2:471158/149‑1 (MQ=255)
                                                         aCCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGaa                                                                                                                <  4:4477774/130‑1 (MQ=255)
                                                          ccTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGc                                                                                            <  1:2991397/149‑1 (MQ=255)
                                                           cTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTc                                                                                              <  4:6870310/146‑1 (MQ=255)
                                                           cTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                                                           <  3:5215868/149‑1 (MQ=255)
                                                           cTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                                                           >  1:3404874/1‑149 (MQ=255)
                                                           cTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                                                           >  2:986660/1‑149 (MQ=255)
                                                            tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                                                          <  4:6581481/149‑1 (MQ=255)
                                                            tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                                                          >  4:3027958/1‑149 (MQ=255)
                                                            tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                                                          >  4:179234/1‑149 (MQ=255)
                                                                        gCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCcagca                                                                                                                           >  3:5060601/1‑104 (MQ=255)
                                                                         cccGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGat                                                                             <  3:6354954/149‑1 (MQ=255)
                                                                         cccGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGat                                                                             <  3:1707856/149‑1 (MQ=255)
                                                                         cccGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGACATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGat                                                                             <  1:5778416/149‑1 (MQ=255)
                                                                          ccgccgCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTa                                                                                                             >  1:5232269/1‑116 (MQ=255)
                                                                          ccgccgCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGata                                                                            >  1:6452959/1‑149 (MQ=255)
                                                                            gccgccAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGACATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATg                                                                          <  2:3042354/149‑1 (MQ=255)
                                                                                  aGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTc                                                                    >  2:3689845/1‑149 (MQ=255)
                                                                                      tatCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGa                                                                <  2:195623/149‑1 (MQ=255)
                                                                                      tatCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGa                                                                <  2:2416104/149‑1 (MQ=255)
                                                                                       atCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAg                                                               >  4:5620299/1‑149 (MQ=255)
                                                                                        tCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGa                                                                                                  <  1:4971075/113‑1 (MQ=255)
                                                                                         cGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTc                                                             <  4:5956633/149‑1 (MQ=255)
                                                                                          ggTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCt                                                            <  4:1157311/149‑1 (MQ=255)
                                                                                           gTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTgg                                                                                                   <  1:1185769/109‑1 (MQ=255)
                                                                                              aGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCAt                                                        <  2:167593/149‑1 (MQ=255)
                                                                                                 cGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGc                                                     <  2:1375673/149‑1 (MQ=255)
                                                                                                   gTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGACATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCa                                                   >  4:2229117/1‑149 (MQ=255)
                                                                                                      aaTGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATgc                                                <  3:1991286/149‑1 (MQ=255)
                                                                                                       aTGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATgcg                                               >  4:4916450/1‑149 (MQ=255)
                                                                                                             ttAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgcgc                                         >  1:2299214/1‑149 (MQ=255)
                                                                                                              tAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGcgcg                                        <  3:2768359/149‑1 (MQ=255)
                                                                                                              tAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGcgcg                                        <  3:665494/149‑1 (MQ=255)
                                                                                                                 cGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTa                                     >  1:6204007/1‑149 (MQ=255)
                                                                                                                   aTACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTc                                                                                              >  4:1603988/1‑90 (MQ=255)
                                                                                                                   aTACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTAtt                                   <  3:1365524/149‑1 (MQ=255)
                                                                                                                   aTACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTAtt                                   >  2:4963384/1‑149 (MQ=255)
                                                                                                                    tACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTg                                  <  1:250632/149‑1 (MQ=255)
                                                                                                                     aCGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGt                                                                         <  1:4551488/109‑1 (MQ=255)
                                                                                                                     aCGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGc                                 <  2:3527481/149‑1 (MQ=255)
                                                                                                                     aCGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGc                                 <  4:5367227/149‑1 (MQ=255)
                                                                                                                      cGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGcc                                >  1:461844/1‑149 (MQ=255)
                                                                                                                       gTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCa                               <  2:5344155/149‑1 (MQ=255)
                                                                                                                       gTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCa                               <  1:2778570/149‑1 (MQ=255)
                                                                                                                        tttCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGGGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAg                              >  2:6343894/1‑149 (MQ=255)
                                                                                                                          tCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCt                            <  4:741169/149‑1 (MQ=255)
                                                                                                                           cAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGcc                                <  4:3541212/144‑1 (MQ=255)
                                                                                                                            aGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTa                                                                                                             >  4:2512644/1‑66 (MQ=255)
                                                                                                                             ggCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACaaa                                                                                                         <  3:334269/69‑1 (MQ=255)
                                                                                                                             ggCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTg                                                                                                    >  4:2806076/1‑74 (MQ=255)
                                                                                                                             ggCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCa                         >  3:44521/1‑149 (MQ=255)
                                                                                                                              gCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTc                                                                                              <  3:2170729/79‑1 (MQ=255)
                                                                                                                                cAACCGGGCCGTGATAAACGACATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCt                      <  2:33721/149‑1 (MQ=255)
                                                                                                                                 aaCCGGGCCNTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCtt                     >  3:4673306/1‑149 (MQ=255)
                                                                                                                                 aaCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCtt                     <  2:256960/149‑1 (MQ=255)
                                                                                                                                  aCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCttt                    <  3:422384/149‑1 (MQ=255)
                                                                                                                                   ccGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCTTTc                   <  4:3362950/149‑1 (MQ=255)
                                                                                                                                            gATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGcc                                <  1:7223927/127‑1 (MQ=255)
                                                                                                                                            gATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCTTTCGATACACCg          >  2:2927136/1‑149 (MQ=255)
                                                                                                                                             aTAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCTTTCGATACACCGa         <  1:7016838/149‑1 (MQ=255)
                                                                                                                                                aaCGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCTTTCGATACACCGataa      >  3:5519551/1‑149 (MQ=255)
                                                                                                                                                aaCGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCTTTCGATACACCGataa      <  1:3659368/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTGTCGGCCACTTCCACACACAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCCAGCTGCATCTTTCGATACACCGATAATACG  >  NC_000913/3047848‑3048144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: