Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,093,463 A→G 7.3% S270S (AGT→AGC rhaD ← rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,093,4630AG 7.3% 397.6 / 15.8 123S270S (AGT→AGCrhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (61/53);  new base G (6/3);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 127 total aligned reads displayed.

GTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCA  >  NC_000913/4093316‑4093611
                                                                                                                                                   |                                                                                                                                                    
gTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCAc                                                                                                                                                     >  4:847906/1‑149 (MQ=255)
gTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCAc                                                                                                                                                     >  1:557538/1‑149 (MQ=255)
gTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGACTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTTCCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCAc                                                                                                                                                     >  3:798784/1‑149 (MQ=255)
  aGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTg                                                                                                                                                   <  2:716738/149‑1 (MQ=255)
  aGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTg                                                                                                                                                   <  4:333701/149‑1 (MQ=255)
     tATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCg                                                                                                                                                >  3:94703/1‑149 (MQ=255)
        gcgcACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGt                                                                                                                                             <  1:175040/149‑1 (MQ=255)
            aCCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCg                                                                                                                                         >  3:275417/1‑149 (MQ=255)
                     tttCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGa                                                                                                                                >  1:41200/1‑149 (MQ=255)
                      ttCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGaa                                                                                                                               >  1:45094/1‑149 (MQ=255)
                       tCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGaaa                                                                                                                              >  3:705126/1‑149 (MQ=255)
                        cAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAAc                                                                                                                             >  4:133135/1‑149 (MQ=255)
                          gCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGc                                                                                                                           >  1:807485/1‑149 (MQ=255)
                            aCTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                         <  4:193897/149‑1 (MQ=255)
                                 ttCCAGTTCTGGCCAGATGGGATTATGCGGAGGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGa                                                                                                                    >  4:566905/1‑149 (MQ=255)
                                 ttCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGa                                                                                                                    >  2:828676/1‑149 (MQ=255)
                                  tCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAg                                                                                                                   >  4:872624/1‑149 (MQ=255)
                                      gTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCt                                                                                                               <  3:87788/149‑1 (MQ=255)
                                            gCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAAc                                                                                                         >  2:107352/1‑149 (MQ=255)
                                            gCCAGATGAGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTTCCTGCATGACAGAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAAc                                                                                                         <  4:568896/149‑1 (MQ=255)
                                              cAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAActc                                                                                                       <  3:551984/149‑1 (MQ=255)
                                                gATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCtt                                                                                                     >  3:203235/1‑149 (MQ=255)
                                                           gCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATg                                                                                          <  2:632066/149‑1 (MQ=255)
                                                            ccGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATgg                                                                                         >  3:220742/1‑149 (MQ=255)
                                                             cGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCa                                                                                                                                                      >  4:278497/1‑87 (MQ=255)
                                                             cGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGt                                                                                        >  4:389475/1‑149 (MQ=255)
                                                                 gCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                                    >  2:275217/1‑149 (MQ=255)
                                                                 gCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                                    >  2:57855/1‑149 (MQ=255)
                                                                   tGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGttt                                                                                  >  3:190490/1‑149 (MQ=255)
                                                                   tGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGttt                                                                                  <  4:319517/149‑1 (MQ=255)
                                                                   tGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGttt                                                                                  >  1:791101/1‑149 (MQ=255)
                                                                     aTACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCa                                                                                >  1:282571/1‑149 (MQ=255)
                                                                      tACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                               <  2:421594/149‑1 (MQ=255)
                                                                         tctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGTGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCt                                                                                                                          >  4:532692/1‑103 (MQ=255)
                                                                         tctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCt                                                                                                                          <  3:532692/103‑1 (MQ=255)
                                                                         tctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgcc                                                                            >  3:917984/1‑149 (MQ=255)
                                                                          ctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgccg                                                                           >  3:218689/1‑149 (MQ=255)
                                                                           ttCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGttt                                                                                  >  4:174664/1‑141 (MQ=255)
                                                                           ttCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGttt                                                                                  <  3:174664/141‑1 (MQ=255)
                                                                                gggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                     <  4:203235/149‑1 (MQ=255)
                                                                                gggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                     >  3:140072/1‑149 (MQ=255)
                                                                                  gCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGa                                                                   <  1:790405/149‑1 (MQ=255)
                                                                                   cGTCGGGGTTTACCTGCATTACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGaa                                                                  <  4:629989/149‑1 (MQ=255)
                                                                                   cGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGaa                                                                  <  3:847906/149‑1 (MQ=255)
                                                                                    gTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAAt                                                                 <  4:545/149‑1 (MQ=255)
                                                                                       ggggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                              >  2:480054/1‑149 (MQ=255)
                                                                                       ggggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                              <  1:89853/149‑1 (MQ=255)
                                                                                           tttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCtt                                                          <  4:179351/149‑1 (MQ=255)
                                                                                            ttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTc                                                         <  2:24940/149‑1 (MQ=255)
                                                                                                cTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTa                                                     >  1:298958/1‑149 (MQ=255)
                                                                                                cTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTa                                                     >  3:243099/1‑149 (MQ=255)
                                                                                                 tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATTAACCTTCACTaa                                                    >  1:739567/1‑149 (MQ=255)
                                                                                                 tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTaa                                                    >  3:852158/1‑149 (MQ=255)
                                                                                                 ggCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTaa                                                    >  2:793377/2‑149 (MQ=255)
                                                                                                  gCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAAt                                                   <  4:140072/149‑1 (MQ=255)
                                                                                                   cATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGaa                                                                  <  3:342148/133‑1 (MQ=255)
                                                                                                   cATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGaa                                                                  >  4:342148/1‑133 (MQ=255)
                                                                                                      gACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACtt                                               >  1:517281/1‑149 (MQ=255)
                                                                                                        cgaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgt                                             >  3:234318/3‑149 (MQ=255)
                                                                                                         aaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgtg                                            <  4:852203/149‑1 (MQ=255)
                                                                                                         aaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgtg                                            >  1:514115/1‑149 (MQ=255)
                                                                                                         aaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgtg                                            <  1:480054/149‑1 (MQ=255)
                                                                                                         aaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgtg                                            >  2:868196/1‑149 (MQ=255)
                                                                                                         aaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgtg                                            >  4:658056/1‑149 (MQ=255)
                                                                                                             gCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGa                                        <  3:133135/149‑1 (MQ=255)
                                                                                                               tttGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                             <  3:322038/126‑1 (MQ=255)
                                                                                                               tttGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                             >  4:322038/1‑126 (MQ=255)
                                                                                                                    ggATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCt                                 <  4:190490/149‑1 (MQ=255)
                                                                                                                    ggATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCt                                 <  3:872624/149‑1 (MQ=255)
                                                                                                                    ggATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCt                                 >  1:41166/1‑149 (MQ=255)
                                                                                                                      aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGc                               >  3:40182/1‑149 (MQ=255)
                                                                                                                       tCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCg                              <  3:180004/149‑1 (MQ=255)
                                                                                                                         aTGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGgt                            <  2:108744/149‑1 (MQ=255)
                                                                                                                           gTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGgtgt                          <  4:664889/149‑1 (MQ=255)
                                                                                                                            tttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGgt                            <  4:197478/146‑1 (MQ=255)
                                                                                                                            tttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTc                         <  3:725867/149‑1 (MQ=255)
                                                                                                                            tttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTc                         >  2:923735/1‑149 (MQ=255)
                                                                                                                              tGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGa                       <  1:868196/149‑1 (MQ=255)
                                                                                                                                cTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGAtt                     >  1:726196/1‑149 (MQ=255)
                                                                                                                                   cTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATa                                                  <  2:638416/117‑1 (MQ=255)
                                                                                                                                   cTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATa                                                  >  1:638416/1‑117 (MQ=255)
                                                                                                                                    ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc                 >  2:130507/1‑149 (MQ=255)
                                                                                                                                      cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa                  <  1:273097/142‑1 (MQ=255)
                                                                                                                                      cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa                  <  4:180894/142‑1 (MQ=255)
                                                                                                                                      cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa                  >  3:180894/5‑146 (MQ=255)
                                                                                                                                      cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa                  >  2:273097/5‑146 (MQ=255)
                                                                                                                                       actcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgt                                             >  3:298575/4‑118 (MQ=255)
                                                                                                                                       actcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgt                                             <  4:298575/115‑1 (MQ=255)
                                                                                                                                        agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa             <  1:408431/149‑1 (MQ=255)
                                                                                                                                         tcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg            >  3:668509/2‑149 (MQ=255)
                                                                                                                                         tcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg            >  3:542516/2‑149 (MQ=255)
                                                                                                                                         gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACACTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg            >  3:245345/1‑149 (MQ=255)
                                                                                                                                           gccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGAt                      >  1:784584/1‑137 (MQ=255)
                                                                                                                                           gccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGAt                      <  2:784584/137‑1 (MQ=255)
                                                                                                                                           gccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGt          <  2:130414/149‑1 (MQ=255)
                                                                                                                                           gccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGt          >  2:687965/1‑149 (MQ=255)
                                                                                                                                            ccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGtt         <  4:705126/149‑1 (MQ=255)
                                                                                                                                              agcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTc       <  2:557538/149‑1 (MQ=255)
                                                                                                                                                cgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCAt     <  2:298958/149‑1 (MQ=255)
                                                                                                                                                   aCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCa  <  2:164618/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCA  >  NC_000913/4093316‑4093611

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