Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 5.4% intergenic (+72/‑33) 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,5900CA 5.4% 791.2 / 24.3 221intergenic (+72/‑33)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (98/110);  new base A (5/7);  total (103/118)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-01

Only 100 of 225 total aligned reads displayed.

ATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAG  >  NC_000913/3156321‑3156850
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                    
          caaTTGCCGGAAACCGCAATTGCTTTGAGCAATTAGGCGTGCCTACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                         <  2:300462/271‑1 (MQ=255)
                            aTTTCTTTGGGCAATTAGGCGGGCCAGCCCACCTCTCCGACTACGGGTTGGACGGCAGCGCCATCCCGGCTTTGCTGAAAAAACTTGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGt                                                                                                                                                                                                                                      <  4:152120/274‑1 (MQ=255)
                                ttttGCGCAATTAGGTGTGCCGACCCACCTCTCCGACTACGGTGTGGACGGCAGCTCCATCCAGGCTTTGCTGAAAAAACTGGAAGGGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                                                                                                                                                                                 <  4:237380/274‑1 (MQ=255)
                                               cGTGCCGACCCACCTCTCCGAGTACGGTCTGGACGGGAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                                                                                                                                                                                  <  2:60736/275‑1 (MQ=255)
                                                     gACCCTACTCTCCGACGACGGTCTGGACGGCAGCTCCATCCCGGTGTTGCTGAAAAAACTGGAAGGGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCATCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                                                                                                                                                                                             <  4:80550/274‑1 (MQ=255)
                                                                 ctactcagggcTGGACGGCAGCTCCAGCCCGGCTTTGCTGAAAAAACTGGAAGGACACGTCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                                                                                                                                 <  1:28468/264‑1 (MQ=255)
                                                                         gTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAggagga                                                                                                                                                                                                                                          >  2:215558/1‑225 (MQ=255)
                                                                                           tCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGcc                                                                                                                                                                                                                                                                      >  1:21447/1‑179 (MQ=255)
                                                                                               ggCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCt                                                                                                                                                                                                                                                                  <  3:21447/179‑1 (MQ=255)
                                                                                                   ttGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATTTCATGCCCCAGCTGGGACTGGGCgccgggaaa                                                                                                                                                              >  4:328296/1‑267 (MQ=255)
                                                                                                                   aaGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAATCCGCCCGCTAAGCTTTTTACGCCTCAAACGTTCGTTTTCGGGCATTTCGTCCAGACGTAAGTTCACAACACCTCAACGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAg                                                                                                                                              <  2:78827/275‑1 (MQ=255)
                                                                                                                              aTGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTc                                                                                                                                                                                                                                                                 >  5:273527/1‑149 (MQ=255)
                                                                                                                                gACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCg                                                                                                                                                                                                                                                               <  7:267365/149‑1 (MQ=255)
                                                                                                                                 aCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCgg                                                                                                                                                                                                                                                              <  7:776323/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)
                                                                                                                                     aaCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                                                                                                                                                                                          <  8:632904/149‑1 (MQ=255)
                                                                                                                                      aCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                         >  7:1340632/1‑149 (MQ=255)
                                                                                                                                          ggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGTGAggtgtt                                                                                                                                                                                                                                                     <  8:1415693/149‑6 (MQ=255)
                                                                                                                                           gCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                                                                                                                                                                                    <  5:765633/149‑1 (MQ=255)
                                                                                                                                           gCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                                                                                                                                                                                    >  6:1360370/1‑149 (MQ=255)
                                                                                                                                             gAAAATCATGAGATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                                                                                                                                                                                  >  7:1145169/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)
                                                                                                                                                         aTTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCCCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                             >  8:1057916/1‑126 (MQ=255)
                                                                                                                                                           tACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGGTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATc                                                                                                       <  2:1822/274‑1 (MQ=255)
                                                                                                                                                            aCGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                                                                                                                                                                                               >  3:46422/1‑153 (MQ=255)
                                                                                                                                                             cGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                                                                                                                                                                                  <  5:1067905/149‑1 (MQ=255)
                                                                                                                                                              gTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTTAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                                                                                                                                                                                 >  8:576135/1‑149 (MQ=255)
                                                                                                                                                                tGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                                                                                                                                                                                               >  7:348790/1‑149 (MQ=255)
                                                                                                                                                                tGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                                                                                                                                                                                           <  1:46422/153‑1 (MQ=255)
                                                                                                                                                                  gATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGtt                                                                                                                                                                                                                                                        >  4:3007/1‑122 (MQ=255)
                                                                                                                                                                      tCAGCCGCCGTATATACGAAGCCGCCAGCTAAGCTTTTTACGCCTCAAACTTTCGTATTCGGGCATTTCGTCCAAACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                                                                                                                                                                                         >  7:904660/1‑149 (MQ=255)
                                                                                                                                                                        aGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                                                                                                                                                                                       <  5:1187446/149‑1 (MQ=255)
                                                                                                                                                                             ccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGGTATCGCTCGATTGATACcgccgc                                                                                    >  1:38843/1‑275 (MQ=255)
                                                                                                                                                                                tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                                                                                                                                                                                               >  5:1151413/1‑149 (MQ=255)
                                                                                                                                                                                      cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACa                                                                                                                                                                                                         <  8:1404890/149‑1 (MQ=255)
                                                                                                                                                                                       gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                                                                                                                                                                                        <  6:1409788/149‑1 (MQ=255)
                                                                                                                                                                                       gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                                                                                                                                                                                        >  7:501514/1‑149 (MQ=255)
                                                                                                                                                                                        aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAAGCACCGCCATTCAAAACGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACCAg                                                                          >  2:268293/1‑274 (MQ=255)
                                                                                                                                                                                           ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATg                                                                                                                                                                                                    <  7:792897/149‑1 (MQ=255)
                                                                                                                                                                                                 gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                                                                                                                                                                                              <  8:348790/149‑1 (MQ=255)
                                                                                                                                                                                                   tAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTc                                                                                                                                                                                            >  6:1123807/1‑149 (MQ=255)
                                                                                                                                                                                                       cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                                                                                                                                                                                        >  5:1482323/1‑149 (MQ=255)
                                                                                                                                                                                                           ttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                                                                                                                                                                                    <  8:445092/149‑1 (MQ=255)
                                                                                                                                                                                                            tACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAg                                                                                                                                                                                   >  5:923526/1‑149 (MQ=255)
                                                                                                                                                                                                               gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                                                                                                                                                                                              <  5:1353116/103‑1 (MQ=255)
                                                                                                                                                                                                                 cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGaa                                                                                                                                                                                                                                         <  5:1111367/90‑1 (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                                                                                                                                                                            >  5:321524/1‑149 (MQ=255)
                                                                                                                                                                                                                   cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc                                                                                                                                                                            <  8:280347/149‑1 (MQ=255)
                                                                                                                                                                                                                    aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                                                                                                                                                                                              >  7:1235577/1‑98 (MQ=255)
                                                                                                                                                                                                                      aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                                                                                                                                         <  6:647907/149‑1 (MQ=255)
                                                                                                                                                                                                                        tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                                                                                                                       >  6:731000/1‑149 (MQ=255)
                                                                                                                                                                                                                            gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt                                                                                                                                                                   <  6:222941/149‑1 (MQ=255)
                                                                                                                                                                                                                             ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                                                                                                                                                                                   <  5:204293/100‑1 (MQ=255)
                                                                                                                                                                                                                             ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                                                                                                                                                  >  5:1460464/1‑149 (MQ=255)
                                                                                                                                                                                                                                  gggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaag                                                                                                                                                             >  6:1009933/1‑149 (MQ=255)
                                                                                                                                                                                                                                   ggCATTTGGTCCAGACTAAAGTTCACAACACCTCACCGGAGCCTGTTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCaa                              <  2:328296/275‑1 (MQ=255)
                                                                                                                                                                                                                                     cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  7:185477/149‑1 (MQ=255)
                                                                                                                                                                                                                                     cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  5:687045/149‑1 (MQ=255)
                                                                                                                                                                                                                                     cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  6:93688/149‑1 (MQ=255)
                                                                                                                                                                                                                                            tCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGa                                                                                                                                                   <  5:137204/149‑1 (MQ=255)
                                                                                                                                                                                                                                            tCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGa                                                                                                                                                   <  6:323458/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                               >  6:775156/1‑149 (MQ=255)
                                                                                                                                                                                                                                                gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                               <  5:477091/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gACTTAACTTCACAACACCTCACCGGAGCCTGCTGCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGATACAGGATGGCAATGTCATGCCCCAGCTGGGAGTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                               <  8:1236186/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                                                                                                                                          <  8:786901/135‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                               <  6:196158/146‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTa                                                                                                                                            >  7:1095075/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        ttCACAACCCCTCAACGGGGCCTGCTCCGGTGAGTTCAGATAAAGGAGGAAGGTATGGCTAATCCAACAGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGAGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATcac         <  2:105273/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          cacaACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCg                                                                                                                                     >  6:1212937/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCa                                                                                                                                  >  7:469099/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCa                                                                                                                                  >  5:348417/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cacCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCg                                                                                                      <  2:131254/175‑1 (MQ=255)
                                                                                                                                                                                                                                                               aacctacccGGAGCCTGTTCCAGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCAAACCGGTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTTGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAg  <  2:318380/268‑1 (MQ=255)
                                                                                                                                                                                                                                                                acCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                                                                                                                       <  6:703973/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                aacTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTc                                                                                                                               >  8:204355/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                  cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATTCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaa                                                                                                                             >  8:484266/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt                                                                                                                                                                   <  8:92865/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaa                                                                                                                             >  8:830040/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaa                                                                                                                             <  6:273527/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAcc                                                                                                                                      >  8:1188939/1‑139 (MQ=255)
                                                                                                                                                                                                                                                                   tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAcc                                                                                                                                      <  7:1188939/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                            <  7:1120776/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                            >  6:298156/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                            >  7:1041704/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaa                                                                                                                           >  5:161382/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaaa                                                                                                                          <  5:1137946/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc                                                                                                                        <  7:487592/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGAGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGt                                                                                                                      <  5:977139/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGtt                                                                                                                     <  5:621989/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAg                                                                                                                   >  6:394362/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                    
ATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAG  >  NC_000913/3156321‑3156850

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