Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,723,266 G→A 5.4% W382* (TGG→TGA)  ydhK → putative efflux protein (PET) component of YdhJK efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,723,2660GA 5.4% 238.7 / 2.4 74W382* (TGG→TGA) ydhKputative efflux protein (PET) component of YdhJK efflux pump
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base A (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GTCTGCAAAGCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTGGTCTGA  >  NC_000913/1723118‑1723411
                                                                                                                                                    |                                                                                                                                                 
gTCTGCAAAGCGGTGTAGATGATCATACCAGACGCCCACGGACATCCGGCCTGGCTCGTCATACCGCTAACGCCGAAGCGCGGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATgg                                                                                                                                                   <  3:187052/149‑1 (MQ=255)
   tGCAAAGCGGTGTAGATGATCATACCAGACTCACACGGACATCCGGCCTGTCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGAt                                                                                                                                                >  2:218761/1‑149 (MQ=255)
     cAAAGCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGa                                                                                                                                                 <  5:394226/146‑1 (MQ=255)
     cAAAGCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATg                                                                                                                                               <  3:172589/148‑1 (MQ=255)
     cAAAGCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGc                                                                                                                                              >  5:276082/1‑149 (MQ=255)
       aaGCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCg                                                                                                                                            <  2:4112/149‑1 (MQ=255)
        aGCGGTGTAGATGATCATACCAGACTCCCCCGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGAGATGCCgg                                                                                                                                           >  3:156273/1‑149 (MQ=255)
        aGCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGCGCGGGCTTCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGAGATGcggt                                                                                                                                           >  3:162244/1‑146 (MQ=255)
         gCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGt                                                                                                                                          <  5:321268/149‑1 (MQ=255)
           ggTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGc                                                                                                                                        <  2:109776/149‑1 (MQ=255)
               tAGATTATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATCCCGATAACGCCGCAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATTGAGTATTGCTTCGCAATGGGATGCCGGTGCCAAt                                                                                                                                    <  1:83230/149‑1 (MQ=255)
                    gatCATACAAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCAt                                                                                                                                >  6:105402/1‑148 (MQ=255)
                        aTACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTAGATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAAc                                                                                                                            <  3:227063/148‑1 (MQ=255)
                            cAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                                                                                                                       <  3:283323/149‑1 (MQ=255)
                            cAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                                                                                                                       >  3:296508/1‑149 (MQ=255)
                              gACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGGATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCTTTAACGCTGgca                                                                                                                     <  4:374877/149‑1 (MQ=255)
                              gACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                                                                                                                       >  8:282958/1‑147 (MQ=255)
                              gACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgca                                                                                                                     <  3:278186/149‑1 (MQ=255)
                               aCTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgcag                                                                                                                    >  8:411967/1‑149 (MQ=255)
                                  cccACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGATGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGAt                                                                                                                                                <  2:423178/118‑1 (MQ=255)
                                  cccACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGAt                                                                                                                                                >  1:423178/1‑118 (MQ=255)
                                   ccACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAt                                                                                                                >  4:132420/1‑149 (MQ=255)
                                   ccACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGGACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAt                                                                                                                >  8:326934/1‑149 (MQ=255)
                                          cATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGc                                                                                                         >  8:40337/1‑149 (MQ=255)
                                              cGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGta                                                                                                      <  3:132420/148‑1 (MQ=255)
                                              cGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGta                                                                                                      <  7:282958/148‑1 (MQ=255)
                                               ggCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGAGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtact                                                                                                    <  4:156273/149‑1 (MQ=255)
                                               ggCCTGGCTCGTCATACCGATAACGACGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGAGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtact                                                                                                    <  4:162244/149‑1 (MQ=255)
                                                       tCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactctact                                                                                               <  2:321311/146‑1 (MQ=255)
                                                       tCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactctact                                                                                               >  1:321311/1‑146 (MQ=255)
                                                              aCCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCg                                                                                      >  3:103820/1‑148 (MQ=255)
                                                               ccGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGGCGTCgca                                                                                    <  6:276082/149‑1 (MQ=255)
                                                                     aCGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCAcc                                                                               <  5:105402/148‑1 (MQ=255)
                                                                        ccGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGtt                                                                            >  7:21898/1‑148 (MQ=255)
                                                                                 tgtgGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTg                                                                  >  4:5899/1‑149 (MQ=255)
                                                                                        cGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTc                                                           <  7:40337/149‑1 (MQ=255)
                                                                                        cGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTc                                                           >  4:317887/1‑149 (MQ=255)
                                                                                         gggCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCt                                                          <  1:218761/149‑1 (MQ=255)
                                                                                           gCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGa                                                        <  7:411967/149‑1 (MQ=255)
                                                                                               cGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATgcg                                                    >  6:39597/1‑149 (MQ=255)
                                                                                                   cATTTTGTACGGTGATGATGATTGGCGCAAGGAGTATTGCTTCGCAATGGGACTCCGGTGCCAATGCATTAACGCTGGCAGCAAATAGCTGCGTACTCTACTCCGCCGTCGCATCACCGTTTAAGTCGTTGTCACTTCTGATGCGCAc                                                 <  4:76783/148‑1 (MQ=255)
                                                                                                    aTTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGc                                               >  7:343506/1‑149 (MQ=255)
                                                                                                           aCGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGgtgt                                         >  4:293602/1‑148 (MQ=255)
                                                                                                             gTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGtt                                                                            >  2:422532/1‑111 (MQ=255)
                                                                                                             gTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGtt                                                                            <  1:422532/111‑1 (MQ=255)
                                                                                                             gTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTAc                                      >  8:397671/1‑149 (MQ=255)
                                                                                                             gTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTAc                                      >  3:293618/1‑149 (MQ=255)
                                                                                                              ttGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGc                                                                                                         >  6:164999/1‑81 (MQ=255)
                                                                                                              ttGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGc                                                                                                         <  5:164999/81‑1 (MQ=255)
                                                                                                                 atgatgatTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTc                                  >  1:171991/1‑149 (MQ=255)
                                                                                                                   gatgatTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGc                                >  1:364814/1‑149 (MQ=255)
                                                                                                                   gatgatTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTATCGc                                >  7:341748/1‑149 (MQ=255)
                                                                                                                    atgatTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCTTTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGc                                <  8:141426/148‑1 (MQ=255)
                                                                                                                      gatTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATTCGCACGCTGGTGTTACTTTCGCTa                              <  4:296508/148‑1 (MQ=255)
                                                                                                                      gatTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTAt                             >  1:66862/1‑149 (MQ=255)
                                                                                                                        ttGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTc                           <  7:326934/149‑1 (MQ=255)
                                                                                                                              cATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCttt                     <  3:317887/149‑1 (MQ=255)
                                                                                                                               aTGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTtg                    <  8:21898/149‑1 (MQ=255)
                                                                                                                                tGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTtgt                   <  3:293602/149‑1 (MQ=255)
                                                                                                                                  gAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTgg                 <  2:66862/149‑1 (MQ=255)
                                                                                                                                     tATTGCTTTGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTc               >  6:26818/1‑148 (MQ=255)
                                                                                                                                     tATTGCTTTGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTc               >  6:32751/1‑148 (MQ=255)
                                                                                                                                     tATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTc               >  4:188394/1‑148 (MQ=255)
                                                                                                                                     tATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTc               >  1:232519/1‑148 (MQ=255)
                                                                                                                                     tATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTc               >  3:362385/1‑148 (MQ=255)
                                                                                                                                       ttGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCaaa            >  1:111977/1‑149 (MQ=255)
                                                                                                                                         gCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAAtt          >  6:226220/1‑149 (MQ=255)
                                                                                                                                          cTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAAttt         >  6:337709/1‑149 (MQ=255)
                                                                                                                                           ttCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTTGTCAAATTTg        >  2:374637/1‑149 (MQ=255)
                                                                                                                                           ttCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTg        >  8:188904/1‑149 (MQ=255)
                                                                                                                                             cGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTg        <  2:171991/147‑1 (MQ=255)
                                                                                                                                              gCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTGGt      <  2:364814/148‑1 (MQ=255)
                                                                                                                                              gCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGGGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTGGTc     <  7:397671/149‑1 (MQ=255)
                                                                                                                                                 aTGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTGGTCTGa  <  8:225575/149‑1 (MQ=255)
                                                                                                                                                 aTGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTGGTCTGa  >  7:225575/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTCTGCAAAGCGGTGTAGATGATCATACCAGACTCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTTAAGTCGTTGTCACTTCTGATGCGCACGCTGGTGTTACTTTCGCTATTCAGCTTTGTGGTCAAATTTGGTCTGA  >  NC_000913/1723118‑1723411

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