Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,093,463 A→G 7.1% 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.1% 309.9 / 11.0 98S270S (AGT→AGCrhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (41/50);  new base G (5/2);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

AGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCA  >  NC_000913/4093318‑4093611
                                                                                                                                                 |                                                                                                                                                    
aGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTg                                                                                                                                                   >  2:653266/1‑149 (MQ=255)
aGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAAGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTg                                                                                                                                                   <  4:378981/149‑1 (MQ=255)
  ttATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGc                                                                                                                                                 <  1:510420/149‑1 (MQ=255)
    aTGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGa                                                                                                                                               >  3:19514/1‑149 (MQ=255)
     tGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAg                                                                                                                                              >  1:391883/1‑149 (MQ=255)
      gcgcACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGt                                                                                                                                             <  4:545281/149‑1 (MQ=255)
         cACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGc                                                                                                                                          <  1:98715/149‑1 (MQ=255)
          aCCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCg                                                                                                                                         >  3:332426/1‑149 (MQ=255)
            cGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTa                                                                                                                                       >  4:714301/1‑149 (MQ=255)
             gTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTaa                                                                                                                                      >  2:405667/1‑149 (MQ=255)
                   tttCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGa                                                                                                                                >  4:243684/1‑149 (MQ=255)
                   tttCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGa                                                                                                                                >  3:578500/1‑149 (MQ=255)
                        gCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGc                                                                                                                           >  2:170621/1‑149 (MQ=255)
                         cACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCt                                                                                                                          <  2:431827/149‑1 (MQ=255)
                          aCTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGGTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACAGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                         <  4:401258/149‑1 (MQ=255)
                          aCTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                         <  1:231768/149‑1 (MQ=255)
                           cTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGaaa                                                                                                                              >  4:473430/1‑143 (MQ=255)
                           cTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGaaa                                                                                                                              <  3:473430/143‑1 (MQ=255)
                           cTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTg                                                                                                                        >  2:651994/1‑149 (MQ=255)
                           cTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTg                                                                                                                        <  2:79/149‑1 (MQ=255)
                            tGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGc                                                                                                                       <  3:564604/149‑1 (MQ=255)
                             gCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGcc                                                                                                                      <  3:38034/149‑1 (MQ=255)
                              cTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCg                                                                                                                     >  4:534816/1‑149 (MQ=255)
                                      tCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTAt                                                                                                             <  1:651994/149‑1 (MQ=255)
                                        tGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCa                                                                                                           >  4:697231/1‑149 (MQ=255)
                                            cAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAActc                                                                                                       <  3:603991/149‑1 (MQ=255)
                                                tGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCa                                                                                                   <  4:182435/149‑1 (MQ=255)
                                                   gATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACgg                                                                                                >  4:128973/1‑149 (MQ=255)
                                                            gACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTc                                                                                       <  4:690152/149‑1 (MQ=255)
                                                             aCGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCt                                                                                      <  2:391883/149‑1 (MQ=255)
                                                               gCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                                    >  2:744848/1‑149 (MQ=255)
                                                                   aTACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCa                                                                                <  4:172077/149‑1 (MQ=255)
                                                                    tACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                               <  3:138267/149‑1 (MQ=255)
                                                                      ctctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgc                                                                             >  4:59643/1‑149 (MQ=255)
                                                                         ttCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTc                                                                                       <  1:752547/136‑1 (MQ=255)
                                                                         ttCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTc                                                                                       >  2:752547/1‑136 (MQ=255)
                                                                           cGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGccc                                                                        <  1:170621/149‑1 (MQ=255)
                                                                             tGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                      <  1:363504/149‑1 (MQ=255)
                                                                             tGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTCTCATGCCGCCCAt                                                                      <  3:239047/149‑1 (MQ=255)
                                                                              gggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                     <  2:404036/149‑1 (MQ=255)
                                                                              gggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                     <  2:61208/149‑1 (MQ=255)
                                                                              gggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                     >  4:121618/1‑149 (MQ=255)
                                                                              gggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                     <  1:43606/149‑1 (MQ=255)
                                                                              gggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                     <  2:298788/149‑1 (MQ=255)
                                                                                 cGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGaa                                                                  <  3:243684/149‑1 (MQ=255)
                                                                                 cGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGaa                                                                  <  1:310655/149‑1 (MQ=255)
                                                                                  gTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGGCGAGCGCTa                                                                                                              >  3:654446/1‑104 (MQ=255)
                                                                                  gTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTa                                                                                                              <  4:654446/104‑1 (MQ=255)
                                                                                     ggggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                              >  1:435201/1‑149 (MQ=255)
                                                                                      gggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                             <  4:91950/149‑1 (MQ=255)
                                                                                      gggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                             >  4:491569/1‑149 (MQ=255)
                                                                                        gTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCt                                                           >  2:587950/1‑149 (MQ=255)
                                                                                         tttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCAGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCtt                                                          <  1:405667/149‑1 (MQ=255)
                                                                                            aCCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGTGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCAc                                                       <  4:185421/149‑1 (MQ=255)
                                                                                             ccTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                               <  3:723845/124‑1 (MQ=255)
                                                                                             ccTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                               >  4:723845/1‑124 (MQ=255)
                                                                                               tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTaa                                                    <  1:386832/149‑1 (MQ=255)
                                                                                                 cATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTTTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATa                                                  >  3:657786/1‑149 (MQ=255)
                                                                                                    gACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACtt                                               >  1:620813/1‑149 (MQ=255)
                                                                                                       aaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgtg                                            >  4:101572/1‑149 (MQ=255)
                                                                                                        aaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGc                                           <  3:705825/149‑1 (MQ=255)
                                                                                                        aaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGc                                           <  1:653266/149‑1 (MQ=255)
                                                                                                        aaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGc                                           >  2:353828/1‑149 (MQ=255)
                                                                                                           gCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGa                                        >  4:439277/1‑149 (MQ=255)
                                                                                                                gCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAttttt                                   <  2:435201/149‑1 (MQ=255)
                                                                                                                   gATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTg                                <  4:695452/149‑1 (MQ=255)
                                                                                                                   gATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGCCTGTTTCATGCCGCCCATCGaa                                                                  <  1:230998/115‑1 (MQ=255)
                                                                                                                   gATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGCCTGTTTCATGCCGCCCATCGaa                                                                  >  2:230998/1‑115 (MQ=255)
                                                                                                                    aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgccgc                                                                          >  3:107349/1‑106 (MQ=255)
                                                                                                                    aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgccgc                                                                          <  4:107349/106‑1 (MQ=255)
                                                                                                                    aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTCATACTTGTGCTGAt                                       >  2:27347/1‑141 (MQ=255)
                                                                                                                    aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAt                                       <  1:27347/141‑1 (MQ=255)
                                                                                                                        tGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGa                                        <  3:78141/136‑1 (MQ=255)
                                                                                                                        tGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGa                                        >  4:78141/1‑136 (MQ=255)
                                                                                                                        tGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGgtg                           >  3:16895/1‑149 (MQ=255)
                                                                                                                            tGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGa                       <  2:110290/149‑1 (MQ=255)
                                                                                                                             gCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGAt                      <  3:121618/149‑1 (MQ=255)
                                                                                                                               tCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTa                    >  2:737224/1‑149 (MQ=255)
                                                                                                                                ccTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaa                   >  1:162751/1‑149 (MQ=255)
                                                                                                                                  ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc                 >  4:145511/1‑149 (MQ=255)
                                                                                                                                  ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc                 >  3:721628/1‑149 (MQ=255)
                                                                                                                                  ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc                 >  1:756762/1‑149 (MQ=255)
                                                                                                                                  gccactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc                 >  1:344684/7‑149 (MQ=255)
                                                                                                                                   tACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAcc                >  3:361745/1‑149 (MQ=255)
                                                                                                                                   ccactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAcc                >  2:388078/6‑149 (MQ=255)
                                                                                                                                     actcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                             >  2:737627/4‑102 (MQ=255)
                                                                                                                                     actcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                             <  1:737627/99‑1 (MQ=255)
                                                                                                                                     actcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGa                                        <  2:446444/120‑1 (MQ=255)
                                                                                                                                     actcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGa                                        >  1:446444/4‑123 (MQ=255)
                                                                                                                                      ctcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa             >  2:264982/3‑149 (MQ=255)
                                                                                                                                      agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCATGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa             <  1:507104/149‑1 (MQ=255)
                                                                                                                                      agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTTTGCTGATTTTTCTGCGGTGTCGATTAACCCGaa             >  2:380922/1‑149 (MQ=255)
                                                                                                                                      agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa             <  3:534816/149‑1 (MQ=255)
                                                                                                                                       gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg            <  1:30991/149‑1 (MQ=255)
                                                                                                                                       gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg            <  3:556095/149‑1 (MQ=255)
                                                                                                                                         gccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGt          >  1:467572/1‑149 (MQ=255)
                                                                                                                                           cagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGttt        <  4:578500/149‑1 (MQ=255)
                                                                                                                                                cACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATcc   <  4:19514/149‑1 (MQ=255)
                                                                                                                                                 aCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCa  <  2:467572/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCA  >  NC_000913/4093318‑4093611

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