Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,008,422 G→T 6.0% intergenic (+45/+17) rhtC → / ← rhtB threonine efflux pump/homoserine, homoserine lactone and S‑methyl‑methionine efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4220GT 6.0% 332.0 / 7.8 115intergenic (+45/+17)rhtC/rhtBthreonine efflux pump/homoserine, homoserine lactone and S‑methyl‑methionine efflux pump
Reads supporting (aligned to +/- strand):  ref base G (52/56);  new base T (4/3);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

Only 100 of 119 total aligned reads displayed.

TGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTAACCG  >  NC_000913/4008274‑4008571
                                                                                                                                                    |                                                                                                                                                     
 gCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCa                                                                                                                                                     >  1:373000/1‑150 (MQ=255)
 gCCGCAAATGCGCCGTGGTTATCAACGGCTGGCGAAGTGGATTGCTGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGGCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCa                                                                                                                                                     <  1:815768/150‑1 (MQ=255)
    gCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                  <  1:1128233/150‑1 (MQ=255)
    gCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                                  <  2:373000/150‑1 (MQ=255)
     cAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                                 >  1:569688/1‑150 (MQ=255)
      aaaTGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                                >  1:1366292/1‑150 (MQ=255)
       aaTGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGAc                                                                                                                                               <  2:965058/150‑1 (MQ=255)
           cgcCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACatta                                                                                                                                           >  2:199291/1‑150 (MQ=255)
           cgcCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACatta                                                                                                                                           >  1:1358656/1‑150 (MQ=255)
           cgcCGTGGTTATCAACGTCTGGCGAAGTGGATTGAAGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGAGTTTACGCCGCATCCGACatta                                                                                                                                           >  2:774892/1‑150 (MQ=255)
            gcCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACattat                                                                                                                                          <  2:250403/150‑1 (MQ=255)
               gTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAtttt                                                                                                                                       <  1:1277670/150‑1 (MQ=255)
                    tATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGCAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACg                                                                                                                                  <  1:1098695/150‑1 (MQ=255)
                      tCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCa                                                                                                                                <  2:1242403/150‑1 (MQ=255)
                           gTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCt                                                                                                                           >  1:1342798/1‑150 (MQ=255)
                             cTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCg                                                                                                                         <  2:367741/150‑1 (MQ=255)
                               ggCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGcc                                                                                                                       >  2:617584/1‑150 (MQ=255)
                               ggCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCCTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGcc                                                                                                                       >  2:88387/1‑150 (MQ=255)
                                gCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCg                                                                                                                      >  2:1511445/1‑150 (MQ=255)
                                  gAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGAt                                                                                                                    >  2:1337481/1‑150 (MQ=255)
                                       ggATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTaa                                                                                                               >  2:1120510/1‑150 (MQ=255)
                                        gATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAAc                                                                                                              <  1:1371074/150‑1 (MQ=255)
                                          ttGATGGTTTTGCCGTTGCGTTATTTTCCGGATTTGGCATTGATTTGATTTTTTCGCGGTGATGCCAGACGCGTCTTCTGAGTAAGTCGGATAAGGGGTTTACGCCGCATACGACATTATTTTTCTCGCATGCCTCGCCGATGCTAACAg                                                                                                            <  2:1096683/150‑1 (MQ=255)
                                           tGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgc                                                                                                           <  1:509934/150‑1 (MQ=255)
                                            gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTAGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGTCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                          >  1:217457/1‑150 (MQ=255)
                                              tGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCt                                                                                                        <  1:18228/150‑1 (MQ=255)
                                               ggTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTc                                                                                                       <  2:935845/150‑1 (MQ=255)
                                                  tttGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCcca                                                                                                    >  2:121366/1‑150 (MQ=255)
                                                    tGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGGCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCccacc                                                                                                  <  2:764821/150‑1 (MQ=255)
                                                       cGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGc                                                                                               <  1:160538/150‑1 (MQ=255)
                                                       cGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGc                                                                                               >  1:1395872/1‑150 (MQ=255)
                                                          ggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATa                                                                                            <  1:810739/150‑1 (MQ=255)
                                                          ggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGAATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATa                                                                                            <  1:1510601/150‑1 (MQ=255)
                                                           gCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATaa                                                                                           >  1:998745/1‑150 (MQ=255)
                                                           gCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATaa                                                                                           >  1:540979/1‑150 (MQ=255)
                                                             gTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaac                                                                                         <  1:1329505/150‑1 (MQ=255)
                                                             gTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaac                                                                                         >  1:493634/1‑150 (MQ=255)
                                                              ttATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaaca                                                                                        <  2:1358656/150‑1 (MQ=255)
                                                               tATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaacaa                                                                                       <  1:418908/150‑1 (MQ=255)
                                                                aTTTTCCGGATTGTGCATTCATTTGTTTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGTATACGGCGTTTACGCCGCATCCGACATTATTTTTCAAGCATGCCTCGCCGATGCTAACAGCGCTCCCACAAGCATAaacaac                                                                                      <  2:1617300/150‑1 (MQ=255)
                                                                aTTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaacaac                                                                                      <  2:541141/150‑1 (MQ=255)
                                                                    gCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                                  <  2:1260416/150‑1 (MQ=255)
                                                                    gCCGGATTTGGCATTCATTTGATTATTTCGCGGGGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                                  <  2:552824/150‑1 (MQ=255)
                                                                     ccGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGcc                                                                                 >  1:827212/1‑150 (MQ=255)
                                                                      cGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCg                                                                                >  1:368817/1‑150 (MQ=255)
                                                                           ttGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAAt                                                                           >  2:818278/1‑150 (MQ=255)
                                                                              gCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCtt                                                                        >  1:1219959/1‑150 (MQ=255)
                                                                              gCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCtt                                                                        >  1:535185/1‑150 (MQ=255)
                                                                               cattcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGAATCCGACATTATTTTTCACGCTTGCCTCGCCGCTGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                                       >  1:829720/1‑150 (MQ=255)
                                                                                  tcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCTGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCCCAAGCATAAACAACGAGCCGAAAATCTTAATc                                                                    >  1:533112/1‑150 (MQ=255)
                                                                                    ttttGATTCTTTCTGGGGTATGCCAGACGCGTCTACATAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGAGATTATTTTTCACGCATTCCGCGCCGCTGCTAACAGCGCGCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAg                                                                  <  1:1154383/149‑1 (MQ=255)
                                                                                    attTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAg                                                                  >  1:1065445/1‑150 (MQ=255)
                                                                                    attTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAg                                                                  <  1:140433/150‑1 (MQ=255)
                                                                                    attTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGACGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAg                                                                  <  2:893836/150‑1 (MQ=255)
                                                                                      ttGATTATTTCGCTGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAATGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAACCGAAAATCTTATTCAgcg                                                                <  1:944575/150‑1 (MQ=255)
                                                                                       tGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGACGAAAATCTTATTCAgcgc                                                               >  2:298700/1‑150 (MQ=255)
                                                                                        gATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCTCCATCATAAACAACGAGCCGAAAATCTTATTCAgcgac                                                              >  1:89117/1‑148 (MQ=255)
                                                                                         attattTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCt                                                             <  2:569688/150‑1 (MQ=255)
                                                                                           tattTCGCGGTGATGCCAGACGAGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAacgcatcc                                                           <  1:199291/150‑5 (MQ=255)
                                                                                               tCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCt                                                       <  1:1511445/150‑1 (MQ=255)
                                                                                                     ggATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTg                                                 <  2:103647/149‑1 (MQ=255)
                                                                                                       aTGCCAGACGCTTCTTCAGAGTAAGTAGGATAAGGCGTTTACGCCGCAGCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGt                                               <  2:1054783/150‑1 (MQ=255)
                                                                                                       aTGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGt                                               >  1:291175/1‑150 (MQ=255)
                                                                                                         gCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTcc                                             >  2:1602882/1‑150 (MQ=255)
                                                                                                            aGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCttt                                          >  1:832867/1‑150 (MQ=255)
                                                                                                            aGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCttt                                          <  2:368817/150‑1 (MQ=255)
                                                                                                             gACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTa                                         <  1:1648694/150‑1 (MQ=255)
                                                                                                              aCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTaa                                        <  2:1366292/150‑1 (MQ=255)
                                                                                                                gcgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATc                                      <  1:774892/150‑1 (MQ=255)
                                                                                                                 cgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTAATACAGCGCCTTCATCTGCTTTGGTCCTATAATcc                                     >  1:339208/1‑150 (MQ=255)
                                                                                                                   tCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCGGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGACTTCATCTGCTTTGGTCCTTTAATCCAt                                   <  2:1638561/150‑1 (MQ=255)
                                                                                                                   tCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCAt                                   >  2:820829/1‑150 (MQ=255)
                                                                                                                   gcTTCAGAGTAAGTCGGATACGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCAt                                   <  2:1631988/149‑1 (MQ=255)
                                                                                                                    cTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGACTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCGTCATCTGCTTTGGTCGTTTAATCCATa                                  >  2:300077/1‑150 (MQ=255)
                                                                                                                     ttCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATag                                 <  1:1602882/150‑1 (MQ=255)
                                                                                                                      tCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATaga                                >  2:150485/1‑150 (MQ=255)
                                                                                                                        agagTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGc                              <  1:1672764/150‑1 (MQ=255)
                                                                                                                         gagTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCa                             >  1:484042/1‑150 (MQ=255)
                                                                                                                         gagTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCa                             >  1:906717/1‑150 (MQ=255)
                                                                                                                         gagTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCa                             <  1:12177/150‑1 (MQ=255)
                                                                                                                         aataatgTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCa                             >  2:177037/7‑150 (MQ=38)
                                                                                                                             aaGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATcc                         >  1:158674/1‑150 (MQ=255)
                                                                                                                               gTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt                       >  1:914135/1‑150 (MQ=255)
                                                                                                                               gTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt                       >  1:460503/1‑150 (MQ=255)
                                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGtt                      >  1:295024/1‑150 (MQ=255)
                                                                                                                                 cGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTg                     >  2:1544680/1‑150 (MQ=255)
                                                                                                                                  ggATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGa                    <  1:1434063/150‑1 (MQ=255)
                                                                                                                                   gATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAg                   >  2:891512/1‑150 (MQ=255)
                                                                                                                                    aTAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGc                  <  2:1585154/150‑1 (MQ=255)
                                                                                                                                     tAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCa                 >  2:185260/1‑150 (MQ=255)
                                                                                                                                     tAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCa                 <  1:1544712/150‑1 (MQ=255)
                                                                                                                                       aGGCGTTTACGCCGCATCCGACATTATTTTTCGCGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAg               <  2:408875/150‑1 (MQ=255)
                                                                                                                                        ggCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAgg              <  1:1338733/150‑1 (MQ=255)
                                                                                                                                            tatACGCCGCATCCGACATTATTTTTCACGCATGCCACGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTgg          <  2:1485551/148‑1 (MQ=255)
                                                                                                                                              tACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCg        <  2:1407496/150‑1 (MQ=255)
                                                                                                                                              tACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCg        <  2:827212/150‑1 (MQ=255)
                                                                                                                                              tACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCg        >  2:1482621/1‑150 (MQ=255)
                                                                                                                                              tACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCg        <  1:617584/150‑1 (MQ=255)
                                                                                                                                                  ccgcATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAGCGAGCCGAAAATCTTATTGAGGGCCTTCATCTGCTTGGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGTCGGAAc    >  2:601131/1‑150 (MQ=255)
                                                                                                                                                    gcATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGCGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCATTCCGTTGAGCAAGGGTGGCGTAACCg  >  2:195180/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTAACCG  >  NC_000913/4008274‑4008571

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: