Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,008,422 G→T 5.7% 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 5.7% 568.6 / 14.0 156intergenic (+45/+17)rhtC/rhtBthreonine efflux pump/homoserine, homoserine lactone and S‑methyl‑methionine efflux pump
Reads supporting (aligned to +/- strand):  ref base G (57/90);  new base T (4/5);  total (61/95)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 160 total aligned reads displayed.

TGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTAA  >  NC_000913/4008274‑4008568
                                                                                                                                                    |                                                                                                                                                  
tGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAcgccgc                                                                                                                                                   >  2:1268254/1‑150 (MQ=255)
      aaaTGCGCCGTGGTTATCAACGTCTGGCGTAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                             <  1:1631203/150‑1 (MQ=255)
      aaaTGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                                             <  1:697827/150‑1 (MQ=255)
      aaaTGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                                             <  1:1239450/150‑1 (MQ=255)
       aaTGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGAc                                                                                                                                            <  2:1832360/150‑1 (MQ=255)
           cgcCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACatta                                                                                                                                        <  2:1964245/150‑1 (MQ=255)
             ccGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACattatt                                                                                                                                      <  2:1478331/150‑1 (MQ=255)
             ccGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACattatt                                                                                                                                      <  2:790526/150‑1 (MQ=255)
             ccGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTGTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACattatt                                                                                                                                      <  2:2397085/150‑1 (MQ=255)
                tGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGGATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAttttt                                                                                                                                   >  1:166515/1‑150 (MQ=255)
                tGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAttttt                                                                                                                                   <  2:1199534/150‑1 (MQ=255)
                gggTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAttttt                                                                                                                                   <  2:1199844/149‑1 (MQ=255)
                  gTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCa                                                                                                                                 >  2:1727908/1‑150 (MQ=255)
                  gTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCa                                                                                                                                 >  2:1859106/1‑150 (MQ=255)
                     aTCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGc                                                                                                                              <  2:1923746/150‑1 (MQ=255)
                      tCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCa                                                                                                                             >  2:1659124/1‑150 (MQ=255)
                      tCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCa                                                                                                                             >  2:907210/1‑150 (MQ=255)
                        aaCGTCTGGCGAAGTGGATTGATGGTTTTGCGTGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATg                                                                                                                           <  1:2375649/150‑1 (MQ=255)
                        aaCGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATg                                                                                                                           <  1:817604/150‑1 (MQ=255)
                        aaCGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATg                                                                                                                           <  1:2375408/150‑1 (MQ=255)
                        aaCGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATg                                                                                                                           <  1:2507/150‑1 (MQ=255)
                        aaCGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATg                                                                                                                           <  1:1346859/150‑1 (MQ=255)
                        aaCGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATg                                                                                                                           <  1:1508502/150‑1 (MQ=255)
                          cGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCGCGCATGcc                                                                                                                         <  1:396929/150‑1 (MQ=255)
                          cGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGcc                                                                                                                         <  1:115411/150‑1 (MQ=255)
                          cGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGcc                                                                                                                         <  1:1359658/150‑1 (MQ=255)
                          cGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGcc                                                                                                                         <  2:2161657/150‑1 (MQ=255)
                           gTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCCGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCAGCCGACATTATTTTTCACGCATGCCt                                                                                                                        >  1:1149401/1‑150 (MQ=255)
                           gTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCt                                                                                                                        >  1:141390/1‑150 (MQ=255)
                            tCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTc                                                                                                                       >  2:1214765/1‑150 (MQ=255)
                              tGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGc                                                                                                                     <  2:1583453/150‑1 (MQ=255)
                                  gAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCGGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGAt                                                                                                                 <  2:937948/150‑1 (MQ=255)
                                  gAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGAt                                                                                                                 <  2:937933/150‑1 (MQ=255)
                                   aaGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATg                                                                                                                <  1:124162/150‑1 (MQ=255)
                                         aTTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACa                                                                                                          <  1:1876066/150‑1 (MQ=255)
                                         aTTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACa                                                                                                          <  1:1876081/150‑1 (MQ=255)
                                          ttGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAg                                                                                                         <  1:387303/150‑1 (MQ=255)
                                          ttGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAg                                                                                                         >  2:810653/1‑150 (MQ=255)
                                            gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                       <  1:2289070/150‑1 (MQ=255)
                                            gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                       <  1:783045/150‑1 (MQ=255)
                                            gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                       <  1:837967/150‑1 (MQ=255)
                                            gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACACGTCTTCAGAGTAAGTCGGATAAGGAGTATACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                       >  2:63952/1‑150 (MQ=255)
                                            gATGGTTGTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                       <  1:1522015/150‑1 (MQ=255)
                                             aTGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTTTTTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAAGAAGCATCCGACATTATTTTTAAAGCATGCCTCGCCGATGCTAACAgagc                                                                                                      >  2:2231283/1‑150 (MQ=255)
                                                 ttttGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTccc                                                                                                  <  2:469419/150‑1 (MQ=255)
                                                 ttttGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTccc                                                                                                  >  2:16019/1‑150 (MQ=255)
                                                 ttttGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTccc                                                                                                  >  2:619106/1‑150 (MQ=255)
                                                        ggggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCa                                                                                           <  2:1634878/150‑1 (MQ=255)
                                                        ggggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCa                                                                                           <  1:1225461/150‑1 (MQ=255)
                                                        ggggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCa                                                                                           <  1:660234/150‑1 (MQ=255)
                                                        ggggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCa                                                                                           <  2:1227640/150‑1 (MQ=255)
                                                        ggggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCa                                                                                           <  2:2000029/150‑1 (MQ=255)
                                                                    gCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                               <  2:2300857/150‑1 (MQ=255)
                                                                    gCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                               <  2:192887/150‑1 (MQ=255)
                                                                          tttGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGaaaa                                                                         <  2:952244/150‑1 (MQ=255)
                                                                          tttGGCATTCATTAGATTATTTCGCGGTGATGCCAGACGCGACTTCAGAGTAAGTCGGATAAGGCGATTACGCCGCATCCGACAATAATTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGaaaa                                                                         <  2:2148359/150‑1 (MQ=255)
                                                                               cattcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGATTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCATCATAAACAACGAGCCGAAAATCTTa                                                                    <  1:308407/150‑1 (MQ=255)
                                                                               cattcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                                    >  2:418186/1‑150 (MQ=255)
                                                                               cattcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                                    <  1:308329/150‑1 (MQ=255)
                                                                                     tttGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAgc                                                              <  2:706902/150‑1 (MQ=255)
                                                                                     tttGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAgc                                                              <  2:1807205/150‑1 (MQ=255)
                                                                                         attattTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCt                                                          >  2:2116971/1‑150 (MQ=255)
                                                                                            attTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCa                                                       >  2:661484/1‑150 (MQ=255)
                                                                                             tttCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCAt                                                      >  2:249180/1‑150 (MQ=255)
                                                                                               tCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCt                                                    >  1:885101/1‑150 (MQ=255)
                                                                                                 gcgGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGc                                                  >  1:1233871/1‑150 (MQ=255)
                                                                                                 gcgGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGc                                                  >  1:5996/1‑150 (MQ=255)
                                                                                                        tGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTc                                           <  2:103007/150‑1 (MQ=255)
                                                                                                          ccAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCt                                         >  1:1058337/1‑150 (MQ=255)
                                                                                                               cgcgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAAt                                    <  1:1069548/150‑1 (MQ=255)
                                                                                                               cgcgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAAt                                    <  2:252375/150‑1 (MQ=255)
                                                                                                               cgagTCTTCAGAGTAAGTCGGATAAGGCTTTTACGCCGCATCAGAAATTATTTTTCACGCATGCCTCGCCTTTGCTAACAGCGCTACCACAAGCATAAACAACGATCCGTAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAAt                                    <  2:1129939/147‑1 (MQ=255)
                                                                                                                 cgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATcc                                  <  1:475151/150‑1 (MQ=255)
                                                                                                                 cgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATcc                                  <  1:2063823/150‑1 (MQ=255)
                                                                                                                 cgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATcc                                  <  1:1637786/150‑1 (MQ=255)
                                                                                                                   tCTTCGGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCAt                                <  1:1520704/150‑1 (MQ=255)
                                                                                                                   tCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCAt                                <  2:1149401/150‑1 (MQ=255)
                                                                                                                    cTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATa                               <  1:1823335/150‑1 (MQ=255)
                                                                                                                     ttgagagTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATag                              >  1:725801/4‑150 (MQ=255)
                                                                                                                     ttCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATag                              >  1:1200410/1‑150 (MQ=255)
                                                                                                                     gtCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATag                              <  2:856411/149‑1 (MQ=255)
                                                                                                                       cAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATagag                            >  2:745991/1‑150 (MQ=255)
                                                                                                                         gagTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCtt                                        <  1:73234/136‑1 (MQ=255)
                                                                                                                             atgTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATcc                      >  1:1402138/3‑150 (MQ=255)
                                                                                                                             atgTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATcc                      >  2:925551/3‑150 (MQ=255)
                                                                                                                             atgTCGGATGCGGCGTAAACGCCGTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTAAGCGCCTTCATCTGCTTTGGTCATTTAAACCATAGAGCAATcc                      >  1:1825347/3‑150 (MQ=255)
                                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTTACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGtt                   <  2:924920/150‑1 (MQ=255)
                                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGtt                   <  2:1556212/150‑1 (MQ=255)
                                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGtt                   <  2:1556210/150‑1 (MQ=255)
                                                                                                                                tCGGATAAGGCGTGGACGCCGCAGCCGACAGTATTTGACACGCCTGCCTCGCCGAGGCTAACCGCGCTCCCACCAGCATAAACAACGTGCCGAAAATCTTCGTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGtt                   <  2:1759662/150‑1 (MQ=255)
                                                                                                                                  ggATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGa                 <  1:1380453/150‑1 (MQ=255)
                                                                                                                                   gATAAGGCGTTTACGCCGCATCCGACATTTTTTTTCACGCTTGCCTCGGCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAg                <  2:2226606/150‑1 (MQ=255)
                                                                                                                                   gATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAg                >  2:712437/1‑150 (MQ=255)
                                                                                                                                   gATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAg                >  2:1680726/1‑150 (MQ=255)
                                                                                                                                    aTAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGc               <  2:1768629/150‑1 (MQ=255)
                                                                                                                                       aGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAg            <  1:671832/150‑1 (MQ=255)
                                                                                                                                       aGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAg            <  1:992962/150‑1 (MQ=255)
                                                                                                                                        ggCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAgg           >  1:540189/1‑150 (MQ=255)
                                                                                                                                                 gccgcATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTaa  <  1:2274568/150‑1 (MQ=255)
                                                                                                                                                 gccgcATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTaa  <  2:1514774/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTAA  >  NC_000913/4008274‑4008568

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: