Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,416,967 T→G S6A (TCA→GCA)  yjfP → acyl CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,416,9670TG100.0% 300.4 / NA 74S6A (TCA→GCA) yjfPacyl CoA esterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (45/29);  total (45/29)

ACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTG  >  NC_000913/4416820‑4417111
                                                                                                                                                   |                                                                                                                                                
aCCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCAc                                                                                                                                               >  2:1232783/1‑151 (MQ=255)
  cTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCAcgc                                                                                                                                             >  2:320099/1‑151 (MQ=255)
      tACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGc                                                                                                                                         >  2:1186712/1‑151 (MQ=255)
      tACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGc                                                                                                                                         >  2:1284280/1‑151 (MQ=255)
      tACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGc                                                                                                                                         >  1:212578/1‑151 (MQ=255)
          cgATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGc                                                                                                                                     >  2:471304/1‑151 (MQ=255)
          cgATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGc                                                                                                                                     >  1:188610/1‑151 (MQ=255)
          cgATTTTCGAATAAAAACGTGACGAAAGCGTAACCCACATCTCAACATTAGGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTgg                                                                                                                                      <  1:1202507/150‑1 (MQ=255)
           gATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGc                                                                                                                                     >  1:748624/1‑150 (MQ=255)
            aTTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGc                                                                                                                                     >  1:1113172/1‑149 (MQ=255)
               ttCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGat                                                                                                                                 <  1:623970/150‑1 (MQ=255)
                    aTGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTccc                                                                                                                           >  2:1207305/1‑151 (MQ=255)
                     tGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCg                                                                                                                          >  2:381184/1‑151 (MQ=255)
                       aaacCGTGACGAAACCTTACCCCACCTCCAAAAATTATGATCAAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGAGTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGtt                                                                                                                        <  1:579029/151‑1 (MQ=255)
                       aaaaCGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGtt                                                                                                                        >  1:553483/1‑151 (MQ=255)
                           cGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGttcttc                                                                                                                    >  2:878019/1‑151 (MQ=255)
                              gACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCa                                                                                                                   >  2:39168/1‑149 (MQ=255)
                              gACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATg                                                                                                                 >  1:1096307/1‑151 (MQ=255)
                               aCGAAAGCGTAACCCACATCTAAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCa                                                                                                                   <  1:815056/148‑1 (MQ=255)
                                cGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCt                                                                                                               >  1:352142/1‑151 (MQ=255)
                                  aacgggaagcccACAGCTCAACATTAGGAACTATAACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATg                                                                                                                 <  1:413203/138‑1 (MQ=255)
                                      cGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCt                                                                                                         >  1:249974/1‑151 (MQ=255)
                                        tAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGt                                                                                                       >  2:1000765/1‑151 (MQ=255)
                                              acaTCTAAACATTATAAACAAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCgg                                                                                                    <  1:945984/148‑1 (MQ=255)
                                               caTCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCa                                                                                                 >  2:178629/1‑150 (MQ=255)
                                                     aaCATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCaaaaa                                                                                             <  2:1262931/148‑1 (MQ=255)
                                                     aaCATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCaaa                                                                                               <  1:744991/146‑1 (MQ=255)
                                                        aTTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATAcc                                                                                       >  1:879343/1‑151 (MQ=255)
                                                         ttATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATAccc                                                                                      >  2:274786/1‑151 (MQ=255)
                                                          tATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATAcccc                                                                                     >  2:545951/1‑151 (MQ=255)
                                                           ataaaCTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATAcccc                                                                                     <  1:845696/147‑1 (MQ=255)
                                                           aTGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCg                                                                                    >  1:482761/1‑151 (MQ=255)
                                                              aaCTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGt                                                                                   <  1:986303/149‑1 (MQ=255)
                                                                  aaTACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGt                                                                             <  1:877178/151‑1 (MQ=255)
                                                                   aTACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTg                                                                            >  1:1089663/1‑151 (MQ=255)
                                                                        gCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTa                                                                        >  1:175252/1‑150 (MQ=255)
                                                                          cagcagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAAtt                                                                     >  1:402814/1‑151 (MQ=255)
                                                                          cagcagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAAtt                                                                     >  2:343719/1‑151 (MQ=255)
                                                                           agcagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAAttt                                                                    >  2:321212/1‑151 (MQ=255)
                                                                           agcagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAAttt                                                                    <  1:691000/151‑1 (MQ=255)
                                                                            gcagAAACGTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACACTAACTGGCAGATATTCCCGTACTGCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAAt                                                                      <  1:917431/148‑1 (MQ=255)
                                                                                aaaCTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTAt                                                               >  1:173644/1‑151 (MQ=255)
                                                                                 aaCTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATc                                                              >  2:860156/1‑151 (MQ=255)
                                                                                    ttGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACACATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACg                                                           >  1:849053/1‑151 (MQ=255)
                                                                                        gCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCt                                                        <  1:1058986/150‑1 (MQ=255)
                                                                                           aaGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCttt                                                      <  1:569854/149‑1 (MQ=255)
                                                                                             ggATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACtt                                                  >  2:539626/1‑151 (MQ=255)
                                                                                             ggATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACtt                                                  <  2:1267959/151‑1 (MQ=255)
                                                                                                tttCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCa                                                <  1:800730/150‑1 (MQ=255)
                                                                                                 ttcttcATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCAt                                               <  2:212113/150‑1 (MQ=255)
                                                                                                  tcttcATGTCATACCAGACAAACAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAGAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATcc                                             >  1:86834/1‑151 (MQ=255)
                                                                                                       aTGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCt                                        <  1:321206/151‑1 (MQ=255)
                                                                                                       aTGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCt                                        >  2:160952/1‑151 (MQ=255)
                                                                                                       aTGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCt                                        <  1:14190/151‑1 (MQ=255)
                                                                                                         gTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTgg                                      >  1:320014/1‑151 (MQ=255)
                                                                                                           cATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGgtg                                    <  1:211066/151‑1 (MQ=255)
                                                                                                           cATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGgtg                                    >  2:617682/1‑151 (MQ=255)
                                                                                                            aTACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGgtgt                                   >  2:1162877/1‑151 (MQ=255)
                                                                                                              aCCAGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGtat                                 >  2:65298/1‑151 (MQ=255)
                                                                                                               ccAGACAAAAAGGAGAGGCAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGta                                  <  1:842055/149‑1 (MQ=255)
                                                                                                                 aGACAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGc                              <  2:1226931/151‑1 (MQ=255)
                                                                                                                   aCAAAAAGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTa                            <  1:177732/151‑1 (MQ=255)
                                                                                                                     aaaaaGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTAtt                          >  2:132476/1‑151 (MQ=255)
                                                                                                                       aaaGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTAtttt                        <  2:1184647/151‑1 (MQ=255)
                                                                                                                        aaGGAGAGACAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTg                       >  1:947949/1‑151 (MQ=255)
                                                                                                                           gagagaCAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCg                    <  1:1159801/151‑1 (MQ=255)
                                                                                                                              agaCAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTg                 <  1:1185982/151‑1 (MQ=255)
                                                                                                                               gaCAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTgc                >  1:1099464/1‑151 (MQ=255)
                                                                                                                                aCAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTgcg               <  2:1129899/151‑1 (MQ=255)
                                                                                                                                 cAGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTgcgc              >  1:1109220/1‑151 (MQ=255)
                                                                                                                                  aGATGATTGAAATAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCt             <  1:380700/151‑1 (MQ=255)
                                                                                                                                           aaaTAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGc    >  2:942132/1‑151 (MQ=255)
                                                                                                                                             aTAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTg  >  1:927630/1‑151 (MQ=255)
                                                                                                                                             aTAGAAGCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTg  >  2:902384/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTG  >  NC_000913/4416820‑4417111

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