Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,585,149 Δ1 bp 5.3% coding (1338/2280 nt) ydeP ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,585,1490G. 5.3% 390.3 / 13.9 113coding (1338/2280 nt)ydePputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (53/54);  new base . (2/4);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

CGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCG  >  NC_000913/1585003‑1585297
                                                                                                                                                  |                                                                                                                                                     
cGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGc                                                                                                                                                     >  3:615763/1‑149 (MQ=255)
  gTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCaa                                                                                                                                                   >  3:340398/1‑149 (MQ=255)
  gTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCaa                                                                                                                                                   >  1:5357/1‑149 (MQ=255)
       gCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgca                                                                                                                                             >  4:286911/1‑149 (MQ=255)
       gCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgca                                                                                                                                             <  2:849670/149‑1 (MQ=255)
           aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACTGTACAGCGCTCGCTTCCCGATCTGGCATTGCCACCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTGGCATGCTGCAATTGCCAGCAt                                                                                                                                         <  2:162291/149‑1 (MQ=255)
           aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCAt                                                                                                                                         <  1:138599/149‑1 (MQ=255)
             tGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGCAATTGCCAGCATgc                                                                                                                                       <  2:844415/149‑2 (MQ=255)
               gCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGt                                                                                                                                       <  3:164364/147‑1 (MQ=255)
               gCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTcc                                                                                                                                     <  3:58238/149‑1 (MQ=255)
                cTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCa                                                                                                                                    >  4:214807/1‑149 (MQ=255)
                 tACTTGTACCTCCAAATGTAATTGCGTTAACGGTACAGCGCTTGTTTCCCGATCTGGCATTGCCAGATCAAAATTGCCCCCCATACAGATCAATGCTGTAGCCTGCCCTGTACATTTTGCTGGAATGCTGGCAATTGCAGCATGTCCAg                                                                                                                                   <  4:447498/149‑1 (MQ=255)
                 tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAg                                                                                                                                   >  4:511731/1‑149 (MQ=255)
                 tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAg                                                                                                                                   <  1:353943/149‑1 (MQ=255)
                 tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAg                                                                                                                                   >  3:540793/1‑149 (MQ=255)
                       tACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCAt                                                                                                                             <  3:452318/149‑1 (MQ=255)
                         ccGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGa                                                                                                                           <  2:540091/149‑1 (MQ=255)
                           gCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAgg                                                                                                                         >  1:890464/1‑149 (MQ=255)
                           gCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAgg                                                                                                                         >  3:82627/1‑149 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGt                                                                                                                                 >  1:471071/1‑140 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGt                                                                                                                                 <  2:471071/140‑1 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAgg                                                                                                                         >  3:613756/1‑148 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAgg                                                                                                                         <  4:613756/148‑1 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAgg                                                                                                                         <  2:370370/148‑1 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAgg                                                                                                                         >  1:370370/1‑148 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                        >  3:180124/1‑149 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                        >  4:754589/1‑149 (MQ=255)
                            ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                        >  1:841770/1‑149 (MQ=255)
                              aaaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTgg                                                                                                                      <  3:286911/149‑1 (MQ=255)
                                aTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTgggg                                                                                                                    <  4:340398/149‑1 (MQ=255)
                                 tCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGt                                                                                                                   <  1:699561/149‑1 (MQ=255)
                                 tCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGt                                                                                                                   >  2:679908/1‑149 (MQ=255)
                                    aaTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGaa                                                                                                                <  2:387327/149‑1 (MQ=255)
                                     aTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                               >  2:670744/1‑149 (MQ=255)
                                     aTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                               >  4:806305/1‑149 (MQ=255)
                                      ttGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                                                              >  4:213112/1‑149 (MQ=255)
                                      ttGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                                                              >  4:552413/1‑149 (MQ=255)
                                          gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCGAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                          <  3:754589/149‑1 (MQ=255)
                                          gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                          <  4:782505/149‑1 (MQ=255)
                                           ttAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAg                                                                                                         >  4:486470/1‑149 (MQ=255)
                                              aCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcgc                                                                                                      <  4:818402/149‑1 (MQ=255)
                                                  tACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                                                                  >  4:716330/1‑149 (MQ=255)
                                                     aGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCa                                                                                               >  4:678184/1‑149 (MQ=255)
                                                              ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTAGATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                      >  4:62163/1‑149 (MQ=255)
                                                              ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                                                              >  3:561911/1‑125 (MQ=255)
                                                              ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                      >  2:892210/1‑149 (MQ=255)
                                                              ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                      <  2:513901/149‑1 (MQ=255)
                                                              ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                      >  2:849610/1‑149 (MQ=255)
                                                                cccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATCGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGa                                                                                    >  4:607939/1‑149 (MQ=255)
                                                                  tgaccTGGCATTGCTAGCGCAAAATTGCCCCTCATGCAGATCAATGCTCGAGCCTGCCCCGTATATATGGCTTGCATGCTGGCAATAGCAGCGTGTCCAGGTGCACGAGGTGCGGTGAAGCCATAGAGCTCACCCAGACGAGGCAGaaa                                                                                  <  3:545797/145‑1 (MQ=255)
                                                                  cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                  <  3:618404/149‑1 (MQ=255)
                                                                  cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                  <  4:146690/149‑1 (MQ=255)
                                                                  cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                  <  1:546625/149‑1 (MQ=255)
                                                                    aTCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAct                                                                                >  2:456277/1‑149 (MQ=255)
                                                                      cTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActct                                                                              <  1:488168/149‑1 (MQ=255)
                                                                      cTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActct                                                                              <  2:583310/149‑1 (MQ=255)
                                                                        ggCATTGCCAGCGCAAAATTGCCCCCCTTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGc                                                                            <  1:679908/149‑1 (MQ=255)
                                                                         gCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa                                                                           <  1:5411/149‑1 (MQ=255)
                                                                           aTTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGa                                                                         >  3:536345/1‑149 (MQ=255)
                                                                                   cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtct                                                                 <  4:513943/149‑1 (MQ=255)
                                                                                     cAAAACTGCCACCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGATTCTCg                                                               <  4:498566/149‑1 (MQ=255)
                                                                                      aaaaTCGACCCCAAAGCAGAACAATTCGCGAGTGTGCCCCGGACATATTGGTTGCATCATGGCAATTGCAGCATGTCCAGGTCCGTGAGGTGGTGCCAAGCCTGAGGGTTAATCCAGATGAGCCGGAAACTCAGCAGACGGTTTCTCgg                                                              <  2:445907/149‑1 (MQ=255)
                                                                                       aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                             <  1:558426/149‑1 (MQ=255)
                                                                                       aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                             <  1:131653/149‑1 (MQ=255)
                                                                                             ccccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATAcc                                                       <  2:841770/149‑1 (MQ=255)
                                                                                                cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGAc                                                    >  1:862081/1‑149 (MQ=255)
                                                                                                cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGAc                                                    <  4:759518/149‑1 (MQ=255)
                                                                                                 ccATGCAGATCTATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACg                                                   >  3:514412/1‑149 (MQ=255)
                                                                                                 ccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACg                                                   >  1:116290/1‑149 (MQ=255)
                                                                                                  cATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACgg                                                  <  1:629986/149‑1 (MQ=255)
                                                                                                        gATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTCCGGt                                            >  4:869895/1‑149 (MQ=255)
                                                                                                            aaTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcc                                                                                                 <  4:64096/92‑1 (MQ=255)
                                                                                                            aaTGCTCGAGCCTGACCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcc                                                                                                 >  3:64096/1‑92 (MQ=255)
                                                                                                             aTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGccc                                       >  4:269884/1‑149 (MQ=255)
                                                                                                                cTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGt                                    >  2:376/1‑149 (MQ=255)
                                                                                                                   gAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACa                                 <  2:115861/149‑1 (MQ=255)
                                                                                                                    aGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTAACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACAt                                <  2:814846/149‑1 (MQ=255)
                                                                                                                     gCCTGCCCCGTACATATTGGTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACAtt                               >  3:696459/1‑149 (MQ=255)
                                                                                                                        tGCCCCGTACATATTGCTTGGATGCGGGCAATTGCAGCATGTCCAGGTGCATGAGGGGGGGTGAAGCCATAGCGCTCACCCATACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            <  4:173247/149‑1 (MQ=255)
                                                                                                                        tGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATTTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            <  3:759870/149‑1 (MQ=255)
                                                                                                                          ccccGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgt                          <  3:213112/149‑1 (MQ=255)
                                                                                                                           cccGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgtg                         >  4:579932/1‑149 (MQ=255)
                                                                                                                            ccGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTg                                     <  2:818958/136‑1 (MQ=255)
                                                                                                                              gTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTcc                      <  2:516956/149‑1 (MQ=255)
                                                                                                                                 cATATTGCTTGCATGCTGGCAATTGCAGTATGTCCAGGTGTATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGCAACACTGCCGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCCCg                   <  4:238140/149‑1 (MQ=255)
                                                                                                                                  atatTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGt                  >  4:156958/1‑149 (MQ=255)
                                                                                                                                     ttgcttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGt               >  2:469229/1‑149 (MQ=255)
                                                                                                                                     ttgcttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGt               <  2:862081/149‑1 (MQ=255)
                                                                                                                                         ttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa                                                           <  1:485912/101‑1 (MQ=255)
                                                                                                                                         ttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGc           >  2:143487/1‑149 (MQ=255)
                                                                                                                                         ggacATGCTGCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            >  3:443452/4‑132 (MQ=255)
                                                                                                                                         ggacATGCTGCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            <  4:443452/129‑1 (MQ=255)
                                                                                                                                          tgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            >  4:636075/1‑131 (MQ=255)
                                                                                                                                          tgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            <  3:636075/131‑1 (MQ=255)
                                                                                                                                            cATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGa        <  3:678184/149‑1 (MQ=255)
                                                                                                                                             aTGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGAt       >  2:188243/1‑149 (MQ=255)
                                                                                                                                             aTGCTGCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                             <  3:107012/95‑1 (MQ=255)
                                                                                                                                             aTGCTGCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                             >  4:107012/1‑95 (MQ=255)
                                                                                                                                              tGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATa      >  2:660243/1‑149 (MQ=255)
                                                                                                                                                  ggCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCg  >  3:785616/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCG  >  NC_000913/1585003‑1585297

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