Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,050,713 A→T 5.6% intergenic (‑46/+402) gcvT ← / ← ubiI aminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,050,7130AT 5.6% 299.7 / 3.6 98intergenic (‑46/+402)gcvT/ubiIaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (46/46);  new base T (4/2);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGC  >  NC_000913/3050564‑3050862
                                                                                                                                                     |                                                                                                                                                      
tcgATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGaa                                                                                                                                                        >  2:1406595/1‑150 (MQ=255)
    tttGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTa                                                                                                                                                    <  1:539179/150‑1 (MQ=255)
      tGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATc                                                                                                                                                  <  1:754146/150‑1 (MQ=255)
      tGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATc                                                                                                                                                  <  2:1704218/150‑1 (MQ=255)
      tGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATc                                                                                                                                                  <  2:1695583/150‑1 (MQ=255)
      tGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATc                                                                                                                                                  <  2:1484935/150‑1 (MQ=255)
           aCCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTTCAAAGGAGTCTGTTGTGCCATCTTGTCCTCATGGAATAAGCGGGGCTGACCACTTTTTCAAGGTGAAAGTATcaccac                                                                                                                                             >  2:1143048/1‑150 (MQ=255)
           aCCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAATTTATcaccac                                                                                                                                             <  2:864273/150‑1 (MQ=255)
            ccGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTTACAACTTTTTCATGGTGAAATTATCACCACg                                                                                                                                            >  2:1228409/1‑150 (MQ=255)
            ccGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACg                                                                                                                                            >  1:573351/1‑150 (MQ=255)
            ccGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACg                                                                                                                                            >  1:347945/1‑150 (MQ=255)
               tAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCAGCTAAGCGTGTGTTGTTCTTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGATAATTTCACCATGaaa                                                                                                                                         <  2:1805868/150‑1 (MQ=255)
                aaTGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGTGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAAc                                                                                                                                        <  1:705560/150‑1 (MQ=255)
                   gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCa                                                                                                                                     <  1:513566/150‑1 (MQ=255)
                   gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCa                                                                                                                                     <  1:506262/150‑1 (MQ=255)
                   gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCa                                                                                                                                     <  1:715825/150‑1 (MQ=255)
                   gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCAACTAGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCa                                                                                                                                     <  1:709051/150‑1 (MQ=255)
                    cAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAg                                                                                                                                    >  1:1110334/1‑150 (MQ=255)
                    cAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAg                                                                                                                                    <  1:1569312/150‑1 (MQ=255)
                    cAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAg                                                                                                                                    >  1:1102972/1‑150 (MQ=255)
                     aGCGGCATCATCCAGCCGTGGAAATCCACCATGCCAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTTTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCGTGGTGATAATTTCACCATGaaaaagttg                                                                                                                                   <  1:377792/150‑7 (MQ=255)
                          cctcatcCAGCCCTGGAAAGCCACCATGCGAGAGACGCAAAGCGTGTGTTGTTCGTACATAGGACTCTGTAGTGCCAGCATGTCCTCATTGAATTAGATGGGCAGACAACTTTTTCATGGTCAAATTTTCACCACGAAACGCAGCATCgg                                                                                                                              >  2:1062519/3‑150 (MQ=255)
                          catcatCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCgg                                                                                                                              <  2:688361/150‑1 (MQ=255)
                            tcatcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAg                                                                                                                            >  2:1757636/1‑150 (MQ=255)
                             catcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTCAAATTATAACAACGAAACCCAGCATCGGAGc                                                                                                                           >  1:1464515/1‑150 (MQ=255)
                             catcCAGACGTGTAAATCCACCATTCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGc                                                                                                                           <  1:291853/150‑1 (MQ=255)
                              atcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCCCGAAACCCAGCATCGGAGcc                                                                                                                          >  2:1529222/1‑150 (MQ=255)
                                ccAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGCACCCCAGCATCGGAGCCAc                                                                                                                        >  2:576033/1‑150 (MQ=255)
                                 cAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAGCTTTTTCATGGTGATATTATCACCACGAAACCCAGCATCGGAGCCCCt                                                                                                                       >  2:1245629/1‑150 (MQ=255)
                                    ccGTGTAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTccc                                                                                                                    <  1:380109/150‑1 (MQ=255)
                                    ccGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGCGCCACTccc                                                                                                                    >  1:1209211/1‑150 (MQ=255)
                                      gTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCgg                                                                                                                  <  1:24047/150‑1 (MQ=255)
                                      gTGGAAATACACCATGCGAGCGCCGCAAATCGTGTGGTGTTCGTACAAAGGAGTCTGTTGTGCCCTCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGTTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCgg                                                                                                                  <  1:182117/150‑1 (MQ=255)
                                        ggAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTc                                                                                                                <  2:868638/150‑1 (MQ=255)
                                          aaaTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACCACTTTTTCATGGTGAAATTATCACCACGGAACCCAGCATCGGCGCCACTCCCGGTccc                                                                                                              >  1:175015/1‑150 (MQ=255)
                                          aaaTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTccc                                                                                                              <  2:1721351/150‑1 (MQ=255)
                                             tCCACCATGCGAGCTCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCaa                                                                                                           <  2:1389191/150‑1 (MQ=255)
                                              ccaccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACACCTTTTTCATGGTGATATTATCACCACGATACCCAGCCTGGGAGCCACTCCCGGTCCCCAAc                                                                                                          >  2:1596227/1‑150 (MQ=255)
                                               caccaTGCGAGCGCCGAAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                                                         >  1:1449031/1‑150 (MQ=255)
                                                       gAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGtt                                                                                                 >  1:1527021/1‑150 (MQ=255)
                                                            ccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtt                                                                                            <  1:1308800/150‑1 (MQ=255)
                                                             cgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCttt                                                                                           <  1:597587/150‑1 (MQ=255)
                                                             cgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCttt                                                                                           <  2:49782/150‑1 (MQ=255)
                                                              gcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtttt                                                                                          >  1:794966/1‑150 (MQ=255)
                                                               cAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTg                                                                                         <  2:1330533/150‑1 (MQ=255)
                                                                 aaGCGTGTGTTGTTCGTACAAAGGAGTCTTTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCAAGAAAACCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGcc                                                                                       <  2:325498/150‑1 (MQ=255)
                                                                 aaGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGcc                                                                                       <  2:1292827/150‑1 (MQ=255)
                                                                  aGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTGTCTTTTGCCt                                                                                      >  2:705936/1‑150 (MQ=255)
                                                                  aGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCt                                                                                      >  2:1615232/1‑150 (MQ=255)
                                                                     gtgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGaa                                                                                   <  2:1400959/150‑1 (MQ=255)
                                                                      ggtgtTGGTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                  <  1:862807/149‑1 (MQ=255)
                                                                         gttgttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTa                                                                               >  2:1139777/1‑150 (MQ=255)
                                                                           tgttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacc                                                                             <  1:47964/150‑1 (MQ=255)
                                                                            gttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacca                                                                            >  1:12963/1‑150 (MQ=255)
                                                                                 tACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGCCAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGaa                                                                       >  2:1402628/1‑150 (MQ=255)
                                                                                  aCAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACTCAATCGTTCTCTTTTGCGTGAACTTACCACCGaaa                                                                      >  2:168674/1‑150 (MQ=255)
                                                                                  aCAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGaaa                                                                      >  2:525468/1‑150 (MQ=255)
                                                                                    aaaGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACa                                                                    <  2:442507/150‑1 (MQ=255)
                                                                                       ggAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                                                                 <  1:812259/150‑1 (MQ=255)
                                                                                         aGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTg                                                               <  2:1076850/150‑1 (MQ=255)
                                                                                          gTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCTGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGt                                                              >  1:79897/1‑150 (MQ=255)
                                                                                           tCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACATCTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCGACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGtt                                                             <  1:1547893/150‑1 (MQ=255)
                                                                                            cTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTa                                                            <  1:361075/150‑1 (MQ=255)
                                                                                              gttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                          <  2:1078649/150‑1 (MQ=255)
                                                                                               ttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAcc                                                         <  2:1208183/150‑1 (MQ=255)
                                                                                               ttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAcc                                                         <  2:1554403/150‑1 (MQ=255)
                                                                                                    ccATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAg                                                    <  2:1449031/150‑1 (MQ=255)
                                                                                                     cATCTTGTCCTCATTTAATAATCGTGGATGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAgg                                                   <  2:257171/150‑1 (MQ=255)
                                                                                                     cATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAgg                                                   <  1:1406595/150‑1 (MQ=255)
                                                                                                           gTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAACCCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAAGCGGTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAt                                             >  1:1798440/1‑150 (MQ=255)
                                                                                                            tCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAtt                                            >  2:857956/1‑150 (MQ=255)
                                                                                                             ccTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTg                                           >  1:678582/1‑150 (MQ=255)
                                                                                                             ccTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTg                                           >  2:850898/1‑150 (MQ=255)
                                                                                                              cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGa                                          >  1:1772903/1‑150 (MQ=255)
                                                                                                                cATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc                                        >  2:1444835/1‑150 (MQ=255)
                                                                                                                    gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGAtcatca                                   >  2:1515410/1‑150 (MQ=255)
                                                                                                                    gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGAtcatca                                   >  2:338297/1‑150 (MQ=255)
                                                                                                                    gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCACCAACGCAATCGGTCTCTTTTGCCTGAACTGACCACAGAAACAGCCTGTTAACCATAAGGTAAAATTGCtcatca                                   >  1:1142466/1‑150 (MQ=255)
                                                                                                                       tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACAt                                <  1:466807/150‑1 (MQ=255)
                                                                                                                        aaGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACATTATGGTAAAATTGATAATCACAtt                               >  2:459982/1‑150 (MQ=255)
                                                                                                                         aGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACtttt                              >  2:1148914/1‑149 (MQ=255)
                                                                                                                          gCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAg                             <  1:683074/150‑1 (MQ=255)
                                                                                                                           cGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGc                            >  2:1460951/1‑150 (MQ=255)
                                                                                                                            tgggCTGACAACTTTTTCATGGTGAAATTATCACCACTAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCt                           <  1:1876278/149‑1 (MQ=255)
                                                                                                                            ggggCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACAACCGAAACAGACTGTTAACCATAAGGTAAAATTGATAATCACATTAGCt                           >  1:276487/1‑150 (MQ=255)
                                                                                                                                cTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATg                       >  1:1798765/1‑150 (MQ=255)
                                                                                                                                cTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGATCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATg                       >  2:1453627/1‑150 (MQ=255)
                                                                                                                                  gACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGt                     >  2:11395/1‑150 (MQ=255)
                                                                                                                                       cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGCAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGATTATGGTTaaaa                >  1:258007/1‑150 (MQ=255)
                                                                                                                                       cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaa                >  1:1784270/1‑150 (MQ=255)
                                                                                                                                         ttttCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGCTAATGTGATGATCAATTTTACCTTATGGTTaacagt              <  2:573351/150‑5 (MQ=25)
                                                                                                                                         gtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTATTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATTAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaaaa              >  1:64886/2‑150 (MQ=255)
                                                                                                                                         gtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCAGGTCCCCAACGCAATCGTTATATTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATTAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaaaa              >  1:996775/2‑150 (MQ=255)
                                                                                                                                           ttCATGGTGAAATTATCACCACGAACCCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTAAGGTTAAAAAATg            >  2:1653337/1‑150 (MQ=255)
                                                                                                                                            tCATGGTTAAATTATCACCACGAAACCCAGCATCGGAGACACTCACGGTCCCCAACGCAATCGTTCTCTTTTGACTTAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGc           <  1:543171/150‑1 (MQ=255)
                                                                                                                                            tCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGc           <  1:1757636/150‑1 (MQ=255)
                                                                                                                                             cATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCa          >  1:7424/1‑150 (MQ=255)
                                                                                                                                               tGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaa        >  1:861886/1‑150 (MQ=255)
                                                                                                                                                ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa       >  1:1587749/1‑150 (MQ=255)
                                                                                                                                                     aaTTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATcgc  <  1:1708481/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGC  >  NC_000913/3050564‑3050862

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: