Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,317 A→T 6.1% intergenic (+59/+23) ycaK → / ← pflA putative NAD(P)H‑dependent oxidoreductase YcaK/pyruvate formate‑lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913950,3170AT 6.1% 291.7 / 6.2 97intergenic (+59/+23)ycaK/pflAputative NAD(P)H‑dependent oxidoreductase YcaK/pyruvate formate‑lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (43/48);  new base T (2/4);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

ACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGCCC  >  NC_000913/950172‑950466
                                                                                                                                                 |                                                                                                                                                     
aCAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGt                                                                                                                                                   >  2:717351/1‑150 (MQ=255)
aCAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGt                                                                                                                                                   <  1:425554/150‑1 (MQ=255)
     ggTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCg                                                                                                                                              <  1:1237902/150‑1 (MQ=255)
       tATTCGATGAGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGCt                                                                                                                                            <  2:1310816/150‑1 (MQ=255)
       gatgccaTGGGGACGAACTTCACGCGAGCTATTATCAGTCGTCATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGGTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCt                                                                                                                                            <  2:1287455/144‑1 (MQ=255)
          tCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTa                                                                                                                                         >  2:728021/1‑150 (MQ=255)
          tCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTa                                                                                                                                         >  1:1656580/1‑150 (MQ=255)
            gATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACg                                                                                                                                       <  1:599315/150‑1 (MQ=255)
            gATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGGGTTAAAAACGCGGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGCtttttt                                                                                                                                       <  1:492972/150‑4 (MQ=255)
             aTGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCtttttta                                                                                                                                      <  1:286007/150‑5 (MQ=255)
             aTGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACgg                                                                                                                                      >  1:201660/1‑150 (MQ=255)
              tGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACggg                                                                                                                                     <  2:1797907/150‑1 (MQ=255)
              tGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACggg                                                                                                                                     <  2:1198893/150‑1 (MQ=255)
               ggggAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACgggg                                                                                                                                    <  1:728021/150‑1 (MQ=255)
                gggAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGa                                                                                                                                   >  2:1774938/1‑150 (MQ=255)
                gggAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGa                                                                                                                                   >  2:169576/1‑150 (MQ=255)
                    agaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAAtt                                                                                                                               >  2:1057013/1‑150 (MQ=255)
                    agaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAAtt                                                                                                                               <  1:554554/150‑1 (MQ=255)
                    agaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAAtt                                                                                                                               >  1:1884205/1‑150 (MQ=255)
                    agaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAAtt                                                                                                                               <  1:1678393/150‑1 (MQ=255)
                       aCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGa                                                                                                                            >  2:1867054/1‑150 (MQ=255)
                       aCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGa                                                                                                                            <  1:717351/150‑1 (MQ=255)
                            aCGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACAtt                                                                                                                       >  2:531216/1‑150 (MQ=255)
                              gcgAGCTATTATCATTCGTTATTATCTCAGGTACTGGATATTGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTAc                                                                                                                     <  2:728650/150‑1 (MQ=255)
                               cgAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTAcc                                                                                                                    <  1:877142/150‑1 (MQ=255)
                                    tattatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATg                                                                                                               >  1:1695318/1‑150 (MQ=255)
                                    tattatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATg                                                                                                               <  1:842721/150‑1 (MQ=255)
                                    tattatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATg                                                                                                               >  1:513571/1‑150 (MQ=255)
                                     attatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCt                                                                                                                   <  2:1153683/145‑1 (MQ=255)
                                     attatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCt                                                                                                                   >  1:1153683/1‑145 (MQ=255)
                                      ttatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGAc                                                                                                             <  1:775826/150‑1 (MQ=255)
                                       tatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGAcc                                                                                                            >  2:738114/1‑150 (MQ=255)
                                       tatCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGAcc                                                                                                            <  1:1306225/150‑1 (MQ=255)
                                          cAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTa                                                                                                         >  2:909113/1‑150 (MQ=255)
                                           aGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTAc                                                                                                        <  1:531216/150‑1 (MQ=255)
                                           aGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTAc                                                                                                        >  1:64090/1‑150 (MQ=255)
                                               gTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCt                                                                                                    <  1:459910/150‑1 (MQ=255)
                                                ttattaTCTCAGGTACGGGATATGGTAGATGCACTATAAGATTTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGTCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTc                                                                                                   <  2:691529/150‑1 (MQ=255)
                                                ttattaTCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTc                                                                                                   >  1:197507/1‑150 (MQ=255)
                                                 tattatCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCa                                                                                                  <  2:695783/150‑1 (MQ=255)
                                                   ttatCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAg                                                                                                <  2:834085/150‑1 (MQ=255)
                                                      tctcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGCAGCGCGGAATAAAAAAGCGGCAACTCAAGAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAAt                                                                                             >  2:1842567/1‑150 (MQ=255)
                                                       ctcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATg                                                                                            <  2:201660/150‑1 (MQ=255)
                                                       ctcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATg                                                                                            >  2:1031336/1‑150 (MQ=255)
                                                       ctcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATg                                                                                            <  1:952095/150‑1 (MQ=255)
                                                        tcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGc                                                                                           >  2:214597/1‑150 (MQ=255)
                                                             tACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTc                                                                                      >  1:514202/1‑150 (MQ=255)
                                                             tACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTc                                                                                      >  2:159762/1‑150 (MQ=255)
                                                                gggATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCAcg                                                                                   >  1:392661/1‑150 (MQ=255)
                                                                     atGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTc                                                                              <  1:528067/150‑1 (MQ=255)
                                                                       ggTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTAc                                                                                                        <  2:1228934/122‑1 (MQ=255)
                                                                       ggTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTAc                                                                                                        >  1:1228934/1‑122 (MQ=255)
                                                                       ggTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGa                                                                                               >  2:1962171/1‑131 (MQ=255)
                                                                       ggTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGa                                                                                               <  1:1962171/131‑1 (MQ=255)
                                                                       ggTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCa                                                                            >  1:1734368/1‑150 (MQ=255)
                                                                          aGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGCCCGTAGTGGGCAAGAATGCGTTTCACGCGTTCCATgg                                                                         >  2:1022490/1‑150 (MQ=255)
                                                                          aGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATgg                                                                         <  1:425130/150‑1 (MQ=255)
                                                                            aTGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGtc                                                                       <  1:1774938/150‑1 (MQ=255)
                                                                               cACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTctct                                                                    >  2:1022206/1‑150 (MQ=255)
                                                                                aCTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCtctt                                                                   <  1:909113/150‑1 (MQ=255)
                                                                                 cTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCAGCTTTACGGGGAAATTAGAACAATACCTTATGAACGAACTGCTCAAGAAAGCCTTTCACGCGTTCCATGGTCtcttt                                                                  >  1:717884/1‑150 (MQ=255)
                                                                                        aTGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGTCTTTCACGCGTTCCATGGTCTCTTTCTTCggt                                                           <  2:1656580/150‑1 (MQ=255)
                                                                                          ttgtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         <  1:463311/149‑1 (MQ=255)
                                                                                          gtgtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         <  2:1695318/150‑1 (MQ=255)
                                                                                              tAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTa                                                     >  1:1475460/1‑150 (MQ=255)
                                                                                               aaaaaCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTaa                                                    >  1:984235/1‑150 (MQ=255)
                                                                                               aaaaaCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTaa                                                    >  1:150344/1‑150 (MQ=255)
                                                                                                       tGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCg                                            <  1:159762/150‑1 (MQ=255)
                                                                                                         tAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAg                                          >  1:1304238/1‑150 (MQ=255)
                                                                                                         tAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAg                                          >  1:492348/1‑150 (MQ=255)
                                                                                                          aGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         >  2:782266/1‑134 (MQ=255)
                                                                                                          aGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         <  1:1295823/134‑1 (MQ=255)
                                                                                                          aGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         <  1:782266/134‑1 (MQ=255)
                                                                                                          aGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         >  2:1295823/1‑134 (MQ=255)
                                                                                                          aGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGt                                         >  2:1161070/1‑150 (MQ=255)
                                                                                                           gCATAATGAAGCGCGGAATAAAAAAGCTTCAACTCAATAAAGTTTCCGCTTTACGTTGAAATTAGAACATTACCTTATGTCCGTACTGCTCAAGAATGCCTTTCACGCGTTCGATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGtt                                        <  2:508588/150‑1 (MQ=255)
                                                                                                                  gAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTActc                                 >  2:1325609/1‑150 (MQ=255)
                                                                                                                   aaGCTCGTAATAAAATAGCGGCAACTCAATAAAGTTGCCGCTTTACGTGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTGCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTActct                                <  2:717884/150‑1 (MQ=255)
                                                                                                                    aGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCtt                               >  1:467909/1‑150 (MQ=255)
                                                                                                                        cgGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCAcc                           <  2:197507/150‑1 (MQ=255)
                                                                                                                        cgGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCAcc                           <  2:1884205/150‑1 (MQ=255)
                                                                                                                          tccccgtaaaGCGGCAACTTTATTGAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCa                         <  1:779138/143‑1 (MQ=39)
                                                                                                                             tAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTg                      <  1:738114/150‑1 (MQ=255)
                                                                                                                                aaaaGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCa                   >  2:66722/1‑150 (MQ=255)
                                                                                                                                  aaGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCAcc                 <  1:214597/150‑1 (MQ=255)
                                                                                                                                     cGGCAACTTTATTGAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCAt              >  2:1926114/1‑150 (MQ=255)
                                                                                                                                     cGGCAACTTTATTGAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCAt              >  2:783446/1‑150 (MQ=255)
                                                                                                                                      ggCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCAtt             <  2:984235/150‑1 (MQ=255)
                                                                                                                                      ggCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCAtt             <  2:1734368/150‑1 (MQ=255)
                                                                                                                                        cAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACATCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtg           <  2:391853/150‑1 (MQ=255)
                                                                                                                                        cAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACATCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtg           >  1:391853/1‑150 (MQ=255)
                                                                                                                                        cAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtg           <  2:392661/150‑1 (MQ=255)
                                                                                                                                        cAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtg           <  1:705855/150‑1 (MQ=255)
                                                                                                                                          aCTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtgtg         <  2:513571/150‑1 (MQ=255)
                                                                                                                                           cTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTgtgt        >  1:848236/1‑150 (MQ=255)
                                                                                                                                             cAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGttt      <  1:1022206/150‑1 (MQ=255)
                                                                                                                                              aaTAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTg     >  1:342156/1‑150 (MQ=255)
                                                                                                                                                 aaaGTTGCCGCTTTACGGGGATATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTGTCGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGccc  <  1:1057013/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGCCC  >  NC_000913/950172‑950466

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