Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,091,487 A→G 5.7% I131V (ATT→GTT)  hisD → bifunctional histidinal dehydrogenase/ histidinol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,091,4870AG 5.7% 220.6 / 3.0 70I131V (ATT→GTT) hisDbifunctional histidinal dehydrogenase/ histidinol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base G (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

AAACAGGCGATGGCGGTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAG  >  NC_000913/2091349‑2091633
                                                                                                                                          |                                                                                                                                                  
aaaCAGGCGATGGCGGTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTggc                                                                                                                                            <  6:74779/147‑1 (MQ=255)
        gATGGCGGTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGccc                                                                                                                                  >  5:31888/1‑149 (MQ=255)
         aTGGGGGTGGCAGTAATGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGGCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGcccc                                                                                                                                 <  4:351409/149‑1 (MQ=255)
          tGGCGGTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCg                                                                                                                                >  3:297067/1‑149 (MQ=255)
           tgCGGTGGCAGTAAAGAATATTGAATCCTTCCACACTGCGCAAAAACTGCCTAAGGTAGATGTAGAACCGCAGCAAGGAGTCCGTTTCCATCAGGTCACGCGTCAGGTAGATTCAGTTGGGTTGTATATTCCTGCCGGCTCCGCCACg                                                                                                                                >  6:191666/2‑148 (MQ=255)
           ggCGGTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGc                                                                                                                               >  5:269466/1‑149 (MQ=255)
               gTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGGCGTTAGATGTAGAATCGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTTGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCtctt                                                                                                                           <  5:191666/149‑1 (MQ=255)
                tGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCtcttc                                                                                                                          >  1:45233/1‑149 (MQ=255)
                  gCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCtcttct                                                                                                                         >  5:108537/1‑148 (MQ=255)
                        aaGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACTCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACgg                                                                                                                   >  7:406921/1‑148 (MQ=255)
                        aaGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACgg                                                                                                                   <  2:380750/148‑1 (MQ=255)
                            atatTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTa                                                                                                              <  6:31888/149‑1 (MQ=255)
                                tGAAACCTTCCCCACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGCTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACTGTATTAATGc                                                                                                          >  1:21181/1‑149 (MQ=255)
                                   aaCCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTgg                                                                                                       >  7:89394/1‑149 (MQ=255)
                                    aCCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGc                                                                                                      <  6:269466/149‑1 (MQ=255)
                                      cTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGa                                                                                                    <  3:217251/149‑1 (MQ=255)
                                          ctcacTGCGAAATAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCAAGCAGTTCACGGGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGTCCCGCTCTTCTCAACGTTATTAATGCTGGAGACTc                                                                                                 <  8:406921/146‑1 (MQ=255)
                                           acacTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCg                                                                                               >  2:79641/1‑149 (MQ=255)
                                            cacTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCgg                                                                                              >  1:18991/1‑149 (MQ=255)
                                             acTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTcc                                                                                                >  4:31353/1‑146 (MQ=255)
                                             acTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGc                                                                                             <  6:177313/149‑1 (MQ=255)
                                             acTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGGTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTcc                                                                                                <  3:31353/146‑1 (MQ=255)
                                                 cgcAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGa                                                                                           <  8:89394/147‑1 (MQ=255)
                                                    aaaaaCTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGctc                                                                                                                             >  1:90408/1‑110 (MQ=255)
                                                    aaaaaCTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGctc                                                                                                                             <  2:90408/110‑1 (MQ=255)
                                                    aaaaaCTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTAtt                                                                                      >  1:256267/1‑149 (MQ=255)
                                                       aaCTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCg                                                                                   >  4:230504/1‑149 (MQ=255)
                                                        aCTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGTTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCgg                                                                                  <  2:18991/149‑1 (MQ=255)
                                                              gccgGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATGTTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGCCTCCGGCGAGTATTGCGGGCTGt                                                                             >  5:232485/1‑148 (MQ=255)
                                                              gccgGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATGTTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGt                                                                             >  5:226790/1‑148 (MQ=255)
                                                               ccgGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTa                                                                            >  4:356376/1‑148 (MQ=255)
                                                                    aGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAg                                                                      >  4:378990/1‑149 (MQ=255)
                                                                      aTGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAgtg                                                                    <  3:378990/149‑1 (MQ=255)
                                                                      aTGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAgtg                                                                    <  1:79641/149‑1 (MQ=255)
                                                                         tAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTaa                                                                                                             <  1:402552/105‑1 (MQ=255)
                                                                         tAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTaa                                                                                                             >  2:402552/1‑105 (MQ=255)
                                                                           gAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgct                                                               <  3:356376/149‑1 (MQ=255)
                                                                            aaaCGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctg                                                              >  5:194563/1‑149 (MQ=255)
                                                                              aCGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctgtg                                                            >  8:61977/1‑149 (MQ=255)
                                                                                gCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATGTTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctgtgct                                                          <  6:226790/149‑1 (MQ=255)
                                                                                gCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgct                                                               <  2:230858/144‑1 (MQ=255)
                                                                                gCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgct                                                               >  1:230858/1‑144 (MQ=255)
                                                                                gCAGCCAGGCGTGCGGTGCCTGGAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATGGTCCTGGCGGCGCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctgtgct                                                          <  6:232485/149‑1 (MQ=255)
                                                                                      aGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCAccgc                                                    >  2:14388/1‑149 (MQ=255)
                                                                                             cgtTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGAt                                             <  8:63592/149‑1 (MQ=255)
                                                                                             cgtTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGAt                                             >  8:408960/1‑149 (MQ=255)
                                                                                               ttGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTg                                           >  6:172602/1‑149 (MQ=255)
                                                                                                tGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCGGTAAAAAAGTGGTGCTGTGCTCCCCGCCGCCGCTTGc                                          >  2:39014/1‑149 (MQ=255)
                                                                                                  ccAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCg                                        <  8:22110/149‑1 (MQ=255)
                                                                                                        ggTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATgaga                                  <  6:108537/149‑1 (MQ=255)
                                                                                                          tCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGGGGGCTGTAAAAAAGTGGTGGTGTGCTCACCGCCGCCGATTGCCGATGAGATc                                <  2:21181/149‑1 (MQ=255)
                                                                                                            aCGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCt                              >  5:253554/1‑149 (MQ=255)
                                                                                                            aCGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCt                              >  1:214661/1‑149 (MQ=255)
                                                                                                              gcgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGa                                              >  1:28053/1‑131 (MQ=255)
                                                                                                              gcgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGa                                              <  2:28053/131‑1 (MQ=255)
                                                                                                               cgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTa                           <  1:14388/149‑1 (MQ=255)
                                                                                                                   cGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATgcg                       >  6:363016/1‑149 (MQ=255)
                                                                                                                         cTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAg                 >  7:305059/1‑149 (MQ=255)
                                                                                                                          ttCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATgaga                                  >  7:88844/1‑131 (MQ=255)
                                                                                                                          ttCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATgaga                                  <  8:88844/131‑1 (MQ=255)
                                                                                                                               ttGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTTTTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGc           >  2:94313/1‑149 (MQ=255)
                                                                                                                               ttGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGc           >  4:203401/1‑149 (MQ=255)
                                                                                                                               ttGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGc           >  4:220299/1‑149 (MQ=255)
                                                                                                                               ttGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGc           >  8:351969/1‑149 (MQ=255)
                                                                                                                                 gggTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCg          <  7:408960/148‑1 (MQ=255)
                                                                                                                                 gggTTGTATATTCCTGGCGGCTACGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGTCTGTAAAAAAGTGGTGGTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCg          <  4:311960/148‑1 (MQ=255)
                                                                                                                                  ggTTGTATATTCATGGCGGCTCCGACCAGCTCTTCTCAACGGTATTAATGCTGGAGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGgt        <  2:214661/149‑1 (MQ=255)
                                                                                                                                    ttGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGgtgt      >  5:14392/1‑149 (MQ=255)
                                                                                                                                    ttGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGgtgt      <  2:45233/149‑1 (MQ=255)
                                                                                                                                      gTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGtgtg     >  8:79503/1‑148 (MQ=255)
                                                                                                                                       tataTTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGtgtg     >  2:191705/1‑147 (MQ=255)
                                                                                                                                       tataTTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGc    <  6:194563/148‑1 (MQ=255)
                                                                                                                                        atatTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAg  >  2:269632/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
AAACAGGCGATGGCGGTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAG  >  NC_000913/2091349‑2091633

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