Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,763 A→T 5.4% intergenic (+29/+19) yeaD → / ← yeaE D‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,864,7630AT 5.4% 323.7 / 5.1 110intergenic (+29/+19)yeaD/yeaED‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (56/48);  new base T (2/4);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 114 total aligned reads displayed.

GCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGCT  >  NC_000913/1864614‑1864911
                                                                                                                                                     |                                                                                                                                                    
gCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCaa                                                                                                                                                      >  1:19056/1‑150 (MQ=255)
gCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCaa                                                                                                                                                      >  1:803920/1‑150 (MQ=255)
       gatgGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                                                                               >  2:1486284/1‑150 (MQ=255)
          ggCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                                                                            <  1:940702/150‑1 (MQ=255)
          ggCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                                                                            <  1:241164/150‑1 (MQ=255)
          ggCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                                                                            <  1:262570/150‑1 (MQ=255)
           gCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGca                                                                                                                                           <  2:971129/150‑1 (MQ=255)
           gCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATACATCCGGca                                                                                                                                           >  1:1722021/1‑150 (MQ=255)
            cTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcac                                                                                                                                          <  1:1794233/150‑1 (MQ=255)
            cTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcac                                                                                                                                          <  1:1672526/150‑1 (MQ=255)
              aCAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacac                                                                                                                                        <  2:520135/150‑1 (MQ=255)
              aCAAAACATTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacaca                                                                                                                                       <  1:491417/150‑1 (MQ=255)
               cAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacaca                                                                                                                                       >  1:1164986/1‑150 (MQ=255)
                 aaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTTCCGGGTGCAATGCATCCGGCACacaac                                                                                                                                     >  2:1393150/1‑150 (MQ=255)
                 aaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaac                                                                                                                                     >  2:21776/1‑150 (MQ=255)
                 aaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACacaac                                                                                                                                     <  1:873767/150‑1 (MQ=255)
                    cATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACACAACAtt                                                                                                                                  <  1:241765/150‑2 (MQ=255)
                     aTTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTTCCGGGTTCAATGCAAACTGCACAAAacaaca                                                                                                                                 >  2:1503541/1‑150 (MQ=255)
                         gtttgtGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACTTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                                                                             >  1:317103/1‑150 (MQ=255)
                         gtttgtGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCAATCGCGTGGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                                                                             <  1:1015404/150‑1 (MQ=255)
                            tgtgTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCata                                                                                                                          >  2:929864/1‑150 (MQ=255)
                            tgtgTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTCCACATCTGGCGCAATCCATTCCCGTTGCGAAACGTTAAATTACGTAAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCata                                                                                                                          <  1:157190/150‑1 (MQ=255)
                             gtgtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCatat                                                                                                                         >  2:1255289/1‑150 (MQ=255)
                              tgtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATc                                                                                                                        >  2:1606693/1‑150 (MQ=255)
                               gtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTTCCGGGTGCAATTCATCCGGCACACAACAACACACCATATcc                                                                                                                       >  2:919488/1‑150 (MQ=255)
                               gtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGGCGGGTGCAAATCAAACGGCCCACAACATCACACCatatac                                                                                                                       >  1:1386288/1‑148 (MQ=255)
                               gtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                                                                       >  2:34360/1‑150 (MQ=255)
                               gtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                                                                       <  1:1288703/150‑1 (MQ=255)
                                 aGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGCCCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAg                                                                                                                     >  2:226901/1‑150 (MQ=255)
                                 aGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAg                                                                                                                     >  2:1223602/1‑150 (MQ=255)
                                  gAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgc                                                                                                                    >  1:486560/1‑150 (MQ=255)
                                   aaaCGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                                                                   >  2:1685593/1‑150 (MQ=255)
                                      cGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGTCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATACGGCACACAACATCACACCATATCCAGCGCAg                                                                                                                >  2:1733899/1‑150 (MQ=255)
                                             cGCTTCAGAAACGAAAAATGTTACCAAATAGAAACCTGCACATCTGGCGCAATCAATTCGCGTTGCTAAACTTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGAACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                                                                         <  1:466185/150‑1 (MQ=255)
                                             cGCTTCACAAACGCAAAATGTGACCAAAGAGAACCGTGCACTTCGGGCGCAATCCATTCGCGTTGCTAAACTTTAATTTACGTTAATTTGGTGTGCCGGGTGAAATGCATCCGGAACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                                                                         <  2:1756196/150‑1 (MQ=255)
                                              gCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCTATCCATTCGCCTTGAGAAACGTTAATTTACGTTAATGTTGTGACCCGGGTGCAATGCATCCGGCACACAACATCACACCATAACCAGCCCAGTTTTTcca                                                                                                        >  1:221210/1‑149 (MQ=255)
                                              gCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACCACATCACACCATATCCAGCGCAGGTTTTccg                                                                                                        >  1:1223065/1‑149 (MQ=255)
                                                  cAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATTCAACCCGCCCACAACCTAACACCATATACAGCGCACTTTTTACtattt                                                                                                    >  1:1640155/1‑149 (MQ=255)
                                                  cAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                                                                    <  1:31493/150‑1 (MQ=255)
                                                   aGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTgg                                                                                                   <  2:879742/150‑1 (MQ=255)
                                                    gAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGt                                                                                                  <  2:1239904/150‑1 (MQ=255)
                                                    gAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGt                                                                                                  >  2:829293/1‑150 (MQ=255)
                                                      aaCGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGc                                                                                                >  2:255120/1‑150 (MQ=255)
                                                        cGCAAAAAGTGACCAATGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGAATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCg                                                                                              <  1:453347/150‑1 (MQ=255)
                                                        cGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCCCCCAACATCCCCCCCTATACCGCGCAGTTTTTCGTTTTGGTGCCg                                                                                              >  2:465621/1‑150 (MQ=255)
                                                        cGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCg                                                                                              >  2:155349/1‑150 (MQ=255)
                                                         gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCgg                                                                                             >  2:1029780/1‑150 (MQ=255)
                                                         gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCAGGCACACAAAATCACACCATATCCAGCGGAGTTTTTCCATTTGGTGCCgg                                                                                             >  1:1345341/1‑150 (MQ=255)
                                                         gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATCCATCCGGCAGACAACATCACACCATATCCCGCGCAGTTTTTCCTTTTGGTGCCgg                                                                                             >  2:1724287/1‑150 (MQ=255)
                                                          cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCAAGAGAAGTTTTTCCTTTTGGTGCCGGa                                                                                            >  1:156002/1‑150 (MQ=255)
                                                          cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTGGTGTGCCGGGGGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                                                                            <  1:249472/150‑1 (MQ=255)
                                                          cAAAAAGAGCCCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGATGCGAAACGTTAATATACGTTAATGTTGTGTGCCGGGTTCAAAGCAACCCGCACACAACATAACACCATAACCACCCCAGTTTTTACTTTTGCTGCCCGa                                                                                            >  1:378131/1‑150 (MQ=255)
                                                           aaaaaGTGACCAAAGAGAAACCTGCACATCTGCCGCAATCCATTCTCGTTGCGAAACGTTAATTTACGTTAATGGTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGat                                                                                           <  1:121256/150‑1 (MQ=255)
                                                               aGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGc                                                                                       >  2:40154/1‑150 (MQ=255)
                                                                    ccAAAGAGAAACCTGCACATCTTGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATACGGCACACAACATCACACCATATCCAGCGGAGTTTTTGATTTTGGTGCCGGATATGCCTTAt                                                                                  >  1:113080/1‑150 (MQ=255)
                                                                      aaaGAGATACCTGCACATCTGGCGCAATCCATTGGCGTTGCGAATCGTTAATTTTCGTTAATGTTGCGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTATGCTTTTGGTGCCGGATATGCCTTATcc                                                                                >  1:1033821/1‑150 (MQ=255)
                                                                      aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                                                                                <  2:174812/150‑1 (MQ=255)
                                                                      aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTAACGTAAATGTTTGGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                                                                                <  1:1607376/150‑1 (MQ=255)
                                                                          agaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCa                                                                            <  2:1227617/150‑1 (MQ=255)
                                                                          agaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGACGGGTGCAATTCATCCGGCACACAACATAACCCCATATACAGCGCAGTTTTTACTTTTGGTGCCGGATATTCCTTTTACAGCa                                                                            >  2:1347572/1‑150 (MQ=255)
                                                                              aCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACTTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGc                                                                        <  1:616286/150‑1 (MQ=255)
                                                                               ccTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCt                                                                       <  2:1598400/150‑1 (MQ=255)
                                                                               ccTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTACTTTTGGTGCCGGATATGCCTTATACAGCATAGCt                                                                       >  2:1826506/1‑150 (MQ=255)
                                                                                  gCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTTCCCGGTGCAATGCAACCGGCACACAACATCACACCAAATACCACGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                                                                    >  2:556223/1‑150 (MQ=255)
                                                                                      aTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCg                                                                >  1:18283/1‑150 (MQ=255)
                                                                                       cctgccccAAGCCATCCCCGGTGCGAAACGGAAATTTACGTAAATCGGTCGTGCCCGGTCCAAGGCAACCGGCACACAACATCACACCATATCCACCGCAGTTTTTCCATATGCTGCCGGATAAGCCTTATCCAGCATAGCTAATACCGc                                                               <  2:683603/143‑1 (MQ=255)
                                                                                        gtGGCGCAATCCATTCGCGTTGCGAAACTTTAATTTACGTTAATGTTGTGTGCCGGGTGAAATGAATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt                                                              <  1:1663265/149‑1 (MQ=255)
                                                                                          ggCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTTCAGGGTGCAATGCATCCGGCACACAACATCCCACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCgagga                                                            >  2:1739839/1‑146 (MQ=255)
                                                                                          ggCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTga                                                            <  1:713647/150‑1 (MQ=255)
                                                                                           gctcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTGGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa                                                           <  1:890781/147‑1 (MQ=255)
                                                                                           gcgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa                                                           >  2:1140324/1‑150 (MQ=255)
                                                                                           gcgaAATCCATTCGCGTTGCTAACCTTTATTTTACGTAATTGTGGTGGGCCGGGTCAAATGCAGCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa                                                           <  1:1370672/150‑1 (MQ=255)
                                                                                            cgaaaTCCATTCGGGTGGCGAAAGGTTAATTTACGTTAATGTTGTGTGCCGGGTGAAATGCATCCGGAACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag                                                          <  2:803920/147‑1 (MQ=255)
                                                                                             gcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTTGTGCCGGATATGCCTTATACAGCATAGCTAATTCCGCTgaaga                                                         >  2:1662915/1‑150 (MQ=255)
                                                                                               aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGaaa                                                       >  2:1682370/1‑150 (MQ=255)
                                                                                                   cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTTCCGGGTGCAAAGCAAACGGCACACAACAAAACAACAAATACACCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATACAGCATAGCTAAATACGCTGAAGCAAGTTc                                                   >  2:1318359/1‑150 (MQ=255)
                                                                                                    aTTCGCGTTGCTAAACTTTAATTTAGTTTATTGTTGTGTGCCGGTTTCAATGCATCCGGAACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCg                                                  <  2:1459081/150‑1 (MQ=255)
                                                                                                     ttCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGa                                                 <  1:34360/150‑1 (MQ=255)
                                                                                                          gTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCa                                            >  1:1783087/1‑150 (MQ=255)
                                                                                                                  tggTTAATTAACGTTATTGTTGTGTGCCGGGTGAAATTCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGc                                    <  1:919488/148‑1 (MQ=255)
                                                                                                                   cGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGcc                                   >  1:698013/1‑150 (MQ=255)
                                                                                                                              gTTAATGTTGTGTGCCGGGTGCAATACATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTttgttg                        <  2:1722021/150‑1 (MQ=255)
                                                                                                                               ttAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTgg                       <  1:117351/150‑1 (MQ=255)
                                                                                                                                      tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACCGCCGCATTTTGTTGGACATggg                >  2:149571/1‑150 (MQ=255)
                                                                                                                                      tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATggg                <  1:963214/150‑1 (MQ=255)
                                                                                                                                      tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATggg                <  1:477308/150‑1 (MQ=255)
                                                                                                                                        agtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTGTATCCTTTTAGTGCCTGATATGACTTATCCAGCATAGCTAAATCCGCTGAAGCAAGTTCGCCCTCACGCACCGCCGCATATTGTTGGACATGGGCa              >  2:1302753/2‑150 (MQ=255)
                                                                                                                                         gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCaa             >  1:1046417/1‑150 (MQ=255)
                                                                                                                                         gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCaa             >  1:1400070/1‑150 (MQ=255)
                                                                                                                                         gtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCaa             >  2:130765/1‑150 (MQ=255)
                                                                                                                                              gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTg        <  2:19056/150‑1 (MQ=255)
                                                                                                                                               ggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTgg       >  1:1339741/1‑150 (MQ=255)
                                                                                                                                                gTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGc      <  2:465906/150‑1 (MQ=255)
                                                                                                                                                gTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGc      <  1:1682370/150‑1 (MQ=255)
                                                                                                                                                gTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGc      <  1:21776/150‑1 (MQ=255)
                                                                                                                                                gTGCAATGCATCCGGAACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGc      <  1:1347572/150‑1 (MQ=255)
                                                                                                                                                   cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGc   >  2:438029/1‑150 (MQ=255)
                                                                                                                                                   cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGc   >  1:412925/1‑150 (MQ=255)
                                                                                                                                                    aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGCt  >  2:1447268/1‑150 (MQ=255)
                                                                                                                                                    aaTGCATCCGGAACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGACGGATATGCCTTATCAAGCATAGCTAATTCCGCTGAAGAAAGTTCGACATCAAGAACAGCAGCATTTTGATGGACATGGGCAATCGTGGACGCt  >  1:874869/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGCT  >  NC_000913/1864614‑1864911

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