Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,866,468 G→A 6.8% intergenic (‑195/+1) yidE ← / ← ibpB putative transporter/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,866,4680GA 6.8% 232.7 / 5.6 73intergenic (‑195/+1)yidE/ibpBputative transporter/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base A (2/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GAATCGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAG  >  NC_000913/3866321‑3866615
                                                                                                                                                   |                                                                                                                                                   
gAATCGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGt                                                                                                                                                    <  2:445552/149‑1 (MQ=255)
    cGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGc                                                                                                                                                <  3:379431/149‑1 (MQ=255)
     ggTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGGTAGTTAGCt                                                                                                                                               <  1:512592/149‑1 (MQ=255)
      gTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTa                                                                                                                                              <  1:310262/149‑1 (MQ=255)
          aGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttct                                                                                                                                          <  4:463031/149‑3 (MQ=255)
           gCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTaa                                                                                                                                         >  1:573255/1‑149 (MQ=255)
           gCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTaa                                                                                                                                         >  2:325919/1‑149 (MQ=255)
            ctgATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAAc                                                                                                                                        <  1:456233/147‑1 (MQ=255)
             aGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttctctt                                                                                                                                       <  2:164458/149‑6 (MQ=255)
              gATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGACGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATCGCTAGTTAGCTATTTAAcgc                                                                                                                                      >  4:563629/1‑149 (MQ=255)
                tttGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCgg                                                                                                                                    <  2:504343/149‑1 (MQ=255)
                  tGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGa                                                                                                                                  >  4:151682/1‑149 (MQ=255)
                     cTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGt                                                                                                                               >  2:174693/1‑149 (MQ=255)
                         taGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGc                                                                                                                           <  4:261433/149‑1 (MQ=255)
                           gcgcACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTg                                                                                                                         >  4:537869/1‑149 (MQ=255)
                           gcacacaTTATCACGGGAAAGAACGGGCCGAAATATCTTAAACAGTCGAAAGTTATCGCATTGTTGAGTAATTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTg                                                                                                                         >  4:491938/4‑149 (MQ=255)
                             gcACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAg                                                                                                                                                     <  1:260893/119‑1 (MQ=255)
                             gcACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAg                                                                                                                                                     >  2:260893/1‑119 (MQ=255)
                               acaTTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAg                                                                                                                     <  2:422607/149‑1 (MQ=255)
                                     tCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATAc                                                                                                               >  2:91018/1‑149 (MQ=255)
                                                aaCGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttctct                                                                                                                                        <  3:141153/113‑5 (MQ=255)
                                                aaCGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttctct                                                                                                                                        >  4:141153/1‑109 (MQ=255)
                                                  cGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGAt                                                                                                  <  1:70727/149‑1 (MQ=255)
                                                      ccGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGt                                                                                              >  2:11876/1‑149 (MQ=255)
                                                         aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           <  2:9218/149‑1 (MQ=255)
                                                         aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           <  2:573255/149‑1 (MQ=255)
                                                         aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           <  4:159886/149‑1 (MQ=255)
                                                         aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           >  2:337389/1‑149 (MQ=255)
                                                              tCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTc                                                                                      >  1:517904/1‑149 (MQ=255)
                                                                ttAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCAt                                                                                    >  3:377447/1‑149 (MQ=255)
                                                                ttAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCAt                                                                                    >  4:2658/1‑149 (MQ=255)
                                                                  aaaCAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTa                                                                                  >  2:105502/1‑149 (MQ=255)
                                                                    aCAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACg                                                                                >  4:486644/1‑149 (MQ=255)
                                                                      aGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGaa                                                                              <  1:109086/149‑1 (MQ=255)
                                                                          gAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTaa                                                                          <  1:91018/149‑1 (MQ=255)
                                                                            aaTTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAAt                                                                        <  4:215824/149‑1 (MQ=255)
                                                                               ttAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCaa                                                                     <  1:174693/149‑1 (MQ=255)
                                                                                  aCGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAatat                                                                  <  1:337389/149‑1 (MQ=255)
                                                                                       ttgttgAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGt                                                             >  1:275955/1‑149 (MQ=255)
                                                                                          ttgAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTaaa                                                          <  3:53302/149‑1 (MQ=255)
                                                                                            gAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAAcc                                                        >  4:442619/1‑149 (MQ=255)
                                                                                              gtaagtaaGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGt                                                      >  2:415574/1‑149 (MQ=255)
                                                                                                aagtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGttt                                                    >  4:63659/1‑149 (MQ=255)
                                                                                                 agtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCa                                                               <  2:354083/137‑1 (MQ=255)
                                                                                                 agtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCa                                                               >  1:354083/1‑137 (MQ=255)
                                                                                                 agtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTa                                                   <  1:15554/149‑1 (MQ=255)
                                                                                                  gtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTAc                                                  <  3:151682/149‑1 (MQ=255)
                                                                                                  gtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTAc                                                  >  1:362432/1‑149 (MQ=255)
                                                                                                  gtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTAc                                                  <  4:226526/149‑1 (MQ=255)
                                                                                                   taagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACg                                                 >  3:262809/1‑149 (MQ=255)
                                                                                                      gACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAg                                              <  3:442619/149‑1 (MQ=255)
                                                                                                      gACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAg                                              <  2:362432/149‑1 (MQ=255)
                                                                                                       aCAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAgg                                             >  2:186681/1‑149 (MQ=255)
                                                                                                        cAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGt                                            >  4:157813/1‑149 (MQ=255)
                                                                                                        cAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGt                                            >  2:414935/1‑149 (MQ=255)
                                                                                                         aaaaaaaGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGtt                                           <  1:186681/149‑1 (MQ=255)
                                                                                                           aaaaaGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTgc                                         >  1:28274/1‑149 (MQ=255)
                                                                                                                    caccaGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGAc                                <  3:305580/149‑1 (MQ=255)
                                                                                                                           aTGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCata                         <  3:537869/149‑1 (MQ=255)
                                                                                                                           aTGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCata                         <  2:275955/149‑1 (MQ=255)
                                                                                                                           aTGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCata                         >  1:503323/1‑149 (MQ=255)
                                                                                                                            tGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCatat                        <  1:536492/149‑1 (MQ=255)
                                                                                                                                ggggCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTc                    <  3:486644/149‑1 (MQ=255)
                                                                                                                                 gggCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTca                   >  4:148895/1‑149 (MQ=255)
                                                                                                                                 gggCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTca                   <  1:105502/149‑1 (MQ=255)
                                                                                                                                     aaaGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTcagcca               >  3:62428/1‑149 (MQ=255)
                                                                                                                                     aaaGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTcagcca               >  3:267414/1‑149 (MQ=255)
                                                                                                                                        ttaaaTAGCTAACTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCg            >  2:271971/4‑149 (MQ=255)
                                                                                                                                         agaATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGt           <  3:2658/149‑1 (MQ=255)
                                                                                                                                             tAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAg       <  1:11876/149‑1 (MQ=255)
                                                                                                                                              aGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGc      >  2:123728/1‑149 (MQ=255)
                                                                                                                                               gCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCt     >  4:110182/1‑149 (MQ=255)
                                                                                                                                                 tAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCa   <  3:252450/149‑1 (MQ=255)
                                                                                                                                                 tAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCa   >  1:76368/1‑149 (MQ=255)
                                                                                                                                                  aGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAg  <  4:76368/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAATCGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAG  >  NC_000913/3866321‑3866615

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