Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,007,105 T→G 23.2% V317G (GTG→GGG)  amyA → cytoplasmic alpha‑amylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,007,1050TG23.2% 121.1 / 2.7 68V317G (GTG→GGG) amyAcytoplasmic alpha‑amylase
Reads supporting (aligned to +/- strand):  ref base T (30/18);  new base G (5/14);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTATAAAGAGTGGATCGAACACGTACAGGAAGTTGCGCCAAAGCCGCTGTTTATTGTGGCGGAGTACTGGTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGA  >  NC_000913/2006885‑2007333
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      
ttttATAAAGAGTGGATCGAACACGTACAGGAAGTTGCGCCAAAGCCGCTGTTTATTGTGGCGGAGTACTGGTCGCATGAAGTTGATAAGCTGCAAAAGGATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCCAGAAGCATCGCGAATGGGGCGCGATTCAGCAATGACGCAGCGTTTTACGGGGAACTTAGGAGAAGCCggacattttc                                                                                                                                                                                                                       >  1:426608/1‑229 (MQ=255)
              gATCGAACACGTACAGGAAGTTGCGCCAAAGCCGCTGTTTATTGTGGCGGAGTACTGGTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCCGATTTTCCCGGGGACATTAGTGGAAGCCGATCCTTTCCACCCCGGGAcccccg                                                                                                                                                                                                       >  1:64565/1‑234 (MQ=255)
              gATCGAACACGTACAGGAAGTTGCGCCAAAGCCGCTGTTTATTGTGGCGGAGTACTGGTCGCATGAAGTTAATAAGCCGCCAACGTATATTGGTCAGGTGGAAGGCAAAACCAAGCTGTTTGATGCACCGCTGAAGCTGAAATTTCATTACACATCGCGCATGGGGCGCGGCTACGACAGTACCCAGACTTTCACGGGGACATTAGGGGAAGACGAACCCTTCCCCCCCGCGAccccgg                                                                                                                                                                                                       >  2:398875/1‑234 (MQ=255)
                    acacGTTAAGGAATTTGCGCCAAAGCCGCTTGTTATTGTGGCGGGTGTCTGGGCGGAAGAAGTTGATAAGCTGCAACCGGTTTTTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCa                                                                                                                                                                                                 <  2:60160/239‑1 (MQ=255)
                        tttcAGGAATTTGCGCCAAATCCGTTTTTTATTGTGGCGGTGTTCGGGTCGCATGAATTTGATAAGCTGCAAACGTTTATTGAGTAGGTGGAAGGGAAAACCATGCTGTTTGATGCGCCGCTGGAGATGAAATCCCATGAAGCGTCGCGCATGGGGCGCGTCTACGACATGACGCAGATTTTCACGGGTACATTAGGGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCa                                                                                                                                                                                             <  2:452363/236‑1 (MQ=255)
                           cgggagatggtgataagggacaggTTTATTGTGGCGGAGTACTGGTCGAATGAAGTTGATAAGCTGAAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACg                                                                                                                                                                                           <  1:379710/215‑1 (MQ=255)
                                  ttGCGCCAAAGCCGCTGTTTATTGTGGCGGAGTACTGGTCGCATGAAGTTGATCAGCTGCAAACGTATATTGATCAGGAGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATTAAATTCCATGAAACATCGCGCATTGGGAGCGACTACGACATGACGAAGCATTTACAGGGGACATTAGGTGAAGCCGAACATTTCCACCCCGTGACGCCCGTTGCAAATAACAACACCCa                                                                                                                                                                                    >  2:548954/1‑238 (MQ=255)
                                  ttGCGCCAAAGCCGCTGTTTATTGTGGCGGAGTACTGGTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCCGATTTTAACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCCCGTTGCCAATCCCGACACCCa                                                                                                                                                                                    >  1:790683/1‑238 (MQ=255)
                                                        tgggggggAGTAATGGTCGGATGAAGTTGAGAAGCTGCAAACGTATATTGACCAGGTGGAAGGCAAAACCATGCTGTTGGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAgcgc                                                                                                                                                             <  1:548954/233‑1 (MQ=255)
                                                           ggCGGAGTACTGGTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGGAACTTAGTGGAAGCCGGTCCTTTCCACGCCGTGGCGCTCGTTGCCAATCACGACAACCAACCGGTGGAAGACCCCGAAGAGCCgg                                                                                                                                                          >  2:545167/1‑239 (MQ=255)
                                                               gAGTACTGGTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGCCTACGACATGACGCAAATTTTCACGGGGACATTAAGGGAAGCGGACCCTTTCCCGGCCAGAACGCCCGGTGCAAATCACGACACCCAACCGTTCCAAGCCCTCGAGACGCCGGTcaa                                                                                                                                                      >  2:171807/1‑237 (MQ=255)
                                                                      ggTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCgggg                                                                                                                                               >  1:396050/1‑239 (MQ=255)
                                                                       gTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCCCGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGt                                                                                                                                              >  1:174380/1‑239 (MQ=255)
                                                                       gTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAAcctcggt                                                                                                                                              >  1:373224/1‑234 (MQ=255)
                                                                          gcttGAATTTAATAAGGGTCAAACGTATATTGATCAGGTGGAAGGCAAACCCAGGCTGTTTGGTGGCCCGCTGAAGTGGAAGTTCCATAAAGCTCGGCGGAGGGGGGGCGCTTTGGCCATGACGGAGTTTTCACCGGGTACATAAGGGGAACCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCCCGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTa                                                                                                                                           <  2:804598/236‑1 (MQ=255)
                                                                           caggAAGTTTAATAGTTGAAAACGTTGATTGATGAGGTTGAAGGCAAAACCGTGCGGTTTGTGGCGCCGTGGCAGATGAAATTCCATGAAGGCTCGCGAGTGGGGCGTGACTTCGACATGACGCAGATTTTCACGGGTCCATTAGTGGAAGCCGATCCTTTCCACGCTGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTa                                                                                                                                           <  2:395534/235‑1 (MQ=255)
                                                                                      tAAGCTCCAAACGTATTTTGAACAGGTGGAAGGCAAAACCATGCTTTTTGAGGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGt                                                                                                                              <  1:545167/239‑1 (MQ=255)
                                                                                          cTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCACCCCTTGCAAGCCCTCGAAGGCCCGGTCGACCCGTGGTTTAAACCCCAGGCGTATg                                                                                                                           >  1:226558/1‑238 (MQ=255)
                                                                                            ggaacggtgttgtggtgaGGGGGAAGGAAAACCCAGGGTGGTTGGCGCGCCGCTGCAGGTGAAGTTCAGTGAAGCATCGCGGAGGGGGCGCGACTGCGACAGGCAGAAGGTGTTCACGGGTACATTAGGGGAAGCCGATCCTTCCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGcc                                                                                                                         <  2:655541/221‑1 (MQ=255)
                                                                                            gCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGCCATGACGCAGATTTTCACGGGTCCATTAGGGGAGACCGATCCTTTCAAGCCCGTGACGCCCGTTGCCAATCACCACACCCAACCGTACAAAGCCCTCCAAGCGCCGGTGGCAACGTGTGTTAAACCGCGGGCGAAAGcca                                                                                                                        >  2:246139/1‑238 (MQ=255)
                                                                                             cAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGGAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCa                                                                                                                                                                           >  2:422786/1‑188 (MQ=255)
                                                                                             cAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCa                                                                                                                                                                           <  1:422786/188‑1 (MQ=255)
                                                                                                     tatTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGATACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCCGGCGTATGCCTTAATTttgt                                                                                                               >  1:251727/1‑239 (MQ=255)
                                                                                                        tGACCAGGTGGAAGGCAAAACCTTCCTGTTGGAGGCGCCCGTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACg                                                                                                                                                                                           <  2:507183/161‑1 (MQ=255)
                                                                                                            cAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGTAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAACGGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGa                                                                                                        >  1:902714/1‑239 (MQ=255)
                                                                                                            cAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGa                                                                                                        >  1:10989/1‑239 (MQ=255)
                                                                                                            cAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGGGGTTTAAACCGCTGGCGTATCCCTTAATTTTTTTGCGGGa                                                                                                        >  2:831211/1‑239 (MQ=255)
                                                                                                            cAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGAGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACCCCGTGACGCTCGTTGCCAATCACGACACCCACCCGTTGCAAGCCCCCCAAACGCCGGTCCAACCCGGGTTTACACCGCCGGCGGATGCCTTAATTTTGTTGGGGGa                                                                                                        >  2:429743/1‑239 (MQ=255)
                                                                                                             aGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGCACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGaa                                                                                                       >  1:810937/1‑239 (MQ=255)
                                                                                                               gTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGACCCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGaaaa                                                                                                     >  1:342281/1‑239 (MQ=255)
                                                                                                               gTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGCATGCCTTAATTTTGTTGCGGGaaaa                                                                                                     >  1:820377/1‑239 (MQ=255)
                                                                                                               gTGGAAGGCAAAACCATGCTGTTTGATGCCCCGCTGCAGATGAAATTCCATGAACCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTATAGGAAGCCGATCCTTTCAACGCCGTGACGCTCGTTGCCAATCACGAAAACCCACCGCTGCAAGCCCTCAAAGCCCAGGTCGAACCGTAGTTTAAACCGCTCGCCTATTCCCTCATTTTTTTTCTGAaaaa                                                                                                     >  2:783318/1‑239 (MQ=255)
                                                                                                                  gAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAAcc                                                                                                                                                                                 >  1:507183/1‑161 (MQ=255)
                                                                                                                    aGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACCCCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGACCCGTGGTTTAAACCCCTGCGCGATGCCTTAATTTTGTTGGGGGAAAAAGGCg                                                                                                >  2:201425/1‑239 (MQ=255)
                                                                                                                        aaacCCAGGCGGTTTGATGCGCCCGTGCAGATGAAATTCCATGAAGTATCGCGCATGGGGGGCGACTACGACAGGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTc                                                                                             <  1:429743/238‑1 (MQ=255)
                                                                                                                        aaaaCCTTGCTTGTTGGTGGCCCGCGGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGAGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTc                                                                                             <  1:758216/238‑1 (MQ=255)
                                                                                                                         aaaCCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGGGTTCCt                                                                                           >  1:72853/1‑239 (MQ=255)
                                                                                                                            ccATGCGTTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCAGGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCg                                                                                        <  1:398875/239‑1 (MQ=255)
                                                                                                                               tGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGGAAATGGCGTTCCTTggggt                                                                                     >  2:1063562/1‑237 (MQ=255)
                                                                                                                                gCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGGGCCGGTCGAACCGTGGGTTAAACCGCGGGCGTATGCCTTAATTTTTTTGCGGGAAAATGGCGATACCTCGGttt                                                                                    >  2:968877/1‑237 (MQ=255)
                                                                                                                                      ttgcgccccccttcaggggAAATTCCATTAAGCCTCCGCCTTGGGGCGCGCCGCCGCAATGTCCGATATTTTCACGGGTACATTGGGGGAACCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCCCGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATc                                                                              <  2:342281/221‑1 (MQ=255)
                                                                                                                                      tgtggctccgtggcaagtGAAATTCCAGGAAGCATGGCGCATGGGGGCCGACTACGACATGACGCAGATTTCCGCGGGTACATTAGTGGAAGCCGTTCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATc                                                                              <  2:72853/222‑1 (MQ=255)
                                                                                                                                        gtgcccccctcgaaaTGAAATCCCTGGAAGCACCGCCCATGGGGCGCCTATCCGCCATGCCGCAGATTTTCACGGGTACATTAGGGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCg                                                                            <  2:820377/225‑1 (MQ=255)
                                                                                                                                        gggccccccgggaagtgtaaTTCCATAAGGCACGGCGCTGGGGGGCAGACTACGCCATGACGCAGATTTTCAGGGTTCAATTAGGGGAGGCCGATCCTGTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCg                                                                            <  2:823726/221‑1 (MQ=255)
                                                                                                                                        aTGCGCCGCTGGAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGTTACATTAGGGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCg                                                                            <  1:691312/239‑1 (MQ=255)
                                                                                                                                         ggccgctcgccagtggAAATTCCTTCGAGGACCGCGCGGGGGGCGCACCGACGACAAGCCTCAGATTTGCAGGGGTACATTGTTGGACGCCCACCCGTCCCACGCCGGGACGCTCGTTGCCAGTCCCGACACCAAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGGTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCg                                                                            <  2:704128/223‑1 (MQ=255)
                                                                                                                                              cgcTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGGATTCTATCCGGAcc                                                                        >  1:699906/1‑237 (MQ=255)
                                                                                                                                                aaaaggaacagatcccATGAAGCCGCGCGGTTGGGGGGCGACTACGACATGACGCAGATTTCCACGGGTACATTAGGGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTa                                                                    <  2:916836/225‑1 (MQ=255)
                                                                                                                                                  gCAGATGAAATTCCATGAAGCATCGCGCATGGGGCCCGACTACGACATGACGCAGATTTGCACGGGTACATTAGTGGAAGCCGATCCTTCCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACg                                                                  <  1:1063562/239‑1 (MQ=255)
                                                                                                                                                       tGAAATTCCATGGAAGCTCGCGCGTGGGGCGCGACTCCGACATGGCGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTgcg                                                             <  1:797585/239‑1 (MQ=255)
                                                                                                                                                          aaTTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGCGGTTTAAACCGCTGGGGTTTGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGAATTCTATCCGGCCCTCTCCGGGGCCCAt                                                          >  2:285921/1‑239 (MQ=255)
                                                                                                                                                                aTGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGag                                                    >  1:27534/1‑238 (MQ=255)
                                                                                                                                                                  gaaacccgcggcgaTGGGGCGCGACTACGAATGGACCCAGATTTTCCCGGGTACATTAGGGGAAGCCCATCCTTGCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAg                                                   <  2:1079523/226‑1 (MQ=255)
                                                                                                                                                                       cATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCAt                                                          >  1:916836/1‑227 (MQ=255)
                                                                                                                                                                        aTCGCCGATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGGGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCg                                             <  1:201425/239‑1 (MQ=255)
                                                                                                                                                                          cgcgCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGGATTCTATCCGGGCCTCTACGGTGCCAATTACAAAAGTGCCggt                                           >  1:1049807/1‑239 (MQ=255)
                                                                                                                                                                            cgcATGGGGCGCTACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAgcgc                                                                                                                                                             <  1:777980/124‑1 (MQ=255)
                                                                                                                                                                            cgcATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAgcgc                                                                                                                                                             >  2:777980/1‑124 (MQ=255)
                                                                                                                                                                               aTGGGGCCGGACTACGACATGACGAAGATTTTCACGGGTACATTAGTGGAAGCCGTTCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGgtggtg                                       <  1:731389/238‑1 (MQ=255)
                                                                                                                                                                                      ggccgactagctcattgccgaagttttCCGGGGTAAATATGGGGAACCGGCTCCTTCCCAACCCGTGACTCTAGTTGCCAATCACGACCCCCAACCGTTGCAAGCCCTGGAAGCGCCGGTCGAACCGTAGTTTAACCCGCGGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACggg                                  <  2:810937/215‑1 (MQ=255)
                                                                                                                                                                                       gctactcgcccATGCCGCATCTTTTCAGGGTTCAATTGTGGGAACCGCATCCTTCCCCCGCCGTGCAGCTCGTGGCCAAGCCAGACATCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGTTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCaa                               <  2:226558/229‑1 (MQ=255)
                                                                                                                                                                                         gactcggACATGAGGGAGTTTTCACGGGGTCCATTAGGGGAAGCCGCTCCTTCCCACGCCGTGACGCTCGTTGCCAATCCCGACCCCCAACCGTGGCAAGCCCTGGAAGTGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAAc                             <  2:174380/233‑1 (MQ=255)
                                                                                                                                                                                         gACTACGACATGACGCATATTTTCACGGGTACATTAGGGGAAGCCGATCCTTGCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAAc                             <  1:18066/239‑1 (MQ=255)
                                                                                                                                                                                         gACTACGACATGACGCAGATTTTCACGGGTACATTAGt                                                                                                                                                                                                                                        <  2:396340/38‑1 (MQ=255)
                                                                                                                                                                                              cGACTTGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGAATTCTATCCGGACCCCTAAGGTGCGCATTACGAAAAGGTCGGTGGGGACGGGCAAACCTAt                         >  1:818997/1‑238 (MQ=255)
                                                                                                                                                                                                ataTGGCGGAGGTTTGCGCGGGGACAGTAGGTGGAGCCGATCCTGTCCACCCCGTGACGCTCGTTGCAAATCACGACACCCAACCGGTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTAGCCGGACCGCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCg                      <  2:909595/237‑1 (MQ=255)
                                                                                                                                                                                                   tGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGAt                                                                                                                                                                                                                             >  1:396340/1‑38 (MQ=255)
                                                                                                                                                                                                            aTTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACctctc                                                                    >  1:987051/1‑181 (MQ=255)
                                                                                                                                                                                                                        gTACATTAGTGGAAGCCGTTCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGa  <  1:831211/238‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      
TTTTATAAAGAGTGGATCGAACACGTACAGGAAGTTGCGCCAAAGCCGCTGTTTATTGTGGCGGAGTACTGGTCGCATGAAGTTGATAAGCTGCAAACGTATATTGATCAGGTGGAAGGCAAAACCATGCTGTTTGATGCGCCGCTGCAGATGAAATTCCATGAAGCATCGCGCATGGGGCGCGACTACGACATGACGCAGATTTTCACGGGTACATTAGTGGAAGCCGATCCTTTCCACGCCGTGACGCTCGTTGCCAATCACGACACCCAACCGTTGCAAGCCCTCGAAGCGCCGGTCGAACCGTGGTTTAAACCGCTGGCGTATGCCTTAATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGA  >  NC_000913/2006885‑2007333

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: