Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,435 G→A 7.1% intergenic (+17/+18) ydfZ → / ← ydfI selenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,629,4350GA 7.1% 525.5 / 18.1 168intergenic (+17/+18)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (88/68);  new base A (5/7);  total (93/75)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 171 total aligned reads displayed.

ACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCA  >  NC_000913/1629287‑1629583
                                                                                                                                                    |                                                                                                                                                    
acTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGTAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGa                                                                                                                                                     <  1:1099124/150‑1 (MQ=255)
  tGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATg                                                                                                                                                   >  2:544593/1‑150 (MQ=255)
     cAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGTCACCACTTGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcg                                                                                                                                                >  2:2044828/1‑150 (MQ=255)
       aGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcggc                                                                                                                                              >  2:2126805/1‑150 (MQ=255)
         aTCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCTTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAGGGCGGCAt                                                                                                                                            >  1:1332671/1‑150 (MQ=255)
          tCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                                           >  1:31050/1‑150 (MQ=255)
          accTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTTGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAtt                                                                                                                                           <  1:462617/149‑1 (MQ=255)
            cTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGc                                                                                                                                         >  1:590339/1‑150 (MQ=255)
             tGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcac                                                                                                                                        <  1:1244609/150‑4 (MQ=255)
               tCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttt                                                                                                                                      <  2:489562/150‑1 (MQ=255)
                cGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcacacg                                                                                                                                     <  1:1304311/150‑7 (MQ=255)
                   ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                  >  1:11088/1‑150 (MQ=255)
                      aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGGGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                               <  2:26842/150‑1 (MQ=255)
                      aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGAATGAACTAAGCTTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                               >  2:709604/1‑150 (MQ=255)
                          tGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTGAAGGTGGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACgg                                                                                                                          <  1:1502038/150‑1 (MQ=255)
                             aGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGaa                                                                                                                        >  2:1708333/1‑150 (MQ=255)
                                tCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGGTTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAatat                                                                                                                     <  2:2107459/150‑1 (MQ=255)
                                tCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAatat                                                                                                                     >  1:998222/1‑150 (MQ=255)
                                 cTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAtt                                                                                                                    >  2:1994242/1‑150 (MQ=255)
                                     ccGCTGAGCAAATCCTCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                <  2:1845500/150‑1 (MQ=255)
                                     ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                >  1:10304/1‑150 (MQ=255)
                                        cagAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGGTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCAc                                                                                                             <  1:373520/148‑1 (MQ=255)
                                          gAGCAAATCCGCCGCTTAAAAACCGTAGTTGTTTTCGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACgg                                                                                                           <  2:293702/150‑1 (MQ=255)
                                            gCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGTCACCACTGGCCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCAGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTc                                                                                                         >  1:114491/1‑150 (MQ=255)
                                              aaaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGAATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTAACGGAATATTTTCCCAAGTTCGc                                                                                                       >  2:1250265/1‑150 (MQ=255)
                                                aTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGGTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCtt                                                                                                     <  1:2036210/150‑1 (MQ=255)
                                                   cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                  >  2:819681/1‑150 (MQ=255)
                                                   cgccgcGGAAAAACCGGAGTTGTTGAAGGTTGTGAAGAGAAACTGTCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                  <  1:2144429/150‑1 (MQ=255)
                                                     ccgcGGAAAAACCGTAGTTGTTGCAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                <  2:492219/150‑1 (MQ=255)
                                                     ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                >  2:732057/1‑150 (MQ=255)
                                                       gcgGAAACACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCCCCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGcc                                                                                              >  1:1464290/1‑150 (MQ=255)
                                                       gcgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGcc                                                                                              >  1:889781/1‑150 (MQ=255)
                                                         ggAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                >  2:924676/1‑146 (MQ=255)
                                                          gAAAACCAGTAGTTGTTGAAGGGTGTGACGAGAAACTGGCCCCCCTGGACATGATTCGGCTCGGCATGCCCTAAGCGTGGGACTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATg                                                                                           <  1:1353971/150‑1 (MQ=255)
                                                          gAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCAATGGACCTGATTCGTCTCGGAATGAACTAAGCCTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTAACGGAATATTTTGCCACGGTCGCTTTCTCGCCATg                                                                                           >  2:534918/1‑150 (MQ=255)
                                                           aaaaaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATgc                                                                                          >  1:1039753/1‑150 (MQ=255)
                                                           aaaaaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATgc                                                                                          <  2:19213/150‑1 (MQ=255)
                                                              aaCCGTAGTTGTTGAAGGTTGTGAAGAGACACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCt                                                                                       <  1:1071910/150‑1 (MQ=255)
                                                              aaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCt                                                                                       <  1:819849/150‑1 (MQ=255)
                                                               aCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCta                                                                                      <  1:645356/150‑1 (MQ=255)
                                                                ccGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCACGCCATGCGCtaa                                                                                     <  2:667959/150‑1 (MQ=255)
                                                                        gttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACa                                                                             >  1:638407/1‑150 (MQ=255)
                                                                          tgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                           <  2:307458/150‑1 (MQ=255)
                                                                           gAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGt                                                                          >  2:1490568/1‑150 (MQ=255)
                                                                           gAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGt                                                                          >  2:281745/1‑150 (MQ=255)
                                                                             aGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTAc                                                                        >  1:396623/1‑150 (MQ=255)
                                                                                  gtgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTc                                                                   <  2:21673/150‑1 (MQ=255)
                                                                                   tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTGCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATATAGACAAGTACGTTTcc                                                                  >  2:1171759/1‑150 (MQ=255)
                                                                                     aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAAGAAAGACAAGTACGTTTCCGt                                                                >  2:1495174/1‑150 (MQ=255)
                                                                                        agaATCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAc                                                             <  1:799240/150‑1 (MQ=255)
                                                                                         gaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATg                                                                                           >  1:294492/1‑119 (MQ=255)
                                                                                         gaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCCct                                                            >  2:1695712/1‑150 (MQ=255)
                                                                                           aaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCCTGCGCTAATAAAGACAAGTACGTTTCCGTCActct                                                          >  1:328447/1‑150 (MQ=255)
                                                                                           aaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAatct                                                          >  2:1902752/1‑150 (MQ=255)
                                                                                           aaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCAACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCActct                                                          >  2:1805280/1‑150 (MQ=255)
                                                                                                gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAg                                                     >  1:742951/1‑150 (MQ=255)
                                                                                                gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGGGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGGCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAg                                                     <  2:409674/150‑1 (MQ=255)
                                                                                                  accacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTa                                                   <  2:1661878/150‑1 (MQ=255)
                                                                                                    cacTGGACCTGATTCGTCTCGGCCTGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                 <  1:941420/150‑1 (MQ=255)
                                                                                                    cacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                 >  2:2143815/1‑150 (MQ=255)
                                                                                                    cacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATACCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTAGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                 >  1:2155486/1‑150 (MQ=255)
                                                                                                           ccTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTa                                          >  1:2079681/1‑150 (MQ=255)
                                                                                                           ccTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCAATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTa                                          >  2:834347/1‑150 (MQ=255)
                                                                                                            cTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAt                                         <  2:769287/150‑1 (MQ=255)
                                                                                                              gATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTg                                       <  1:1502323/150‑1 (MQ=255)
                                                                                                              gATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTg                                       >  1:914839/1‑150 (MQ=255)
                                                                                                              gATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTg                                       >  1:5312/1‑150 (MQ=255)
                                                                                                                  cGTCTCGGCATGAACTAAGCTTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTg                                   <  2:755538/150‑1 (MQ=255)
                                                                                                                  cGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAACTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTg                                   >  2:1839380/1‑150 (MQ=255)
                                                                                                                   gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                  <  2:316324/150‑1 (MQ=255)
                                                                                                                   gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                  <  1:1994242/150‑1 (MQ=255)
                                                                                                                    tctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGc                                 <  2:1742557/150‑1 (MQ=255)
                                                                                                                    tctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGc                                 <  2:138443/150‑1 (MQ=255)
                                                                                                                    tctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGc                                 <  1:709604/150‑1 (MQ=255)
                                                                                                                    tctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGc                                 <  2:1089283/150‑1 (MQ=255)
                                                                                                                     ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                <  1:1826616/150‑1 (MQ=255)
                                                                                                                     ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                <  2:1386816/150‑1 (MQ=255)
                                                                                                                     ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                >  1:1909430/1‑150 (MQ=255)
                                                                                                                      tcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCTGTCACTCTTGCAGTAAACAACCTATTGTCTGGCaa                               >  2:320440/1‑150 (MQ=255)
                                                                                                                      tcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCaa                               >  1:1756053/1‑150 (MQ=255)
                                                                                                                       cGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCaaa                              >  2:1281933/1‑150 (MQ=255)
                                                                                                                          cATGAACTAAGCGTGTGAATGTCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTCTCTGGCAAatca                           >  2:1198084/1‑150 (MQ=255)
                                                                                                                                   aGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAg                  <  2:1148263/150‑1 (MQ=255)
                                                                                                                                       tgtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGATTATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGCCAAGTACGGTTCCGTCACTGTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACAAAAGATCg              >  1:503680/1‑150 (MQ=255)
                                                                                                                                       tgtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCg              >  2:342686/1‑150 (MQ=255)
                                                                                                                                       tgtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCg              >  2:506547/1‑150 (MQ=255)
                                                                                                                                       aagcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCg              >  1:1542946/5‑150 (MQ=255)
                                                                                                                                        agcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc             >  2:802369/4‑150 (MQ=255)
                                                                                                                                         tgtaTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc            <  2:1111243/147‑1 (MQ=255)
                                                                                                                                         tgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc            <  1:257976/150‑1 (MQ=255)
                                                                                                                                         tgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc            >  1:193453/1‑150 (MQ=255)
                                                                                                                                         gcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc            >  1:261491/3‑150 (MQ=255)
                                                                                                                                          gAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCt           >  2:2162311/1‑150 (MQ=255)
                                                                                                                                          caaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCt           >  1:1199555/2‑150 (MQ=255)
                                                                                                                                             tGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGAGCGCCTTaa        >  2:1460458/1‑150 (MQ=255)
                                                                                                                                             tGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACGCTTGCAGTAAACAAACTATTGTCTGGCAAATCCTCACCAAAGATCGCGTTaa        >  2:198019/1‑150 (MQ=255)
                                                                                                                                              gccgccGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTTTTGTCTGGCAAATCATCACCAAAGATCGCCTTAAt       <  1:1380474/150‑1 (MQ=255)
                                                                                                                                                 gccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGGCTTGCGCGAATAAAGACAAGTACGGTTCCGTCACTCTTGCAGGAAACAAACTATTGGCTGGCGAATCATCACCAAAGATCGCCTTAATCGc    >  2:229392/1‑150 (MQ=255)
                                                                                                                                                  ccgATGGCGGCATTGCTTTTTTACTTCACTGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGcc   <  1:779357/150‑1 (MQ=255)
                                                                                                                                                   cgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCa  <  1:2047689/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCA  >  NC_000913/1629287‑1629583

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: