Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,591 T→G 5.4% intergenic (+73/‑32) yqhD → / → dkgA aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5910TG 5.4% 788.6 / 29.9 223intergenic (+73/‑32)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base T (100/111);  new base G (5/7);  total (105/118)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

Only 100 of 224 total aligned reads displayed.

GTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAG  >  NC_000913/3156319‑3156850
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                   
gaaTTGACGCGGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCTCCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCGCCGGAGCCTg                                                                                                                                                                                                                                                                    <  2:145338/272‑1 (MQ=255)
                        cccGGAATTTCTTGGAGCAATGAGGGGGGCCGACCCACCTGTCCGACTACTGGCTGGACGGCAGCTCCATCCCGGCTTTGATGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                                                                                                                                                                                            <  2:321102/274‑1 (MQ=255)
                              aTTTCTTTGGGCAATTAGGCGGGCCAGCCCACCTCTCCGACTACGGGTTGGACGGCAGCGCCATCCCGGCTTTGCTGAAAAAACTTGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGt                                                                                                                                                                                                                                      <  4:152120/274‑1 (MQ=255)
                                                 cgccccccccccccTCCCCGGCNNCGTTCNGAACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTTGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                                                                                                                                                                                  <  3:16933/267‑1 (MQ=255)
                                                 cGTGCCGACCCACCTCTCCGAGTACGGTCTGGACGGGAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                                                                                                                                                                                  <  2:60736/275‑1 (MQ=255)
                                                       gACCCTACTCTCCGACGACGGTCTGGACGGCAGCTCCATCCCGGTGTTGCTGAAAAAACTGGAAGGGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCATCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                                                                                                                                                                                             <  4:80550/274‑1 (MQ=255)
                                                        gcccACCTGTCCGACTACGGTCTGGACGGCAGCTCCAGCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTa                                                                                                                                                                                                           <  2:208245/274‑1 (MQ=255)
                                                                  cccaataggctgtGGACGGCAGCTCCACNCNGGGTTTGCTGAAAAAACGGGAAGAGCACGGCATGACCCAACTGGGNGAAAATAATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCNCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                                                                                                                                 <  3:32472/263‑1 (MQ=255)
                                                                                gACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTAGGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACATCACCGGAACCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCCATGTCCTGCCCCa                                                                                                                                                                                    >  4:302946/1‑274 (MQ=255)
                                                                                                                                aTGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTc                                                                                                                                                                                                                                                                 >  5:273527/1‑149 (MQ=255)
                                                                                                                                   aCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCgg                                                                                                                                                                                                                                                              <  7:776323/149‑1 (MQ=255)
                                                                                                                                    cccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                             <  5:1142750/149‑1 (MQ=255)
                                                                                                                                     ccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                                                                                                                                                                                            >  7:1437633/1‑149 (MQ=255)
                                                                                                                                     ccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                                                                                                                                                                                            >  6:1455747/1‑149 (MQ=255)
                                                                                                                                      cAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGa                                                                                                                                                                                                                                                           <  8:689166/149‑1 (MQ=255)
                                                                                                                                      cAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                                                                                                                                                                                           <  5:1238228/149‑1 (MQ=255)
                                                                                                                                       aaCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                                                                                                                                                                                          <  8:632904/149‑1 (MQ=255)
                                                                                                                                        aCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                         >  7:1340632/1‑149 (MQ=255)
                                                                                                                                         cTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGtt                                                                                                                                                                                                                                                        <  5:990580/149‑1 (MQ=255)
                                                                                                                                            ggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGTGAggtgtt                                                                                                                                                                                                                                                     <  8:1415693/149‑6 (MQ=255)
                                                                                                                                             gCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                                                                                                                                                                                    <  5:765633/149‑1 (MQ=255)
                                                                                                                                               gAAAATCATGAGATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                                                                                                                                                                                  >  7:1145169/1‑149 (MQ=255)
                                                                                                                                               gAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                                                                                                                                                                                  >  6:109068/1‑149 (MQ=255)
                                                                                                                                                 aaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                                                                                                                                                                                                <  8:553279/149‑1 (MQ=255)
                                                                                                                                                 aaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                                                                                                                                                                                                >  7:1163369/1‑149 (MQ=255)
                                                                                                                                                    tCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                                                                                                                                             >  6:1137946/1‑149 (MQ=255)
                                                                                                                                                     cATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                                                                                                                                                                                            <  5:570724/149‑1 (MQ=255)
                                                                                                                                                      aTGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                                                                                                                                                                                           <  8:1163369/149‑1 (MQ=255)
                                                                                                                                                         aCATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                                                                                                                                                                                        <  5:1360370/149‑1 (MQ=255)
                                                                                                                                                           aTTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCCCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                             >  8:1057916/1‑126 (MQ=255)
                                                                                                                                                               cGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                                                                                                                                                                                  <  6:224464/149‑1 (MQ=255)
                                                                                                                                                                gTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTTAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                                                                                                                                                                                 >  8:576135/1‑149 (MQ=255)
                                                                                                                                                                  tGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                                                                                                                                                                                               >  7:348790/1‑149 (MQ=255)
                                                                                                                                                                        tCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                                                                                                                                                                                    <  2:3007/122‑1 (MQ=255)
                                                                                                                                                                        tCAGCCGCCGTATATACGAAGCCGCCAGCTAAGCTTTTTACGCCTCAAACTTTCGTATTCGGGCATTTCGTCCAAACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                                                                                                                                                                                         >  7:904660/1‑149 (MQ=255)
                                                                                                                                                                          aGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                                                                                                                                                                                       <  5:169937/149‑1 (MQ=255)
                                                                                                                                                                           gccgccGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCg                                                                                                                                                                                                                      >  6:621989/1‑149 (MQ=255)
                                                                                                                                                                              tccgTATATACCACGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAATTTCACAACACCTCACCGGAGCCTGCTCCGTTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATAccgccg                                                                                     <  1:37573/274‑1 (MQ=255)
                                                                                                                                                                                  tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                                                                                                                                                                                               >  5:1151413/1‑149 (MQ=255)
                                                                                                                                                                                        cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACa                                                                                                                                                                                                         <  8:1404890/149‑1 (MQ=255)
                                                                                                                                                                                         gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                                                                                                                                                                                        <  6:1409788/149‑1 (MQ=255)
                                                                                                                                                                                          aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                                                                                                                                                                                                >  2:351273/1‑124 (MQ=255)
                                                                                                                                                                                          aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAAGCACCGCCATTCAAAACGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACCAg                                                                          >  2:268293/1‑274 (MQ=255)
                                                                                                                                                                                            gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTACGCTACAGGATGGCCATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATATCCGCCATTCAAACAGCGTTAGAAGTGGGTTATCGCTCGCTTACTACCGCCGCGGCCTACAAGac                                                                        >  2:22432/1‑273 (MQ=255)
                                                                                                                                                                                            gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAAt                                                                                                                                                                                               >  1:52961/1‑155 (MQ=255)
                                                                                                                                                                                             ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATg                                                                                                                                                                                                    <  7:792897/149‑1 (MQ=255)
                                                                                                                                                                                              cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGTGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTTCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGAAAGTAATGAGGAATTAATCACCGACATTAAAAAAGCGTTAGAAGTGGGTTATCGCGCGATTGATACCGCCGCGTCCTACAAGAACga                                                                     >  1:43564/1‑275 (MQ=255)
                                                                                                                                                                                              cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATc                                                                                                                                                                                                                            <  4:351273/124‑1 (MQ=255)
                                                                                                                                                                                                cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCa                                                                                                                                                                                           <  3:52961/155‑1 (MQ=255)
                                                                                                                                                                                                   gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                                                                                                                                                                                              <  8:348790/149‑1 (MQ=255)
                                                                                                                                                                                                     tAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTc                                                                                                                                                                                            >  6:1123807/1‑149 (MQ=255)
                                                                                                                                                                                                     tAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTc                                                                                                                                                                                            >  5:647907/1‑149 (MQ=255)
                                                                                                                                                                                                         cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                                                                                                                                                                                        >  5:149020/1‑149 (MQ=255)
                                                                                                                                                                                                             ttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                                                                                                                                                                                    <  8:445092/149‑1 (MQ=255)
                                                                                                                                                                                                              tACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAg                                                                                                                                                                                   >  5:923526/1‑149 (MQ=255)
                                                                                                                                                                                                                 gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                                                                                                                                                                                              >  6:1353116/1‑103 (MQ=255)
                                                                                                                                                                                                                 gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg                                                                                                                                                                                >  6:1371911/1‑149 (MQ=255)
                                                                                                                                                                                                                   cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg                                                                                                                                                                              >  5:1427463/1‑149 (MQ=255)
                                                                                                                                                                                                                    tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                                                                                                                                                                                        <  8:720790/90‑1 (MQ=255)
                                                                                                                                                                                                                    tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                                                                                                                                                                                        >  7:720790/1‑90 (MQ=255)
                                                                                                                                                                                                                    tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                                                                                                                                             >  8:1329077/1‑149 (MQ=255)
                                                                                                                                                                                                                     cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc                                                                                                                                                                            <  8:280347/149‑1 (MQ=255)
                                                                                                                                                                                                                        aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                                                                                                                                         <  6:647907/149‑1 (MQ=255)
                                                                                                                                                                                                                          tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                                                                                                                       >  6:731000/1‑149 (MQ=255)
                                                                                                                                                                                                                              gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt                                                                                                                                                                   <  5:329365/149‑1 (MQ=255)
                                                                                                                                                                                                                               ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                                                                                                                                                                                   >  6:204293/1‑100 (MQ=255)
                                                                                                                                                                                                                               ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                                                                                                                                                                                   <  5:204293/100‑1 (MQ=255)
                                                                                                                                                                                                                               ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                                                                                                                                                  >  5:1460464/1‑149 (MQ=255)
                                                                                                                                                                                                                                      gCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGt                                                                                                                                                                     >  7:1263233/1‑139 (MQ=255)
                                                                                                                                                                                                                                       cATTTCGTCCAGACTTTAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  7:576135/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  5:113510/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  7:185477/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  6:93688/149‑1 (MQ=255)
                                                                                                                                                                                                                                           tCGTCCAGACTTAAGTTCACAACACCTCACTGGAGCTTGCTCCGGTGAGTTCATATAGAGGAAGAGCGTATGGCTAATCCATGCGTTATTAAGCTACAGGATGGCAATGTCATGCGCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTaa                                                                                                                                                      <  6:887643/149‑1 (MQ=255)
                                                                                                                                                                                                                                              tCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGa                                                                                                                                                   <  6:323458/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                               <  5:477091/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  gACTTAACTTCACAACACCTCACCGGAGCCTGCTGCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGATACAGGATGGCAATGTCATGCCCCAGCTGGGAGTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                               <  8:1236186/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  8:786901/135‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTa                                                                                                                                            >  7:1095075/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCa                                                                                                                                        >  7:1130981/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCa                                                                                                                                        >  7:307302/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ttCACAACCCCTCAACGGGGCCTGCTCCGGTGAGTTCAGATAAAGGAGGAAGGTATGGCTAATCCAACAGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGAGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATcac         <  2:105273/275‑1 (MQ=255)
                                                                                                                                                                                                                                                               aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCa                                                                                                                                  >  7:469099/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  agctgacgggggacTGGTCGGGTGAGGTAAGAGAAAGGGGGAACGTATGGCGAAGCCAACCGTTATGAAGAGAAAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGGTAGAAGTGGGTTATCGCTCGATTGAGACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACt     <  4:268293/258‑1 (MQ=255)
                                                                                                                                                                                                                                                                  acCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                                                                                                                       <  6:703973/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                  acCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                                                                                                                       >  5:703973/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                  aacTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTc                                                                                                                               >  8:204355/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                    cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATTCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaa                                                                                                                             >  8:484266/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt                                                                                                                                                                   <  8:92865/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaa                                                                                                                             <  6:273527/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaa                                                                                                                             >  8:830040/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAcc                                                                                                                                      <  7:1188939/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAcc                                                                                                                                      >  8:1188939/1‑139 (MQ=255)
                                                                                                                                                                                                                                                                     tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                            >  6:298156/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                            >  7:1041704/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGTATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc                                                                                                                        <  5:356585/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc                                                                                                                        <  7:487592/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGt                                                                                                                      >  5:1397935/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             gCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTa                                                                                                                    <  5:109068/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                                                                                                                                                                                               >  1:26328/1‑41 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                   
GTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAG  >  NC_000913/3156319‑3156850

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