Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,024,038 G→C 5.1% G908A (GGC→GCC)  xdhD → putative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,024,0380GC 5.1% 355.2 / 8.1 98G908A (GGC→GCC) xdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (39/54);  new base C (2/3);  total (41/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

AAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGC  >  NC_000913/3023795‑3024275
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
aaaTCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACTCCGGTCAATCCAGAGTTAGCGTTGGGAAAAATCTACGGTGCCACCCTGCGTGATATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTAACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGTGAAATCGATCTCGGAAATc                                                                                                                                                                                                                <  2:119635/275‑1 (MQ=255)
  aTCCTCCTGGTTAAATTCTACGCCATGCTGGACTGCGGTACACCGGTTAATCCAGATTTAGGTTTGGGACAAAGCGACGGTGCCACCCTGGGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTGATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGGGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATc                                                                                                                                                                                                                <  4:212588/273‑1 (MQ=255)
         tGGATAAATTCTACGCCTTGCTGGACTGGGGTCAACCGGTAAATCCCGAGTTAGCGTTGGGACAAATCTACGGGTCCACCCTGCGTGCTATCGGCAACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTAACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTaaa                                                                                                                                                                                                        <  4:233597/274‑1 (MQ=255)
             tAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGgcg                                                                                                                                                                                                   <  2:256895/275‑1 (MQ=255)
              aaaTTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGgcg                                                                                                                                                                                                   <  2:81051/274‑1 (MQ=255)
                     aCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGACCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCg                                                                                                                                                                                           >  2:221074/1‑275 (MQ=255)
                                                       agagTTAGCGTTGGGGCAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATACTTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGc                                                                                                                                                          <  4:143215/274‑1 (MQ=255)
                                                                       aaaatcgagcgggCCACCCTGGGGGGGTTCGGCCAAAGTATGAGGGGAGGGATCAATTATTAAGCGTAAGGTCACNCGTTAACGCGTTTTTTACGCAGTTACGGCGCACCGAAAGTTGGTGACATTCCGCGTGAGTTCCGCGCTGGGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGaa                                                                                                                                           <  3:132837/261‑1 (MQ=255)
                                                                        aaaTCTACGGGGCCACCCTGCGTGCTATCGGCCACAGTATGGGCGAAGAGATCATTTATGACGCCGAAGGTCTCCCGTTAACGCGGGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGAGACGATAAAGTCGGCCCGTTCGGGGCTAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATg                                                                                                                                         <  2:220311/274‑1 (MQ=255)
                                                                                    agaCCCGGCGGGCTATCGGCCACAGTATGGGGAAAGGGATCATTGAGGGGGGCGAAGGTCACGAGTTCACGCGGGAGTTACGCAGTTAGGGCGAACCGAAAATAGGGGAAATGCCGGGTGATTTCCGCGCTGTGCTGGTGCCGAGCGAAGATAAAGTCGGCCCGTTCGGGGTGAAATCGATCTCGGAAAGCGGTGTAAATGGAGCAGGTCCGGCGATTGCTACCGAAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACAcc                                                                                                                             <  4:329634/272‑1 (MQ=255)
                                                                                     cgccacgCGTGCTGTGCGCCACAGTTAGGNNNAAGAGCACATTGCATACGCTGAGGGTGGCCCGTTGAGGCTTAAGTTACGCCGTTACGGCGGCCCGAACATTGCTGACAGTCCGCGTGATTTCCGCGTTGTGCTGGTGCCGGGCGACGATAAAGTCGGCCCGTTCGGCGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTTCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCgg                                                                                                                           <  3:146328/269‑1 (MQ=255)
                                                                                        ccggCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCgg                                                                                                                                                                                          <  2:268146/206‑1 (MQ=255)
                                                                                            cGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGGAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTca                                                                                                                                 >  2:72241/1‑262 (MQ=255)
                                                                                                cttcaggCCACAGTATGGGGGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTAGGCAGTTAGGGCGCACCGAAAATTGGGGACATTCGGCGTGATTTCCGCGCTGTGCTGGTGCCGGGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGGAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACAcc                                                                                                                             <  4:72240/257‑1 (MQ=255)
                                                                                                  agcGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCa                                                                                                               <  2:450477/272‑1 (MQ=255)
                                                                                                          caGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                                                                                                                                                                                                                                       <  5:589149/146‑1 (MQ=255)
                                                                                                          caGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                                                                                                                                                                                                                                       >  6:589149/1‑146 (MQ=255)
                                                                                                           aGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGgcgc                                                                                                                                                                                                  >  1:139502/1‑182 (MQ=255)
                                                                                                               tGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCt                                                                                                                                                                                              <  3:139502/182‑1 (MQ=255)
                                                                                                                         gttCATTTATGCCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGgag                                                                                                                         <  2:251063/239‑1 (MQ=255)
                                                                                                                          aTCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg                                                                                                                                                                                                                     >  7:364091/1‑148 (MQ=255)
                                                                                                                           tCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCg                                                                                                                                                                                                                                  <  7:29654/134‑1 (MQ=255)
                                                                                                                            cATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                                                                                                                                                                                                                                       >  7:595519/1‑128 (MQ=255)
                                                                                                                            cATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                                                                                                                                                                                                                                       <  8:595519/128‑1 (MQ=255)
                                                                                                                              tttATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt                                                                                                                                                                                                                 >  5:84682/1‑148 (MQ=255)
                                                                                                                                 aTGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                                                                                                                              <  8:364091/148‑1 (MQ=255)
                                                                                                                                  tGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCg                                                                                                                                                                                                               >  5:197829/1‑146 (MQ=255)
                                                                                                                                  tGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCg                                                                                                                                                                                                               <  6:197829/146‑1 (MQ=255)
                                                                                                                                    aCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                                                                                                                 >  8:1110640/1‑126 (MQ=255)
                                                                                                                                    aCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                                                                                                                 <  7:1110640/126‑1 (MQ=255)
                                                                                                                                       ccGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTgcgc                                                                                                         >  4:114728/1‑243 (MQ=255)
                                                                                                                                        cGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGgtg                                                                                                                                                                                                            >  7:326912/1‑143 (MQ=255)
                                                                                                                                        cGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGgtg                                                                                                                                                                                                            <  8:326912/143‑1 (MQ=255)
                                                                                                                                          aaGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATgg                                                                                                                                                                                                     >  7:994681/1‑148 (MQ=255)
                                                                                                                                           aGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGa                                                                                                     <  2:114729/243‑1 (MQ=255)
                                                                                                                                               cACGCGTTAACGGGTGATTTACGCCGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc                                                                                                                                                                                                                       >  5:1015935/1‑125 (MQ=255)
                                                                                                                                               cACGCGTTAACGGGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc                                                                                                                                                                                                                       <  6:1015935/125‑1 (MQ=255)
                                                                                                                                               cACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg                                                                                                                                                                                                                            >  5:8274/1‑120 (MQ=255)
                                                                                                                                               cACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg                                                                                                                                                                                                                            <  6:8274/120‑1 (MQ=255)
                                                                                                                                                 cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg                                                                                                                                                                                                                     >  8:166550/1‑125 (MQ=255)
                                                                                                                                                 cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg                                                                                                                                                                                                                     <  7:166550/125‑1 (MQ=255)
                                                                                                                                                 cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt                                                                                                                                                                                                                 >  8:312867/1‑129 (MQ=255)
                                                                                                                                                 cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt                                                                                                                                                                                                                 <  7:312867/129‑1 (MQ=255)
                                                                                                                                                     ttAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCgg                                                                                                                                                                                          >  8:543087/1‑148 (MQ=255)
                                                                                                                                                        aCGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGa                                                                                                                                                                                       <  6:84682/148‑1 (MQ=255)
                                                                                                                                                         cgcgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGAt                                                                                                                                                                                      <  8:994681/148‑1 (MQ=255)
                                                                                                                                                           cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                                                                                                                                                                                                                               >  6:620603/1‑105 (MQ=255)
                                                                                                                                                           cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                                                                                                                                                                                                                               <  5:620603/105‑1 (MQ=255)
                                                                                                                                                           cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTg                                                                                                                                                                                    >  8:839423/1‑148 (MQ=255)
                                                                                                                                                           cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                                                                                                                 <  7:80054/103‑1 (MQ=255)
                                                                                                                                                           cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                                                                                                                 >  8:80054/1‑103 (MQ=255)
                                                                                                                                                             tGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGc                                                                                                                                                                                   >  7:1120779/1‑147 (MQ=255)
                                                                                                                                                             tGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGc                                                                                                                                                                                   <  8:1120779/147‑1 (MQ=255)
                                                                                                                                                                tttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATg                                                                                                                                                                                                      >  8:418478/1‑125 (MQ=255)
                                                                                                                                                                tttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATg                                                                                                                                                                                                      <  7:418478/125‑1 (MQ=255)
                                                                                                                                                                 ttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCgg                                                                                                                                                                                                                                      <  6:78812/92‑1 (MQ=255)
                                                                                                                                                                 ttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCgg                                                                                                                                                                                                                                      >  5:78812/1‑92 (MQ=255)
                                                                                                                                                                 ttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGa                                               <  4:221073/275‑1 (MQ=255)
                                                                                                                                                                 ttACGCAGTGACGGCGAACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGa                                               <  1:115781/275‑1 (MQ=255)
                                                                                                                                                                         ttACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCAc                                                                                                                                                                      >  8:83105/1‑148 (MQ=255)
                                                                                                                                                                          tACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAc                                                                                                                                                                                <  7:152941/137‑1 (MQ=255)
                                                                                                                                                                          tACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAc                                                                                                                                                                                >  8:152941/1‑137 (MQ=255)
                                                                                                                                                                                 cACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg                                                                                                                                                                                                                     >  4:305933/1‑93 (MQ=255)
                                                                                                                                                                                 cACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGatgcatg                                                                                                                                                              >  7:908940/1‑148 (MQ=255)
                                                                                                                                                                                     gAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt                                                                                                                                                                                                                 <  2:305934/93‑1 (MQ=255)
                                                                                                                                                                                      aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                                                                                                                              <  7:78226/95‑1 (MQ=255)
                                                                                                                                                                                      aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                                                                                                                              <  5:359205/95‑1 (MQ=255)
                                                                                                                                                                                      aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                                                                                                                              >  8:78226/1‑95 (MQ=255)
                                                                                                                                                                                      aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                                                                                                                              >  6:359205/1‑95 (MQ=255)
                                                                                                                                                                                      aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAcc                                                                                                                                                                               >  7:22369/1‑126 (MQ=255)
                                                                                                                                                                                      aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAcc                                                                                                                                                                               <  8:22369/126‑1 (MQ=255)
                                                                                                                                                                                      aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCa                                                                                                                                                         <  7:543087/148‑1 (MQ=255)
                                                                                                                                                                                          ttGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATCCATGCGGCATCTg                                                                                                                                                     <  7:83105/148‑1 (MQ=255)
                                                                                                                                                                                           tGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTgg                                                                                                                                                    <  7:839423/148‑1 (MQ=255)
                                                                                                                                                                                            ggTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGt                                                                                                                                                   <  8:908940/148‑1 (MQ=255)
                                                                                                                                                                                                  aTTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTca                                                                                                                                 >  2:244421/1‑160 (MQ=255)
                                                                                                                                                                                                   ttCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                                                                                                                     <  8:351292/123‑1 (MQ=255)
                                                                                                                                                                                                   ttCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                                                                                                                     >  7:351292/1‑123 (MQ=255)
                                                                                                                                                                                                      cgcgTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACAcc                                                                                                                             <  4:244420/160‑1 (MQ=255)
                                                                                                                                                                                                         gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc                                                                                                                                                                                                                         <  8:1090877/65‑1 (MQ=255)
                                                                                                                                                                                                         gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc                                                                                                                                                                                                                         >  7:1090877/1‑65 (MQ=255)
                                                                                                                                                                                                         gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTg                                                                                                                                                                                    <  5:1002227/102‑1 (MQ=255)
                                                                                                                                                                                                         gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTg                                                                                                                                                                                    >  6:1002227/1‑102 (MQ=255)
                                                                                                                                                                                                               tCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGc  <  2:6079/274‑1 (MQ=255)
                                                                                                                                                                                                                           tGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAAt                                                                                                                    >  6:16565/1‑148 (MQ=255)
                                                                                                                                                                                                                              tGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTgg                                                                                                                                                    >  5:187459/1‑113 (MQ=255)
                                                                                                                                                                                                                              tGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTgg                                                                                                                                                    <  6:187459/113‑1 (MQ=255)
                                                                                                                                                                                                                               gCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGatgcatg                                                                                                                                                              >  7:70721/1‑102 (MQ=255)
                                                                                                                                                                                                                               gCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTCCTACCGCAATTCACGatgcatg                                                                                                                                                              <  8:70721/102‑1 (MQ=255)
                                                                                                                                                                                                                                  gAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTcc                                                                                                                                                                                            >  3:29039/1‑69 (MQ=255)
                                                                                                                                                                                                                                   aGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTg                                                                                                            <  5:16565/148‑1 (MQ=255)
                                                                                                                                                                                                                                   aGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTg                                                                                                            >  8:590988/1‑148 (MQ=255)
                                                                                                                                                                                                                                     cgacgaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTgcg                                                                                                          <  7:590988/148‑1 (MQ=255)
                                                                                                                                                                                                                                      gacgaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCg                                                                                                                                                                                        <  1:29039/69‑1 (MQ=255)
                                                                                                                                                                                                                                         gaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATGTCACACCGGAGAAAATACTCACTgcgc                                                                                                         >  7:866495/1‑145 (MQ=255)
                                                                                                                                                                                                                                         gaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATGTCACACCGGAGAAAATACTCACTgcgc                                                                                                         <  8:866495/145‑1 (MQ=255)
                                                                                                                                                                                                                                            aaaGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTg                                                                                                       >  5:1043576/1‑144 (MQ=255)
                                                                                                                                                                                                                                            aaaGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTg                                                                                                       <  6:1043576/144‑1 (MQ=255)
                                                                                                                                                                                                                                                   gCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGaaaaa                                                                                                 >  5:295580/1‑143 (MQ=255)
                                                                                                                                                                                                                                                   gCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGaaaaa                                                                                                 <  6:295580/143‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
AAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGC  >  NC_000913/3023795‑3024275

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: